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I.  INTRODUCTION 


Digital  breast  tomosynthesis  (DBT)  is  a  3D  imaging  modality  in  which  tomographic  sections  of  the 
breast  are  generated  from  a  limited  range  of  x-ray  projections.  Preliminary  studies  have  indicated  that  DBT  has 
improved  sensitivity  and  specificity  for  early  cancer  detection  relative  to  conventional  2D  digital 
mammography  (DM).1,  2  In  my  dissertation  research  in  Dr.  Andrew  Maidment’s  lab  at  the  University  of 
Pennsylvania,  I  have  worked  to  develop  analytical  models  of  the  transfer  functions  of  DBT  for  assessment  of 
image  quality.  Although  the  transfer  functions  of  DBT  have  been  modeled  by  other  researchers,  their  approach 
is  limited  by  approximating  the  angle  of  incidence  as  normal  to  the  detector  at  all  points  in  the  projection 
images.  In  modeling  the  anisotropy  of  the  incident  angle  at  each  point  of  the  detector,  I  have  developed 
rigorous  models  of  image  quality  in  DBT  that  would  not  be  possible  under  the  formulation  of  previous  authors. 

II.  BODY 

II.  A.  Resolution  Loss  Due  to  Oblique  X-Ray  Incidence 

In  the  winter  and  early  spring  of  2011,  much  of  my  research  effort  was  devoted  to  writing  a  peer- 
reviewed  publication  analyzing  the  effect  of  oblique  x-ray  incidence  on  the  resolution  and  noise  properties  of 
projection  images  in  DBT.  It  ultimately  felt  rewarding  to  see  this  work  published  in  the  November  2011  issue 
of  the  journal  Medical  Physics .4 

In  my  Medical  Physics  publication,  I  have  extended  Swank’s  calculations  of  the  transfer  functions  of 
turbid  granular  phosphors5  to  oblique  x-ray  incidence  using  the  diffusion  equation6  to  model  the  spread  of 
visible  light  in  a  scintillator.  My  work  successfully  demonstrated  degradation  in  the  modulation  transfer 
function  (MTF)  and  detective  quantum  efficiency  (DQE)  with  oblique  x-ray  incidence.  The  model  is  set  apart 
from  prior  studies  '  on  oblique  x-ray  incidence  in  being  derived  from  first  principles  and  in  closed  form 
without  making  the  simplifying  assumption  that  the  point  spread  function  (PSF)  of  normal  incidence  is  a  delta 
function.10 

In  the  first  specific  aim  (SA1)  of  my  predoctoral  training  grant  through  the  Department  of  Defense 
(DOD)  Breast  Cancer  Research  Program,  I  proposed  to  demonstrate  the  spatial  anisotropy  of  the  transfer 
functions  of  DBT  by  carefully  modeling  the  incident  angle  at  each  point  on  the  detector.  To  illustrate  the 
anisotropy  of  one  of  the  transfer  functions  in  an  individual  projection,  DQE  was  plotted  versus  position  on  the 
detector  at  a  fixed  spatial  frequency  in  my  Medical  Physics  publication  (Figure  7,  p.  6195).4  I  showed  that  an 
oblique  projection  has  greater  spatial  variation  in  DQE  than  the  central  projection. 

In  addition,  in  SA2  of  my  predoctoral  DOD  grant,  I  proposed  to  investigate  strategies  for  optimization  of 
image  quality  in  DBT.  One  important  result  published  in  my  Medical  Physics  paper  was  the  calculation  of  the 
projection  angle  dependence  of  the  optimal  phosphor  thickness  which  maximizes  DQE  (Figures  8-9,  p.  6196).4 
My  work  demonstrated  that  the  optimal  phosphor  thickness  is  projection  angle  dependent,  favoring  smaller 
thicknesses  with  higher  incident  angles. 

II.  B.  Super-Resolution  in  DBT  Image  Reconstructions 

Throughout  the  past  year,  I  have  also  pursued  a  project  investigating  the  existence  of  super-resolution  in 
DBT.  A  full  draft  of  a  paper  on  super-resolution  will  soon  be  sent  to  the  journal  Medical  Physics  for 
consideration  as  a  peer-reviewed  research  publication.11  This  work  expands  upon  a  proceedings  manuscript 
presented  in  the  2011  Physics  of  Medical  Imaging  Conference  of  the  Society  of  Photo-optical  Instrumentation 
Engineers  (SPIE).12  The  proceedings  manuscript,  which  was  supported  by  a  different  funding  mechanism  than 
my  predoctoral  DOD  grant,  earned  the  Best  Student  Paper  Award  at  the  conference  in  February  of  201 1  in  Lake 
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Buena  Vista,  Florida.  My  award  ultimately  received  coverage  in  an  article  written  in  March  of  2011  on  the 
website  medicalphysicsweb.org.13 

In  DBT,  the  image  of  an  object  is  translated  in  sub-pixel  detector  element  increments  with  each 
increasing  projection  angle.  To  my  knowledge,  my  work  is  the  first  to  demonstrate  the  potential  for  super¬ 
resolution  ( i.e .,  sub-pixel  resolution)  in  DBT  reconstructions  as  a  result  of  this  property.  In  order  to  develop  an 
analytical  model  of  super-resolution,  I  calculate  the  reconstruction  of  a  thin  input  object  whose  attenuation 
coefficient  varies  sinusoidally  with  position  within  a  rectangular  prism  parallel  to  the  breast  support.  This  input 
object  is  termed  a  “sine  plate”  in  my  manuscript  submission  for  the  journal  Medical  Physics. 11  To  investigate 
super-resolution,  the  frequency  of  the  sinusoidal  attenuation  coefficient  is  specified  to  be  higher  than  the  alias 
frequency  of  the  detector.  Using  analytical  modeling,  I  demonstrate  that  the  central  projection  represents  this 
high  frequency  input  as  if  it  were  a  lower  frequency,  and  that  its  Fourier  transform  is  maximized  at  a  lower 
frequency  than  the  input  as  evidence  of  aliasing.  By  contrast,  the  high  frequency  input  can  be  resolved  by 
performing  a  simple  backprojection  (SBP)  reconstruction  on  a  grid  whose  pixel  size  is  much  smaller  than  that 
of  the  detector  elements.  Applying  filters  to  the  reconstruction  smoothens  pixilation  artifacts  in  the  spatial 
domain  and  reduces  spectral  leakage  in  the  Fourier  domain. 

My  analytical  model  of  super-resolution  illustrates  the  importance  of  carefully  calculating  the  incident 
angle  at  each  point  on  the  detector  as  proposed  in  my  DOD  predoctoral  grant.  In  order  to  determine  the  sub¬ 
pixel  shifts  in  the  image  of  the  sine  plate  between  projections,  it  is  necessary  to  model  the  spatial  anisotropy  of 
the  incident  angle  on  the  detector.  Because  prior  authors  have  approximated  the  incident  angle  as  normal  to  the 
detector  in  all  projections,3  they  are  unable  to  derive  the  sub-pixel  shifts  in  the  image  of  the  object  between 
projections,  and  hence  are  not  able  to  model  the  presence  of  super-resolution  in  the  reconstruction. 

Using  a  bar  pattern  phantom,  the  existence  of  super-resolution  in  DBT  was  verified  experimentally. 
Individual  projections  of  the  bar  pattern  phantom  exhibited  classical  signs  of  aliasing,14  such  as  representing 
high  frequency  patterns  at  an  erroneous  orientation  with  less  line  pairs  than  expected.  By  contrast,  bar  patterns 
with  frequencies  higher  than  the  alias  frequency  of  the  detector  were  successfully  resolved  in  reconstructions 
performed  on  very  fine  grids.  In  select  clinical  examples,  super-resolution  was  found  to  improve  the  visibility 
of  microcalcifications,  which  are  subtle  signs  of  breast  cancer  in  many  women.15  By  reconstructing  on  grids 
whose  pixel  sizes  are  much  smaller  than  the  detector  elements,  the  morphology  of  the  calcifications  was  seen 
more  clearly,  and  small  calcifications  that  would  not  otherwise  be  visible  became  apparent  due  to  super¬ 
resolution.  For  these  reasons,  super-resolution  satisfies  a  central  goal  in  SA2  of  my  predoctoral  DOD  grant, 
which  was  to  develop  techniques  to  improve  image  quality. 

Finally,  my  research  has  demonstrated  spatial  anisotropies  in  the  existence  of  super-resolution 
depending  upon  the  directionality  of  the  input  frequency.  Although  super-resolution  is  achievable  over  a  broad 
range  of  positions  for  an  input  frequency  parallel  to  the  chest  wall  side  of  the  breast  support,  super-resolution  is 
feasible  at  fewer  positions  for  an  input  frequency  perpendicular  to  the  chest  wall.  Super-resolution  along  the 
latter  directionality  is  achievable  only  at  positions  sufficiently  displaced  from  the  plane  of  the  chest  wall  and 
from  the  mid  plane  perpendicular  to  the  breast  support  and  to  the  chest  wall.  The  anisotropy  of  super-resolution 
has  been  demonstrated  analytically  from  the  calculation  of  the  reconstruction  of  a  sine  plate,  as  well  as 
experimentally  from  the  reconstruction  of  bar  patterns.  This  finding  complements  the  central  aim  of  SA1  in  my 
predoctoral  DOD  grant,  which  was  to  demonstrate  spatial  anisotropy  in  the  image  quality  of  DBT  systems. 

II.  C.  Oblique  Reconstructions  in  DBT 

My  research  on  super-resolution  initially  assumed  a  reconstruction  plane  parallel  to  the  breast  support. 
In  a  subsequent  study,  I  demonstrated  that  super-resolution  is  also  achievable  in  an  obliquely  pitched 
reconstruction  plane.  This  work  was  accepted  for  oral  presentation  at  the  Annual  Meeting  of  the  American 
Association  of  Physicists  in  Medicine  (AAPM)  in  Vancouver,  British  Columbia  on  August  2,  2011. 16  In 
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addition,  similar  work  was  accepted  as  a  poster  presentation  at  the  Era  of  Hope  conference  in  Orlando,  Florida 
from  August  2,  2011  to  August  5,  2011.  Because  the  Era  of  Hope  conference  presented  a  conflict  with  the 
Annual  AAPM  Meeting,  my  graduate  supervisor,  Dr.  Maidment,  presented  this  topic  at  Vancouver  in  my 
absence. 

In  order  to  demonstrate  the  existence  of  super-resolution  in 
an  oblique  reconstruction  plane,  a  theoretical  framework  was 
developed  in  which  the  reconstruction  of  an  obliquely  pitched  sine 
plate  was  calculated  (Figure  1).  With  this  model,  I  have  studied 
differences  in  the  output  profiles  of  individual  projections  and 
reconstructions.  A  Selenia  Dimensions  DBT  system  (Hologic  Inc., 

Bedford  MA)  with  140  pm  detector  elements  and  15  projections 
acquired  over  a  15°  arc  was  simulated.  Also,  the  input  object  was 
oriented  at  a  pitch  (£)  of  30°  relative  to  the  breast  support.  To 
investigate  the  potential  for  super-resolution,  the  input  frequency 
(5.0  lp/mm)  was  specified  to  be  higher  than  the  detector  alias 
frequency,  3.6  lp/mm. 

Figure  2  illustrates  that  the  central  projection  cannot 
resolve  the  high  frequency  input  object.  As  evidence  of  aliasing, 
the  major  peak  in  its  Fourier  transform  occurs  at  a  lower 
frequency  (1.37  lp/mm)  than  the  input  frequency,  5.0  lp/mm.  By 
contrast,  SBP  reconstruction  can  resolve  the  input  and  its  major 
Fourier  peak  correctly  occurs  at  5.0  lp/mm.  Adding  ramp  (RA)  and  spectrum  apodization  (SA)  filters3  to  the 
reconstruction  smoothens  pixilation  artifacts  seen  with  SBP.  These  innovative  results  are  unexpected  based  on 
the  conventional  interpretation  of  the  Central  Slice  Theorem.3  According  to  that  theorem,  Fourier  space  should 
be  sampled  only  within  double-napped  cones  (DNCs)  whose  opening  angle  matches  the  angular  range  of  the 
DBT  scan.  The  analytical  modeling  illustrates  that  super-resolution  is  possible  well  outside  the  DNCs  whose 
opening  angle  spans  -7.5°  to  +7.5°  for  the  Selenia  Dimensions  system. 

The  feasibility  of  super-resolution  in  oblique  reconstructions  was  later  verified  experimentally  using  a 
lead  bar  pattern  phantom.  The  phantom  was  placed  on  a  goniometry  stand  at  a  height  of  7.6  cm  above  the 
breast  support  of  a  Selenia  Dimensions  DBT  system.  The  goniometer  was  adjusted  to  vary  the  pitch  of  the  bar 
patterns  during  the  acquisition  of  15  projections  spanning  a  15°  arc.  Reconstruction  was  performed  along  the 
obliquely  pitched  plane  of  the  bar  patterns  using  a  commercial  prototype  reconstruction  algorithm  (Briona™, 
Real  Time  Tomography,  Villanova,  PA).  In  the  plane  of  the  reconstruction,  the  pixel  size  of  the  reconstruction 
grid  (1 1.7  pm)  was  much  smaller  than  that  of  the  detector  elements  (140  pm),  so  that  the  alias  frequency  of  the 
reconstruction  grid  (42.7  lp/mm)  was  considerably  higher  than  the  alias  frequency  of  the  detector  (3.57  lp/mm). 

In  the  central  projection  acquired  with  bar  patterns  along  a  0°  pitch,  frequencies  up  to  3.0  lp/mm  are 
resolved  (Figure  3).  This  finding  is  expected,  since  3.0  lp/mm  is  below  the  detector  alias  frequency,  3.6  lp/mm 
for  140  pm  detector  elements.  At  4.0  lp/mm,  however,  a  broad  Moire  band14  is  present  perpendicular  to  the  bar 
patterns,  and  at  5.0  lp/mm,  the  bar  patterns  are  represented  as  if  they  were  a  lower  frequency  (~3.0  lp/mm).  By 
contrast,  the  reconstruction  clearly  shows  spatial  frequencies  up  to  5.0  lp/mm.  Ultimately,  the  existence  of 
super-resolution  using  the  30°  pitch  is  significant  because  the  DNCs  in  frequency  space  have  an  opening  angle 
spanning  approximately  -7.5°  to  +7.5°  for  measurements  made  near  the  midpoint  of  the  chest  wall.  At  these 
locations  along  the  detector,  the  divergence  of  the  x-ray  beam  from  the  nominal  projection  angles  is  minimal, 
and  the  30°  pitch  is  therefore  well  outside  the  DNCs  of  frequency  space. 


Figure  1.  A  schematic  diagram  of  the  DBT  acquisition 
geometry  is  shown  (figure  not  to  scale).  The  input  object 
(a  sine  plate)  has  an  attenuation  coefficient  which  varies 
sinusoidally  along  the  angular  pitch  f. 
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Figure  2.  For  a  thin  input  object  with  a  frequency  of  5.0  lp/mm  and  an  oblique  pitch  (0  of  30°,  the  central  projection,  SBP,  and  FBP  reconstructions 
are  compared  in  both  the  spatial  and  Fourier  domains.  The  central  height  of  the  input  object  above  the  detector  is  2.5  cm.  A  Selenia  Dimensions 
DBT  system  (Hologic  Inc.,  Bedford,  MA)  was  simulated  with  140  pm  detector  elements  and  15  projections  acquired  over  a  15°  arc. 


Central  Projection  (0°  Pitch)  Reconstruction  (0°  Pitch)  Reconstruction  (30°  Pitch) 


Figure  3.  A  bar  pattern  phantom  was  imaged  with  the  Selenia  Dimensions  x-ray  unit  on  a  goniometry  stand  at  0°  and  30°  pitches.  The 
reconstructions  along  both  pitches  can  resolve  higher  frequencies  than  a  single  projection,  providing  experimental  evidence  of  super-resolution. 
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To  demonstrate  the  feasibility  of  super-resolution  in  oblique 
reconstruction  planes  using  clinical  images,  reconstructions  of  microcalcifications 
were  performed  on  grids  with  the  pixel  size  28  pm  using  the  Selenia  Dimensions 
system.  As  shown  in  Figure  4,  the  visibility  of  microcalcifications  in  the  lower 
half  of  the  reconstruction  does  not  differ  considerably  using  either  a  0°  or  30° 
pitch.  Although  the  upper  right  calcifications  are  not  visible  with  the  30°  pitch, 
visibility  can  be  improved  by  simply  translating  the  plane  of  reconstruction  until 
they  come  into  focus.  Having  a  clear  image  of  the  structural  features  of 
microcalcifications  at  oblique  pitches  is  useful  because  there  may  be  cases  whose 
morphological  properties  are  best  assessed  along  these  directions. 

II.  D.  Optimization  of  Scan  Time  in  DBT 

In  SA2  of  my  predoctoral  DOD  grant,  I  proposed  to  investigate  differences  in  image  quality  comparing 
systems  with  continuous  tube  motion  and  step-and-shoot  motion.  Systems  with  continuous  tube  motion  have 
the  advantage  of  shorter  scan  time  and  hence  less  patient  motion  than  systems  with  step-and-shoot  motion. 
Their  drawback,  however,  is  blurring  due  to  focal  spot  motion.3, 18, 19 

In  this  work,  I  have  used  the  concept  of  a  sine  plate  to  investigate  subtleties  of  the  imaging  system  that 
were  not  modeled  in  my  super-resolution  research.  Unlike  my  super-resolution  research  which  assumed  a 
stationary  sine  plate,  this  project  simulates  the  motion  of  the  sine  plate  at  a  fixed  velocity  during  the  DBT  scan 
in  order  to  simulate  patient  motion.  The  reconstruction  of  the  sine  plate  is  calculated  for  various  input 
frequencies  and  for  either  continuous  tube  motion  or  step-and-shoot  motion.  My  work  determines  the  optimal 
scan  time  which  minimizes  the  trade-offs  between  patient  motion  and  focal  spot  motion  by  maximizing  the 
modulation  of  the  reconstruction  in  a  system  with  continuous  tube  motion. 

My  work  demonstrates  that  a  system  with  step-and-shoot  motion  will  always  have  superior  modulation 
relative  to  a  system  with  continuous  tube  motion  and  with  otherwise  similar  acquisition  parameters. 
Consequently,  in  order  to  optimize  a  step-and-shoot  system,  the  scan  time  yielding  the  same  modulation  as  the 
highest  achievable  with  continuous  tube  motion  is  calculated.  This  scan  time  provides  a  threshold  below  which 
the  benefits  of  step-and-shoot  motion  are  not  outweighed  by  its  drawbacks. 

This  analysis  of  image  quality  in  systems  with  patient  motion  and  either  continuous  tube  motion  or  step- 
and-shoot  motion  is  presented  in  an  upcoming  proceedings  manuscript  for  the  SPIE  Physics  of  Medical  Imaging 
conference  in  February,  2012  in  San  Diego,  California.20  Because  no  other  researcher  has  modeled  both  tube 
motion  and  patient  motion  simultaneously,  my  work  is  the  first  to  demonstrate  a  platform  for  optimizing  scan 
time  in  DBT. 

III.  KEY  RESEARCH  ACCOMPLISHMENTS 

•  In  a  recently  published  article  for  the  journal  Medical  Physics ,4  I  have  proposed  the  first  closed  form  model 
of  the  effect  of  oblique  x-ray  incidence  on  the  transfer  functions  of  DBT  projections  without  making  the 
simplifying  assumption  that  the  PSF  for  normal  incidence  is  a  delta  function.  Following  an  aim  established 
in  SA1  of  my  predoctoral  DOD  grant,  this  model  has  allowed  me  to  demonstrate  the  spatial  anisotropy  of 
the  transfer  functions  in  each  projection  by  carefully  calculating  the  incident  angle  at  each  point  on  the 
detector. 

•  My  research  has  been  the  first  to  show  that  DBT  is  capable  of  super-resolution  ( i.e .,  sub-pixel  resolution) 
using  analytical  modeling  and  experimental  reconstructions  of  bar  patterns.11  The  clinical  impact  of  super¬ 
resolution  is  optimizing  the  visibility  of  microcalcifications  in  reconstructions  performed  on  grids  with 


0°  Pitch  30°  Pitch 


Figure  4.  Two  reconstructions 
of  microcalcifications  are  shown. 
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much  smaller  pixilation  than  the  detector.  The  existence  of  super-resolution  was  later  generalized  to 

i  /r  i  n 

reconstruction  planes  with  a  broad  range  of  pitches  relative  to  the  breast  support.  ’ 

•  Much  of  my  research  effort  has  also  focused  on  the  development  of  optimization  strategies  for  DBT,  as 
described  in  SA2  of  my  predoctoral  DOD  grant.  To  date,  two  optimization  strategies  have  been  proposed. 
First,  the  phosphor  thickness  which  optimizes  DQE  has  been  calculated  for  various  projection  angles.4 

Second,  the  optimal  scan  time  has  been  evaluated  in  systems  with  patient  motion  and  either  continuous  tube 

20 

motion  or  patient  motion. 

IV.  REPORTABLE  OUTCOMES 

•  Acciavatti  RJ,  Maidment  ADA.  Optimization  of  phosphor-based  detector  design  for  oblique  x-ray  incidence 
in  digital  breast  tomosynthesis.  Medical  Physics.  201 1;38(1 1):6188-202. 

•  Acciavatti  RJ,  Maidment  ADA.  Optimization  of  In-Plane  Resolution  in  Digital  Breast  Tomosynthesis  Using 
Super-Resolution  Image  Reconstruction.  Medical  Physics.  2012:  to  be  submitted  for  review. 

•  Acciavatti  RJ,  Maidment  ADA.  TU-A-301-07:  Experimental  and  Theoretical  Validation  of  Breast 
Tomosynthesis  Reconstructions  along  Oblique  Planes.  Medical  Physics.  201 1;38(6):3746. 

•  Acciavatti  RJ,  Maidment  ADA.  An  Analysis  of  Super-Resolution  in  Oblique  Reconstructions  for  Digital 
Breast  Tomosynthesis.  Era  of  Hope  Conference;  2011  (conference  abstract);  Orlando,  FL;  2011  (conference 
abstract). 

•  Acciavatti  RJ,  Maidment  ADA.  Optimization  of  Continuous  Tube  Motion  and  Step-and-Shoot  Motion  in 
Digital  Breast  Tomosynthesis  Systems  with  Patient  Motion.  In:  Pelc  NJ,  Nishikawa  RM,  Whiting  BR, 
editors.  SPIE;  2012;  San  Diego,  CA:  SPIE;  2012.  p.  (accepted). 

V.  CONCLUSION 

My  research  on  the  resolution  loss  due  to  oblique  x-ray  incidence  in  DBT  projections  (Section  II.  A.) 
can  effectively  be  considered  complete,  since  it  has  been  published  in  the  November  2011  issue  of  the  journal 
Medical  Physics ,4  As  I  proposed  in  SA1  of  my  predoctoral  DOD  grant,  this  article  demonstrates  the  spatial 
anisotropy  of  transfer  functions  in  individual  projections  by  carefully  modeling  the  incident  angle  at  each  point 
on  the  detector.  To  explore  optimization  strategies  for  DBT  as  I  proposed  in  SA2,  this  work  shows  that  the 
phosphor  thickness  which  maximizes  DQE  is  projection  angle  dependent. 

In  addition,  I  have  demonstrated  the  existence  of  super-resolution  in  DBT  (Section  II.  B.)  by  calculating 
the  reconstruction  of  a  sine  plate  whose  attenuation  coefficient  varies  with  position  at  a  frequency  higher  than 
the  detector  alias  frequency.  The  results  of  my  analytical  modeling  were  supported  by  evidence  of  super¬ 
resolution  in  experimental  images  of  bar  patterns  and  select  clinical  images  of  microcalcifications.  A  full 
manuscript  on  super-resolution  will  soon  be  submitted  to  the  Medical  Physics  journal  for  review.11  My  future 
work  on  this  project  will  involve  responding  to  referee  criticisms  if  the  manuscript  is  accepted  for  publication. 

My  most  recent  projects  on  oblique  reconstructions  (Section  II.  C.)  and  optimization  of  scan  time 
(Section  II.  D.)  have  been  accepted  for  presentation  at  conferences,  ’  ’  but  have  not  yet  been  formulated  as 
manuscript  submissions  to  a  journal.  In  the  coming  months,  it  will  be  important  to  expand  upon  the 
fundamental  ideas  of  these  projects  in  order  to  develop  peer-reviewed  manuscript  submissions.  Although  my 
current  work  on  oblique  DBT  reconstructions  studies  the  effect  of  orienting  the  reconstruction  plane  at  various 
angles  relative  to  the  chest  wall  side  of  the  breast  support,  it  will  be  useful  to  model  the  effect  of  orienting  the 
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plane  at  various  angles  relative  to  the  chest  wall-to-nipple  direction  as  well.  In  my  work  on  scan  time 
optimization,  I  plan  to  model  patient  motion  velocities  that  vary  with  time,  such  as  the  pulsatile  motion  of 
objects  near  a  blood  vessel,  in  addition  to  the  constant  velocities  that  are  simulated  in  my  SPIE  proceedings 
manuscript.20 

I  have  completed  all  non-thesis  requirements  for  my  doctoral  degree.  By  the  end  of  2012,  I  plan  to 
complete  my  thesis.  This  time  frame  was  chosen  in  coordination  with  the  termination  of  this  DOD  predoctoral 
training  grant.  My  thesis  committee  met  on  January  5,  2012  and  agreed  with  this  time  frame. 
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Purpose:  In  digital  breast  tomosynthesis  (DBT),  a  volumetric  reconstruction  of  the  breast  is  gener¬ 
ated  from  a  limited  range  of  x-ray  projections.  One  trade-off  of  DBT  is  resolution  loss  in  the  projec¬ 
tions  due  to  non-normal  (i.e.,  oblique)  x-ray  incidence.  Although  degradation  in  image  quality  due 
to  oblique  incidence  has  been  studied  using  empirical  data  and  Monte  Carlo  simulations,  a  theoreti¬ 
cal  treatment  has  been  lacking.  The  purpose  of  this  work  is  to  extend  Swank’s  calculations  of  the 
transfer  functions  of  turbid  granular  phosphors  to  oblique  incidence.  The  model  is  ultimately  used 
as  a  tool  for  optimizing  the  design  of  DBT  detectors. 

Methods:  A  quantum-limited  system  and  20  keV  x-rays  are  considered.  Under  these  assumptions, 
the  modulation  transfer  function  (MTF)  and  noise  power  spectra  (NPS)  are  derived  using  the  diffu¬ 
sion  approximation  to  the  Boltzmann  equation  to  model  optical  scatter  within  the  phosphor.  This 
approach  is  applicable  to  a  nonstructured  scintillator  such  as  gadolinium  oxysulhde  doped  with  ter¬ 
bium  (Gd202S:Tb),  which  is  commonly  used  in  breast  imaging  and  which  can  reasonably  approxi¬ 
mate  other  detector  materials.  The  detective  quantum  efficiency  (DQE)  is  then  determined  from  the 
Nishikawa  formulation,  where  it  is  written  as  the  product  of  the  x-ray  quantum  detection  efficiency, 
the  Swank  factor,  and  the  Lubberts  fraction.  Transfer  functions  are  calculated  for  both  front-  and 
back-screen  configurations,  which  differ  by  positioning  the  photocathode  at  the  exit  or  entrance 
point  of  the  x-ray  beam,  respectively. 

Results:  In  the  front-screen  configuration,  MTF  and  DQE  are  found  to  have  considerable  angular 
dependence,  while  NPS  is  shown  to  vary  minimally  with  projection  angle.  As  expected,  the  high 
frequency  MTF  and  DQE  are  degraded  substantially  at  large  angles.  By  contrast,  ah  transfer  func¬ 
tions  for  the  back-screen  configuration  have  the  advantage  of  significantly  less  angular  dependence. 

Using  these  models,  we  investigated  the  possibility  for  optimizing  the  design  of  DBT  detectors.  As 
an  example  optimization  strategy,  the  phosphor  thickness  which  maximizes  the  DQE  at  a  fixed  fre¬ 
quency  is  analyzed.  This  work  demonstrates  that  the  optimal  phosphor  thickness  for  the  front- 
screen  is  angularly  dependent,  shifting  to  lower  thickness  at  higher  angles.  Conversely,  the  back- 
screen  is  not  optimized  by  a  single  thickness  but  instead  attains  reasonably  high  DQE  values  over  a 
large  range  of  thicknesses.  Although  the  back-screen  configuration  is  not  suited  for  current  detec¬ 
tors  using  a  glass  substrate,  it  may  prove  to  be  preferred  in  future  detectors  using  newly  proposed 
plastic  thin-hlm  transistor  (TFT)  substrates. 

Conclusions:  Using  the  diffusion  approximation  to  the  Boltzmann  equation  to  model  the  spread  of 
light  in  a  scintillator,  this  paper  develops  an  analytical  model  of  MTF,  NPS,  and  DQE  for  a  phos¬ 
phor  irradiated  obliquely.  The  model  is  set  apart  from  other  studies  on  oblique  incidence  in  being 
derived  from  first  principles.  This  work  has  applications  in  the  optimization  of  DBT  detector 
design.  ©2011  American  Association  of  Physicists  in  Medicine.  [DOI:  10.1118/1.3639999] 

Key  words:  oblique  x-ray  incidence,  digital  breast  tomosynthesis  (DBT),  optical  transfer  function 
(OTF),  noise  power  spectra  (NPS),  detective  quantum  efficiency  (DQE) 


I.  INTRODUCTION 

In  many  radiographic  studies,  non-normally  (i.e.,  obliquely) 
incident  x-rays  provide  a  source  of  blurring  at  the  periphery 
of  the  detector  due  to  the  divergence  of  the  x-ray  beam  emit¬ 
ted  from  the  focal  spot.  Que  and  Rowlands  proposed  an  ana¬ 
lytical  model  of  the  resolution  loss  due  to  oblique  incidence 
by  deriving  an  expression  for  the  modulation  transfer  func¬ 
tion  (MTF)  of  amorphous  selenium  (<2-Se)  detectors  from 


first  principles.1  Their  work  assumes  that  the  detector  is 
operated  in  drift  mode,  so  that  the  point  spread  function 
(PSF)  for  normal  incidence  is  a  delta  function  and  hence  the 
MTF  for  normal  incidence  is  unity  at  ah  frequencies.2 

Oblique  incidence  is  more  readily  observed  in  digital 
mammography  (DM)  than  many  other  imaging  studies.  A 
DM  detector  is  placed  closer  to  the  focal  spot  than  most 
modalities3  to  counteract  the  loss  in  x-ray  penetration  result¬ 
ing  from  the  use  of  relatively  low  energies  (^20  keV).4-6 
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The  drawback  of  decreasing  the  source-to-detector  distance 
is  increasing  the  angle  of  incidence  relative  to  the  normal, 
especially  at  the  edges  of  the  detector.  For  example,  the 
maximum  angle  of  incidence  is  25°  for  a  DM  system  with  a 
detector  field-of-view  (FOV)  of  24  x  30  cm  and  a  source-to- 
image  distance  of  70  cm  measured  at  the  midpoint  of  the 
chest  wall. 

In  digital  breast  tomosynthesis  (DBT),  low  dose  x-ray 
projection  images  are  acquired  over  a  limited  range  of  angles 
around  the  breast,  and  sharply  in-focus  slices  at  all  depths  of 
the  breast  volume  are  generated  using  image  reconstruction 
techniques.  Preliminary  studies  indicate  that  DBT  provides 
increased  sensitivity  and  specificity  relative  to  DM  for  the 
early  detection  of  breast  cancer  in  women.7  One  shortcoming 
of  DBT,  however,  is  that  it  is  more  directly  impacted  by  the 
resolution  loss  due  to  oblique  incidence  than  DM.  The  maxi¬ 
mum  projection  angle  in  DBT  can  be  as  large  as  20°  or  30°, 
and  the  angle  of  incidence  at  the  edges  of  the  detector  is 
even  higher  if  one  takes  into  account  the  divergent  x-ray 
beam  geometry.  While  some  DBT  systems  incorporate  a 
rotating  detector  to  counteract  changes  in  obliquity,  many 
systems  employ  a  stationary  detector. 

Mainprize  et  al.  experimentally  demonstrated  the  resolu¬ 
tion  loss  due  to  oblique  incidence  in  cesium  iodide  doped 
with  thallium  (CsI:Tl),  a  structured  phosphor-based  detector, 
using  the  slanted  edge  technique  to  measure  MTF.8  The 
authors  showed  that  at  10°  incidence,  the  MTF  degradation 
becomes  comparable  to  the  resolution  loss  associated  with 
other  common  sources  of  image  blurring,  such  as  the  blur¬ 
ring  of  the  focal  spot  and  the  lateral  spread  of  visible  light 
within  the  scintillator.  At  40°  incidence,  the  MTF  is  reduced 
considerably;  for  example,  at  5  line  pairs  per  mm  (lp/mm), 
the  MTF  is  degraded  by  35%^10%  over  a  broad  range  of 
kVp  and  target-filter  combinations. 

While  Mainprize  et  al.  did  not  measure  noise  power  spec¬ 
tra  (NPS),  Hajdok  and  Cunningham  have  calculated  NPS 
using  Monte  Carlo  simulations  of  a- Se.3  Their  work  demon¬ 
strated  that  unlike  MTF,  NPS  has  minimal  angular  depend¬ 
ence.  Since  the  detective  quantum  efficiency  (DQE)  is 
dependent  upon  the  square  of  the  MTF,  Hajdok  and  Cun¬ 
ningham  have  shown  that  the  DQE  degradation  with  oblique 
incidence  at  high  frequencies  is  more  pronounced  than  the 
MTF  degradation. 

Although  oblique  incidence  has  been  studied  experimen¬ 
tally  and  using  Monte  Carlo  simulations,  a  theoretical  treat¬ 
ment  has  been  lacking.  For  this  reason,  the  purpose  of  this 
paper  is  to  extend  Swank’s  calculations9  of  the  transfer  func¬ 
tions  of  turbid  granular  phosphors  to  oblique  incidence. 
Building  off  our  previous  work  on  oblique  incidence,10  we 
analyze  the  light  diffusion  equations  in  a  nonstructured  scin¬ 
tillator  such  as  gadolinium  oxysulhde  doped  with  terbium 
(Gd202S:Tb),  which  is  commonly  used  in  breast  imaging 
and  which  can  reasonably  approximate  other  detector  mate¬ 
rials.  The  theoretical  formulation  of  this  work  differs  from 
the  one  proposed  by  Que  and  Rowlands  in  not  making  the 
assumption  that  the  PSF  for  normal  incidence  is  a  delta  func¬ 
tion.  Ultimately,  the  analytical  model  is  used  as  a  tool  for 
optimizing  the  design  of  the  phosphor  for  oblique  incidence. 


II.  METHODS 

II.A.  Transfer  functions  for  front-screen  configuration 

The  optical  transfer  function  (OTF),  NPS,  and  DQE  of  a 
turbid  granular  phosphor  are  derived  here  from  first  princi¬ 
ples  for  all  angles  of  incidence.  The  Boltzmann  transport 
equation  may  be  used  to  model  the  spread  of  visible  light  in 
a  turbid  phosphor.  A  first-order,  steady  state  solution  to  the 
Boltzmann  transport  equation  is  a  diffusion  equation  of  the 
form11 

-  V2cf>( r)  +  <r2<Kr)  =  S(r),  (1) 

where  </>( r)  is  the  product  of  the  density  of  the  secondary  car¬ 
riers  (i.e.,  the  optical  photons)  with  the  diffusion  constant,  a 
is  the  reciprocal  of  the  mean  diffusion  length  of  the  second¬ 
ary  carriers,  and  S(r)  is  the  source  function.  The  diffusion 
equation  is  a  useful  approximation  to  the  Boltzmann  trans¬ 
port  equation  provided  that  three  criteria  are  met:  (1)  solu¬ 
tions  for  (/)(?)  are  determined  far  from  the  x-ray  source  S( r) 
relative  to  the  mean  free  path  of  optical  scatter;  (2)  the  opti¬ 
cal  properties  of  the  phosphor  possess  no  discontinuities; 
and  (3)  the  probability  of  optical  absorption  is  small  com¬ 
pared  against  the  probability  of  optical  scatter.12  This  model 
has  been  used  by  Swank  for  normal  incidence  9  and  has 
shown  good  agreement  with  experimental  data.12,13 

As  shown  in  Fig.  1,  the  source  function  S( r)  in  Eq.  (1) 
may  be  modeled  as  the  point  (z0  tan  0,  0,  z0)  along  the  x-ray 
path  length,  where  z0  is  depth  within  the  scintillator  of  total 
thickness  T  and  where  0  is  the  angle  of  incidence  relative  to 
the  normal.  In  terms  of  delta  functions,  S( r)  can  be  written 
as 


S'(r)  =  5(x  —  z0  tan0)<5(y)<5(z  —  z0).  (2) 

Using  the  Fourier  representation  of  the  delta  function,14  the 
source  function  can  equivalently  be  expressed  as  the  integral 


S(r) 


»oo 

e2ni{(x-z0  tan  _  (3) 

J  — oo 


Defining  v  as  the  2D  spatial  frequency  vector  with  compo¬ 
nents  vx  and  vy,  solutions  to  Eq.  (1)  can  be  written  in  the  form 


Photocathode 


Fig.  1 .  In  terms  of  delta  functions,  the  source  function  S(r)  at  the  depth  z0  of 
the  phosphor  is  found  from  trigonometry  to  be  S(x  —  z0  tan  6)S(y)S(z  —  z0), 
where  6  is  the  angle  of  x-ray  incidence  relative  to  the  normal.  The  figure 
assumes  a  front-screen  configuration  in  which  x-rays  are  incident  on  the 
backing  at  z  =  0  before  striking  the  photocathode  at  z  =  T.  Reversing  the 
direction  of  the  arrowhead  of  the  x-ray  beam  converts  the  front-screen  con¬ 
figuration  to  a  back-screen  configuration. 
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(f>(x,y,z)  = 


(4) 


Substituting  Eqs.  (3)  and  (4)  into  Eq.  (1),  it  can  be  shown  that 

d2\l/] 


±  +  q2^  =  e-<Wan ^(z  _  Z())} 


dz2 

where 

q2  =  a2+k2x+k2y, 


(5) 


(6) 


k  =  27tv. 


(7) 


To  solve  Eq.  (5)  for  i/rk(z),  one  can  apply  integral  transform 
techniques.15,16  Denoting  the  Laplace  transform  of  i/rk(z)  as 
£k(p),  the  transform  of  the  differential  equation  is 

(- p 2  +  q 2)  •  C^(p)  +  Cip  +  C2  =  e-ik*z°tw0e-pz°,  (8) 

where  C\  and  C2  are  the  constants  of  integration.  Solving  for 
(p)  and  taking  the  inverse  transform  generates  the  fol¬ 
lowing  piece-wise  expression  for  i/rk(z): 


C  2 

Ci  cosh(gz)  H - sinh(gz),  0  <  z  <  zo 

'K(zi  =  {  c  e-ikxzo  tan  0 

C i  cosh(gz)  H — -sinh(gz) - sinh[g(z  —  zq)],  zq  <  z  <T 

C[  C[ 


(9) 


The  constants  C\  and  C2  can  now  be  determined  from 
boundary  conditions  concerning  secondary  carrier  currents 
directed  toward  the  planes  at  z  =  0  and  z  =  T.  In  terms  of  the 
inverse  relaxation  length,  t,  the  secondary  carrier  currents 
across  any  plane  of  constant  z  are 


plane  z  =  0,  as  diagrammed  schematically  in  Fig.  1.  The 
OTF  of  the  scattering  process,  G(v,  z0),  is  then  determined 
for  a  point  source  from  the  expression 

Pi 


,  .  d(j> 

0T  +  — 
dz 


(10) 


"  dcji 

<px  -  — 

(H) 

"  (q  +  xpo)e^~ikxtm0^zo  +  {q  —  xpo)e~^+ikxtaa0^z°  ' 

dz 

=  XPl 

id  +  ?Po)  (d  +  zpi  )e«T  -(q-  xp0)(q  -  xp1  )e~c‘T _ 

yieft(z)  **4 

y  right  (z)  =  — 


In  the  right-hand  side  of  the  two  equations,  the  first  term 
models  the  effusion  current,  while  the  second  term  comes 
from  Fick’s  law.  The  first  boundary  condition  is  determined 
by  the  reflectivity  r0  of  the  plane  at  z  =  0.  Noting  that 
7right(0)  =  fWieft(O),  one  finds 

d<j) 
dz 

where 


G(v,z0)  = 

Hence, 
G(v,z0) 


1  +  Pi 


#  k 


(16) 


z=T 


=  zPo^lz=0> 


:=0 


Pi)  = 


1  ~  rp 

1  +  V 


(12) 


(13) 


The  second  boundary  condition  is  determined  from  the 
reflectivity  rx  of  the  boundary  at  z  =  T,  as  stipulated  by  the 
expression  j\eft(T)  =  rx  7right(C).  Defining  px  similar  to  p0  and 
noting  that  the  boundary  conditions  hold  for  each  Fourier 
component  i/rk  of  </>,  it  can  be  shown  that 


Ci  = 


{q  +  Tpx)eq(T  z°)  +  (q  —  Tpx)e  q^T  z° ) 


Id  +  t p0)(q  +  r p^7  -{q-  xp0){q  -  x px)e~‘iT 

—ikxZQ  tan  0 


x  e 


C2  —  TpqCi. 


(14) 

(15) 


Consistent  with  Swank’s  approach,  the  photocathode  is 
defined  by  the  plane  z  =  T  and  the  backing  is  defined  by  the 


(17) 

To  calculate  the  OTF  of  the  entire  phosphor,  one  multiplies 
Eq.  (17)  by  the  relative  x-ray  signal  as  a  function  of  the  depth  z0 

^-jtzosec  esecg 


Nf(z 0)  =  : 


(18) 


Y sec  ^  5 

where  p  is  the  linear  attenuation  coefficient  of  the  phosphor, 
and  then  integrates  over  the  phosphor  thickness.  Assuming  a 
front- screen  (F)  configuration  in  which  x-rays  are  first  inci¬ 
dent  on  the  backing  at  z  =  0  before  striking  the  photocathode 
at  z  —  r,  the  OTF  is  thus 

rT 


Gf(v)  = 


NF(z0)G(v,z0)dz0 


(19) 


pp  sec  8 
Y _ g—pr sec  o 


\q  +  xpo)(e^y-~ikxtm0^T  —  1) 


y_  —  ikx  tan  9 


(q  -  xpo){e-^+ik*tw0'>T  -  1) 


y+ 


where 

p= 


ikx  tan  9 


xpi 


(20) 


(21) 


(<?  +  zpo) (<?  +  TPi)e“T  ~  (d -  TPo)(<7 -  *Pi  )e~“T ’ 
y±=q±psQc9.  (22) 

The  MTF  is  found  from  the  normalized  modulus  of  the  OTF.17 

In  the  absence  of  outside  noise  sources,  the  quantum  NPS 
or  W/7(v)  is  calculated  by  integrating  the  product  of  NF(z0) 
with  |G(v,  z0)  | 2  from  z0  =  0  to  z0  =  T. 
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NF(z0)\G(v,z0)\2dz0 


(23) 


P  p  sec  8 

1  _ g~nT  sec  6 


(q  +  Tp0)2(e^q+y-^T  —  1)  2{q2  —  t2pI)(1  —  e  ^rsec0)  (g  —  Tp0)2(l  —  e  ^ +y+)r ) 


<7  +  7- 


/isec  6 


q  +  y+ 


(24) 


With  the  OTF  and  quantum  NPS  known,  it  is  now  possible 
to  determine  the  DQE.  From  the  work  of  Nishikawa,  DQE 
can  be  formulated  as  the  product  of  four  terms12 

DQE(v)  =  AqAsRc(v)Rn(v):  (25) 


where  Aq  is  the  x-ray  quantum  detection  efficiency  (QDE) 
determined  by  the  Lambert-Beer  Law  as 

Aq  =  l-e-^Tsec8,  (26) 


As  is  the  Swank  information  factor 


cm 

W>(0)  ’ 


(27) 


/?c(v)  is  the  Lubberts  fraction 


Rc(v) 


J_  Igf(v)I2 

As'  WF{v)  ’ 


(28) 


and  Rjy(y)  is  the  ratio  of  the  x-ray  quantum  noise  power  to 
the  total  noise  power.  Assuming  a  quantum- limited  imaging 
system,  R^iv)  is  taken  to  be  unity  in  this  work. 


II. B.  Transfer  functions  for  back-screen  configuration 

In  a  similar  fashion,  the  transfer  functions  for  a  back- 
screen  ( B )  configuration  can  be  calculated.  Unlike  the 
front-screen  configuration,  x-rays  first  strike  the  photoca¬ 
thode  at  z  =  T  before  passing  through  the  backing  at 
z  =  0.  This  modification  reverses  the  direction  of  the 
x-ray  beam  in  Fig.  1  without  further  altering  the  diagram. 
Hence 


Nb(z0) 


^-niT-z^secB  secQ 
Y _ sec  ^ 


(29) 


so  that  the  OTF  and  quantum  NPS  are  found  to  be 


/)/( sec  0  (q  +  Tp0)(e<y-  lkx tan  0,T  —  1 ) 
B  v  eiiT  sec  o  _  i  y+  —  ikx  tan  9 

(q-  Tp0)(e-^-+ik*tm8'>T  -  1)' 
y  _  +  ikx  tan  9 


(30) 


Wb(v) 


jg2/tsec  6  (q  +  Tp0)2(e^+y+'>T  -  1)  2(q2  -  x2pl)(eBTsece  -  l)  (q  -  rp0)2(  1  -  g-(g+y-)r) 

enT seed  _  i  q  +  y+  psecO  1  +  y_ 


(31) 


Equations  (30)  and  (31)  follow  from  Eqs.  (19)  and  (23);  the 
subscript  “5”  is  used  to  denote  a  back-screen. 

III.  RESULTS  FOR  A  MODEL  DETECTOR 

III.A.  Transfer  Functions  for  Front-  and  Back-Screen 
Configurations 

The  OTF,  NPS,  and  DQE  calculations  are  now  illustrated 
for  a  phosphor  with  a  reflective  backing  (r0  =  1),  a  nonreflec- 
tive  photocathode  ( rx  =  0),  and  optical  scatter  at  the  diffusion 
limit  (t  — ►  oo).  In  view  of  the  limitations  of  Swank’s  model, 
a  large  value  of  z  has  been  chosen.  As  Swank  demonstrated 
in  his  original  paper,9  the  MTF  of  a  scattering  phosphor 
(t  >  0)  is  always  higher  than  the  MTF  a  nonscattering  phos¬ 
phor  (t  =  0)  at  low  frequencies,  but  the  opposite  trend  arises 
at  high  frequencies.  The  crossover  point  of  the  scattering  and 
nonscattering  MTF  curves  corresponds  to  the  frequency 
beyond  which  Swank’s  model  becomes  inaccurate.  For  very 
small  values  of  t,  the  crossover  point  occurs  at  relatively 
low  frequencies.  However,  for  infinite  t,  the  scattering  MTF 


curve  never  crosses  the  nonscattering  MTF  curve  between  0 
and  10  lp/mm.  Hence,  Swank’s  model  is  increasingly  accu¬ 
rate  in  approaching  the  diffusion  limit. 

In  calculating  the  transfer  functions,  we  assume  20  keV 
monoenergetic  x-rays4-6  and  a  porous,  100  jum  thick 
Gd202S:Tb  phosphor  with  50%  packing  density.  The  attenu¬ 
ation  coefficient  p  for  the  porous  phosphor  is  determined  by 
halving  the  value  for  a  crystalline  phosphor.12  Since  crystal¬ 
line  Gd202S  has  a  mass  density  of  7.34  g/cm3  and  a  mass 
attenuation  coefficient  of  36.9  cm2/g  for  20  keV  x-rays,18 
the  attenuation  coefficient  for  porous  Gd202S  is  13.5  mm-1. 

In  Fig.  2,  cross  sections  of  the  MTF  surface  are  plotted 
versus  frequency  at  two  polar  angles  (a)  of  the  frequency 
vector  (0°  and  90°)  for  multiple  angles  of  incidence  and  two 
optical  absorption  parameters.  The  value  of  the  high  optical 
absorption  parameter  (a  =  20  mm-1)  was  chosen  to  match 
Swank’s  example9  in  which  oT  =  2.  In  practice,  the  optical 
absorption  can  be  increased  by  adding  an  optical  dye  to  the 
phosphor.  Following  convention,  the  polar  angle  is  defined 
as  the  angle  of  the  frequency  vector  relative  to  the  x  axis,  so 
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Spatial  Frequency  (Ip/mm)  Spatial  Frequency  (Ip/mm) 

Fig.  2.  The  modulation  transfer  function  (MTF)  of  a  porous  Gd202S:Tb  phosphor  is  plotted  versus  frequency  measured  along  0°  and  90°  polar  angles  for  mul¬ 
tiple  incident  angles  (0  =  0°,  10°,  20°,  30°)  and  two  optical  absorption  parameters  (c  =  0,  20  mm-1).  The  scintillator  possesses  a  reflective  backing,  a  nonre- 
flective  photocathode,  optical  scatter  at  the  diffusion  limit,  and  quantum-limited  noise.  Also,  the  phosphor  thickness  is  100  /m i,  and  the  incident  x-ray  energy 
is  20  keV.  As  shown,  the  front-screen  configuration  has  considerably  more  angular  dependence  than  the  back-screen  configuration. 


that  the  0°  polar  angle  is  only  perpendicular  to  the  x-ray 
beam  for  normal  incidence  and  the  90°  polar  angle  is  always 
perpendicular  to  the  x-ray  beam  (Fig.  1). 

Consistent  with  Swank’s  work  at  normal  incidence, 
Fig.  2  demonstrates  that  increasing  the  optical  absorption 
increases  the  MTF.  Figure  2(a)  indicates  that  increasing  the 
angle  of  incidence  decreases  the  MTF,  giving  rise  to  poorer 
spatial  resolution  in  the  front-screen  configuration.  For 
example,  comparing  30°  incidence  to  normal  incidence  at 
5.0  lp/mm  (a  =  0°),  the  MTF  decreases  by  17%  in  a  phos¬ 
phor  with  no  optical  absorption  and  by  15%  in  a  phosphor 
with  high  optical  absorption.  As  expected,  the  MTF  has 
minimal  angular  dependence  orthogonal  to  the  ray  of  inci¬ 
dence  [Fig.  2(c)]. 

Figures  2(b)  and  2(d)  demonstrate  that  the  back-screen 
configuration  has  superior  MTF  to  the  front- screen  configu¬ 
ration  for  all  projection  angles.  This  result  has  been  well- 
established  for  normal  incidence.19,20  More  significantly, 
Figures  2(b)  and  2(d)  further  show  that  the  angular  depend¬ 
ence  of  the  MTF  is  much  less  pronounced  in  the  back- screen 
configuration  than  in  the  front-screen  configuration.  For 
example,  comparing  30°  incidence  to  normal  incidence  at 


5.0  lp/mm  along  a  0°  polar  angle,  the  back-screen  MTF 
decreases  by  a  mere  3%.  Unlike  the  front-screen,  the  back- 
screen  MTF  increases  slightly  with  projection  angle  for 
measurements  orthogonal  to  the  incident  ray. 

In  Fig.  3,  normalized  NPS  (NNPS)  is  plotted  versus  fre¬ 
quency  for  the  same  scintillator.  Like  MTF,  NNPS  increases 
with  increasing  optical  absorption  for  all  angles  of  incidence. 
Unlike  MTF,  NNPS  is  independent  of  the  directionality  of 
the  frequency  vector.  Over  projection  angles  typical  of  DBT, 
the  angular  dependence  of  the  NNPS  is  minimal.  For 
example,  comparing  30°  incidence  to  normal  incidence  at 
5.0  lp/mm  in  a  front-screen  configuration,  NNPS  decreases 
by  8%  in  a  phosphor  with  no  optical  absorption  and  by  4% 
in  a  phosphor  with  high  optical  absorption.  In  a  back-screen 
configuration,  NNPS  increases  slightly  by  5%  and  2%, 
respectively. 

Figure  4  shows  DQE  versus  frequency.  In  both  a  front- 
and  back-screen  configuration,  Figures  4(a)  and  4(b)  demon¬ 
strate  that  for  measurements  made  along  a  0°  polar  angle, 
DQE  increases  with  projection  angle  at  low  frequencies  and 
decreases  with  projection  angle  at  high  frequencies.  At  low 
frequencies,  the  angular  dependence  of  the  x-ray  quantum 
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Front- Screen 


Spatial  Frequency  (Ip/mm) 


Fig.  3.  Normalized  noise  power  spectra  (NNPS)  is  plotted  versus  frequency.  NNPS  is  independent  of  the  directionality  of  the  frequency  vector.  While  front- 
screen  NNPS  decreases  with  projection  angle,  back-screen  NNPS  increases  slightly  with  projection  angle.  The  plots  implicitly  share  a  legend  with  Fig.  2. 


detection  efficiency  (Aq)  is  responsible  for  the  DQE 
increase.  In  calculating  Aq  using  Eq.  (26),  the  x-ray  path 
length  increases  from  the  phosphor  thickness  T  with  normal 


incidence  to  7sec0  with  oblique  incidence;  hence  a  greater 
number  of  x-rays  are  converted  to  visible  light.  At  high  fre¬ 
quencies,  the  degradation  in  DQE  with  increasing  projection 
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Fig.  4.  Detective  quantum  efficiency  (DQE)  is  plotted  versus  frequency.  Along  a  0°  polar  angle,  DQE  increases  with  projection  angle  at  low  frequencies  and 
only  decreases  with  projection  angle  at  high  frequencies.  The  front-screen  has  much  more  angular  dependence  than  the  back-screen  at  very  high  frequencies. 
For  measurements  orthogonal  to  the  incident  ray  (90°  polar  angle),  DQE  increases  with  projection  angle  over  a  very  broad  range  of  frequencies  in  both  config¬ 
urations.  The  plots  implicitly  share  a  legend  with  Fig.  2. 
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angle  arises  from  the  combined  angular  dependencies  of  the 
OTF  and  NPS.  The  high  frequency  DQE  for  the  front-screen 
configuration  is  lower  than  that  of  the  back- screen  configura¬ 
tion,  and  its  degradation  with  projection  angle  is  much  more 
pronounced.  For  example,  at  5.0  lp/mm  in  a  front-screen 
irradiated  at  a  30°  angle,  the  DQE  decreases  by  20%  relative 
to  normal  incidence.  In  the  back- screen  configuration,  the 
relative  decrease  in  DQE  is  less  than  5%. 

In  the  direction  orthogonal  to  the  incident  ray,  DQE 
increases  with  projection  angle  over  a  very  broad  range 
of  frequencies  [Figs.  4(c)  and  4(d)].  Comparing  30°  inci¬ 
dence  with  normal  incidence  at  5.0  lp/mm  in  a  front- 
screen  configuration,  DQE  increases  by  6%  in  a  phosphor 
with  no  optical  absorption  and  by  4%  in  a  phosphor  with 
high  optical  absorption.  Using  a  back-screen  configura¬ 
tion,  the  relative  increase  in  DQE  is  approximately  twice 
as  high. 

In  Fig.  5,  the  angular  dependence  of  the  Swank  factor 
(As)  used  for  calculating  the  DQE  is  studied.  Swank  has 
shown  that  As  provides  a  measure  of  the  fluctuation  in  sig¬ 
nal  generated  from  each  x-ray  photon  due  to  variability  in 
the  absorbed  energy  of  each  interacting  x-ray  and  in  the 
number  of  secondary  carriers  generated  from  each  interact¬ 
ing  x-ray.21  Figure  5  demonstrates  that  in  a  phosphor  with 
no  optical  absorption,  the  Swank  factor  is  unity  at  all  projec¬ 
tion  angles  for  either  the  front-  or  back-screen  configuration. 
By  contrast,  in  a  phosphor  with  high  optical  absorption,  the 
Swank  factor  has  slight  angular  dependence  over  projection 
angles  typical  of  DBT.  For  example,  comparing  30°  inci¬ 
dence  to  normal  incidence,  As  increases  by  0.3%  in  a  front- 
screen  configuration  and  by  0.5%  in  a  back-screen  configu- 


- Ap  From-  or  Back- Screen,  0  =  0 

- Av  Front-Screen,  o  =  20  mm1 

- Ap  Back-Screen,  o  =  20  mm  1 


DQL(0),  Front-  or  Back -Screen,  a  -  0 
— —— ■  DQE(0),  Front-Screen,  <r  -  20  mm  1 
*  PQR(O).  Rack-Screen,  c  =  20  mnr1 


ration.  At  very  oblique  angles  approaching  shearing  inci¬ 
dence  ( 9  =  90°),  the  Swank  factor  increases  sharply  to  unity. 

Unlike  As,  DQE(0)  is  projection  angle  dependent  for  all 
possible  optical  absorption  parameters.  For  both  configura¬ 
tions,  the  relative  increase  in  DQE(0)  comparing  30°  inci¬ 
dence  to  normal  incidence  is  6%  in  a  phosphor  with  no 
optical  absorption  and  7%  in  a  phosphor  with  high  optical 
absorption.  The  angular  dependence  of  DQE(0)  is  therefore 
more  pronounced  than  the  Swank  factor. 


III.B.  Anisotropy  of  the  transfer  functions  over  the 
detector  area 

Because  the  focal  spot  of  a  DBT  system  emits  x-rays  in 
all  directions,  the  angle  of  incidence  is  spatially  variant  at 
each  point  on  the  detector.  Assuming  a  stationary  detector 
whose  center-of-rotation  (COR)  coincides  with  the  midpoint 
of  the  chest  wall,  the  angle  of  incidence  relative  to  the  nor¬ 
mal  at  each  point  (x,  y)  on  the  detector  may  be  determined 
from  the  expression 


9  =  arctan 


—  h  sin  A)2  +  y 2 
h  cos  A 


(32) 


where  h  is  the  source-to-COR  distance  and  A  is  the  nominal 
projection  angle  (i.e.,  the  angle  of  the  x-ray  tube  arc  relative 
to  the  normal  at  the  COR).  In  deriving  this  result,  the  phos¬ 
phor  thickness  (T)  is  taken  to  be  negligible  compared  against 
the  source-to-COR  distance  (h).  Also,  the  chest  wall  defines 
the  v  axis  of  the  detector  and  its  midpoint  the  origin.  For  a 
24  x  30  cm  field-of-view  (FOV)  and  a  source-to-COR  dis¬ 
tance  of  70  cm,  the  angle  9  is  plotted  versus  the  length  (x) 
and  the  width  (y)  of  the  detector  for  the  central  projection 
(A  =  0°)  and  an  oblique  projection  (A  =  20°)  in  Fig.  6.  In  the 
central  projection  [Fig.  6(a)],  the  angle  of  incidence  relative 
to  the  normal  varies  between  0°  and  25°,  and  in  the  oblique 
projection  [Fig.  6(b)],  it  varies  between  10°  and  35°.  For  ei¬ 
ther  case,  the  maximum  angle  is  found  at  a  corner  of  the 
FOV  opposite  the  chest  wall. 

Fike  the  angle  of  incidence,  the  transfer  functions  of  the 
phosphor  are  spatially  variant  across  the  detector  area 
(Fig.  7).  To  illustrate  the  spatial  anisotropy  of  one  of  the 
transfer  functions,  a  surface  plot  of  front- screen  DQE  versus 
position  along  the  detector  is  shown  for  the  frequency 
5.0  lp/mm  in  a  phosphor  with  high  optical  absorption 
assuming  a  0°  polar  angle  for  the  frequency  vector.  In  the 
central  projection  [Fig.  7(a)],  the  front-screen  DQE  varies 
between  0.34  and  0.30  (12%  decrease),  and  in  the  oblique 
projection  [Fig.  7(b)],  it  varies  between  0.33  and  0.24  (27% 
decrease).  Although  not  plotted  in  the  figure,  one  can  show 
that  the  back-screen  DQE  has  much  less  variation  over  the 
detector  area  in  either  projection  (<4%  decrease). 


Fig.  5.  The  Swank  information  factor  (As)  has  no  angular  dependence  in  a 
phosphor  lacking  optical  absorption  and  slight  angular  dependence  over  typ¬ 
ical  incident  angles  in  a  phosphor  with  high  optical  absorption.  DQE(0)  has 
greater  relative  variation  with  incident  angle.  Both  As  and  DQE(0)  increase 
sharply  to  unity  at  angles  approaching  shearing  incidence  ( 6  =  90°). 


III.C.  Optimization  of  phosphor  thickness  for  oblique 
incidence 

In  addition  to  illustrating  the  impact  of  oblique  incidence 
on  the  transfer  functions  of  a  phosphor,  the  analytical  models 
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(a)  Central  Projection 


(b}  Oblique  Projection 


I 
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27* 
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Fig.  6.  The  angle  of  incidence  relative  to  the  normal  is  plotted  versus  position  along  the  detector  for  (a)  the  central  projection  and  (b)  an  oblique  projection. 
The  DBT  system  has  a  source-to-COR  distance  of  70  cm,  and  the  nominal  projection  angle  in  (b)  is  20°. 


developed  in  this  work  can  be  used  as  a  platform  for  opti¬ 
mizing  detector  design  over  the  range  of  projection  angles 
used  in  DBT.  One  important  element  in  the  design  of  a  phos¬ 
phor  is  its  x-ray  quantum  detection  efficiency  (QDE).  In 
Figs.  8  and  9,  DQE  at  a  fixed  frequency  is  plotted  versus 
QDE  at  normal  incidence  to  investigate  whether  DQE  can  be 
maximized  by  varying  QDE.  Both  figures  have  been  gener¬ 
ated  using  the  same  phosphor  parameters  analyzed  in 
Figs.  2-5,  except  the  phosphor  thickness  T  is  now  left  as  a 
variable  which  allows  QDE  at  normal  incidence  to  vary. 


In  Eq.  (33),  AQ0  denotes  the  QDE  at  normal  incidence. 

Figure  8  illustrates  that  in  a  front-  or  back- screen  configu¬ 
ration  with  no  optical  absorption,  DQE(O)  can  be  optimized 
by  manufacturing  a  phosphor  with  a  very  large  thickness 
(100%  QDE).  By  contrast,  in  a  phosphor  with  high  optical 
absorption,  the  dependency  of  DQE(O)  on  QDE  is  quite 
different  for  the  two  configurations.  In  a  front- screen, 
DQE(O)  is  maximized  by  an  intermediate  QDE  which  is  pro¬ 
jection  angle  dependent,  favoring  smaller  thicknesses  at 


larger  angles.  For  example,  in  the  0°  and  30°  projections,  the 
optimal  QDE  at  normal  incidence  are  0.73  and  0.72  corre¬ 
sponding  to  97  and  94  fim  thicknesses,  respectively.  By  con¬ 
trast,  the  back-screen  DQE(0)  attains  relatively  high  values 
over  a  broader  range  of  QDE.  With  100%  QDE,  DQE(0)  pla¬ 
teaus  to  0.64  and  0.68  for  the  0°  and  30°  projections, 
respectively. 

In  Fig.  9,  the  dependency  of  DQE  on  QDE  is  analyzed  at 
a  higher  frequency  (5.0  lp/mm).  The  high  frequency  DQE 
for  the  front- screen  is  maximized  at  an  intermediate  QDE 
for  both  optical  absorption  parameters.  For  measurements 
made  along  the  0°  polar  angle  in  a  front-screen  with  no  opti¬ 
cal  absorption,  the  optimal  QDE  values  at  normal  incidence 
are  0.60  for  the  0°  projection  and  0.54  for  the  30°  projection, 
corresponding  to  68  and  58  fim  thicknesses,  respectively 
[Fig.  9(a)].  With  high  optical  absorption,  the  respective  QDE 
optima  are  0.55  and  0.51  (59  and  52  fim  thicknesses).  In  the 
direction  orthogonal  to  the  incident  ray,  the  optimal  QDE 
have  less  projection  angle  dependence  [Fig.  9(c)]. 

Unlike  the  front-screen  configuration,  the  back-screen 
configuration  supports  relatively  high  DQE  over  large  QDE 
values  at  5.0  lp/mm.  For  measurements  made  along  the  0° 
polar  angle  [Fig.  9(b)],  the  back-screen  DQE  plateaus  to 


(a)  Genual  Projection,  a  =  0°  (b) 


Oblique  Projection,  a  =  0° 
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Fig.  7.  The  anisotropy  of  the  transfer  functions  over  the  detector  is  illustrated  by  plotting  the  front-screen  DQE  at  a  fixed  spatial  frequency  (5.0  lp/mm)  versus 
areal  position,  assuming  high  optical  absorption  (a  =  20  mm-1)  and  frequency  measurements  along  a  0°  polar  angle.  The  oblique  projection  has  greater  varia¬ 
tion  in  DQE  over  the  detector  area  than  the  central  projection.  The  directionalities  of  the  x  and  y  axes  are  flipped  relative  to  Fig.  6  to  improve  visualization  of 
the  surfaces. 
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Fig.  8.  For  four  angles  of  incidence  and  two  optical  absorption  parameters,  the  dependency  of  DQE(O)  on  QDE  at  normal  incidence  is  analyzed.  With  no  opti¬ 
cal  absorption,  DQE(O)  for  both  front-  and  back-screen  configurations  is  optimized  by  large  QDE.  With  high  optical  absorption,  the  maximum  DQE(O)  for  the 
front-screen  occurs  at  intermediate  QDE,  while  relatively  high  DQE(O)  for  the  back-screen  occurs  over  a  broad  range  of  large  QDE.  The  optimal  QDE  for  the 
front-screen  is  angularly  dependent. 


0.51  and  0.48  for  the  0°  and  30°  projections,  respectively,  in 
a  phosphor  with  no  optical  absorption  (0.46  and  0.45  with 
high  optical  absorption).  For  measurements  made  along  the 


orthogonal  direction,  the  back- screen  DQE  attains  a  higher 
plateau;  also,  the  asymptote  actually  increases  slightly  with 
projection  angle  [Fig.  9(d)]. 


QDE  at  Normal  Incidence 


QDE  at  Nonna  I  Incidence 


Fig.  9.  At  5.0  lp/mm,  DQE  is  plotted  versus  QDE  at  normal  incidence.  In  the  front-screen  configuration,  DQE  is  optimized  by  an  intermediate  QDE.  The  opti¬ 
mal  QDE  is  projection  angle  dependent,  shifting  to  lower  values  (thinner  phosphors)  at  larger  angles.  By  contrast,  in  the  back-screen  configuration,  relatively 
high  DQE  is  supported  over  large  QDE  values  for  all  projection  angles.  The  plots  implicitly  share  a  legend  with  Fig.  8. 
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Fig.  10.  A  surface  plot  shows  the  dependence  of  DQE  on  both  frequency  and  QDE  at  normal  incidence,  assuming  6  =  30°,  a  =  20  mm-1,  and  a  0°  polar  angle 
for  the  frequency  vector.  At  all  frequencies,  the  front-screen  DQE  is  optimized  by  an  intermediate  QDE.  By  contrast,  the  back-screen  DQE  attains  relatively 
high  values  over  a  broad  range  of  large  QDE  (thick  phosphors). 


To  illustrate  the  dependence  of  DQE  on  both  frequency 
and  QDE,  a  surface  plot  is  shown  in  Fig.  10  assuming  30° 
incidence,  high  optical  absorption,  and  frequency  measure¬ 
ments  along  the  0°  polar  angle.  The  curvature  of  the  front- 
screen  surface  demonstrates  that  the  value  of  QDE  which 
maximizes  DQE  decreases  with  frequency.  By  contrast,  the 
back- screen  surface  shows  that  DQE  is  not  optimized  by  a 
single  QDE  value.  Instead,  at  all  frequencies,  the  back- 
screen  supports  relatively  high  DQE  over  very  large  QDE. 

Figures  8-10  demonstrate  that  the  optimal  thickness  of  a 
front- screen  is  both  projection  angle  dependent  and  fre¬ 
quency  dependent.  In  Fig.  11,  the  combined  dependence  is 
shown  explicitly  in  a  surface  plot,  assuming  a  0°  polar  angle 
and  high  optical  absorption.  The  graph  was  generated  in  mat- 
lab  R2010b  by  discretizing  a  grid  (60  x  60)  of  incident 
angles  and  frequencies  from  0°  to  45°  and  0  to  10  lp/mm, 
respectively.  The  optimal  thickness  which  maximizes  DQE 
was  determined  by  the  zero  of  the  first  partial  derivative  of 
DQE  with  respect  to  phosphor  thickness.  Because  the  zeros 
of  the  first  DQE  derivative  cannot  be  easily  solved  in  closed 


Front-Screen,  a  -  0° 


Spatial  Frequency  Incident 
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Fig.  11.  For  a  front-screen  configuration,  the  dependence  of  the  optimal 
phosphor  thickness  on  both  the  angle  of  incidence  and  frequency  is  ana¬ 
lyzed,  assuming  a  — 20  mm-1  and  a  0°  polar  angle  for  the  frequency  vector. 


form,  Newton’s  method  was  implemented  to  find  the  zeros 
numerically 


Tn+ 1  =  Tn  — 


dDp 

~d2DF 

dT 

T=Tn 

dT2 

T=Tn 

n  G  Jf, 


(34) 


where  DF  is  the  symbolic  abbreviation  for  front-screen 
DQE.  For  all  projection  angles  and  frequencies  investigated 
in  the  plot,  the  initial  guess  (25  fim)  and  the  number  of  itera¬ 
tions  (9)  provided  convergence  exceeding  ten  decimal  pla¬ 
ces.  Figure  11  demonstrates  that  the  optimal  phosphor 
thickness  for  the  front-screen  configuration  is  a  decreasing 
function  of  both  incident  angle  and  frequency,  ranging  from 
97  fim  (0°  incidence,  0  lp/mm)  to  25  fim  (45°  incidence, 
10  lp/mm). 


IV.  COMPARISON  WITH  RESULTS  IN  THE 
LITERATURE 

This  paper  extends  Swank’s  calculations9  of  the  transfer 
functions  of  turbid  phosphors  to  oblique  x-ray  incidence.  In 
the  limiting  case  of  normal  incidence,  the  formulas  presented 
in  this  paper  exactly  reduce  to  Swank’s  results.  Our  work  is 
unique  in  modeling  the  transfer  functions  for  oblique  inci¬ 
dence  in  closed  form  without  making  the  assumption  that  the 
PSF  of  normal  incidence  is  a  delta  function.1  One  benefit  of 
this  approach  is  unifying  many  prior  results  on  oblique  inci¬ 
dence  under  one  model.  For  example,  in  a  front-screen  con¬ 
figuration,  we  demonstrate  that  oblique  incidence  degrades 
the  MTF,  and  that  the  resultant  loss  in  resolution  becomes 
more  pronounced  with  increasing  frequency  and  increasing 
angle.  Although  these  findings  are  derived  for  a  turbid  phos¬ 
phor  such  as  Gd202S:Tb,  they  are  consistent  with  experi¬ 
mental  data  on  CsI:Tl  presented  by  Mainprize  et  al .8  as  well 
as  analytical  modeling  of  a- Se  proposed  by  Que  and  Row¬ 
lands.1  On  a  similar  note,  we  have  observed  that  NPS  is 
degraded  with  increasing  projection  angle,  though  to  a  lesser 
degree  than  MTF.  This  Ending  is  concordant  with  the  prior 
work  of  Hajdok  and  Cunningham,3  who  performed  Monte 
Carlo  simulations  of  a- Se.  As  a  final  point,  we  have  shown 
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that  DQE  increases  with  incident  angle  at  low  frequencies 
but  only  decreases  with  incident  angle  at  high  frequencies. 
Consistent  with  the  findings  of  Hajdok  and  Cunningham,  the 
DQE  degradation  with  projection  angle  at  high  frequencies 
is  more  pronounced  than  the  MTF  degradation,  reflecting  the 
dependency  of  DQE  on  the  square  of  MTF. 

In  this  work,  it  has  been  observed  that  the  Swank  factor  is 
angularly  dependent,  but  that  its  variation  is  small  over  pro¬ 
jection  angles  typical  of  DBT.  In  particular,  it  has  been 
shown  that  the  Swank  factor  changes  by  no  more  than  0.5% 
comparing  0°  and  30°  incidence.  This  observation  is  consist¬ 
ent  with  Monte  Carlo  simulations  of  CsI:Tl  phosphors  con¬ 
ducted  by  Badano  et  al .,  who  demonstrated  that  the  variation 
in  the  Swank  factor  over  projection  angles  typical  of  DBT  is 
minimal.22  While  the  relative  change  in  the  Swank  factor 
with  projection  angle  is  small,  the  relative  increase  in 
DQE(O)  is  more  substantial,  as  it  includes  the  effect  of 
increasing  x-ray  quantum  detection  efficiency  with  increas¬ 
ing  projection  angle. 

Although  our  work  demonstrates  consistency  with  other 
studies  on  oblique  incidence,  it  is  important  to  identify  fun¬ 
damental  differences  between  the  detectors  addressed  in  the 
comparison.  This  work  models  a  turbid  phosphor  in  which 
visible  light  spreads  by  optical  scatter.  By  contrast,  prior 
studies  on  CsI:Tl  (Refs.  8  and  22)  assume  a  structured  phos¬ 
phor  in  which  needlelike  crystals  approximately  10  /urn  in  di¬ 
ameter  transmit  the  optical  photons  to  the  photocathode  by 
total  internal  reflection,  thereby  minimizing  the  lateral 
spread  of  visible  light.19,23  At  a  given  incident  angle,  struc¬ 
tured  phosphors  should  have  higher  MTF  than  turbid  phos¬ 
phors  for  this  reason.  Although  the  transfer  functions  of 
turbid  phosphors  are  different  from  structured  phosphors, 
this  work  demonstrates  that  their  angular  dependence  fol¬ 
lows  comparable  trends.  On  a  similar  note,  our  model  has 
shown  concordance  with  prior  studies  on  oblique  incidence 
in  <2-Se.1,3  In  a- Se,  an  absorbed  x-ray  ionizes  a  Se  atom,  free¬ 
ing  an  electron  and  a  hole  which  migrate  to  different  ends  of 
the  detector  due  to  an  applied  electric  field.19  In  drift  mode, 
the  electric  field  is  small  enough  so  that  the  electron  and 
hole  do  not  have  sufficient  kinetic  energy  to  ionize  Se  atoms 
and  create  an  avalanche  of  electron-hole  pairs.  Because  the 
electron  and  hole  migrate  in  a  nearly  perfect  orthogonal  path 
to  opposite  ends  of  the  detector,  the  MTF  of  a- Se  at  normal 
incidence  is  approximately  unity  for  all  frequencies.2  At 
higher  incident  angles,  the  MTF  of  a- Se  decreases  with  fre¬ 
quency,  but  is  expected  to  be  superior  to  a  turbid  phosphor 
since  there  is  no  lateral  spread  of  visible  light.  The  analytical 
model  of  a- Se  developed  by  Que  and  Rowlands1  can  effec¬ 
tively  be  derived  by  using  Eq.  (2)  for  the  source  function  but 
by  eliminating  Eq.  (1)  for  the  diffusion  of  secondary 
carriers. 

To  our  knowledge,  this  paper  is  the  first  to  investigate  the 
angular  dependence  of  the  transfer  functions  of  the  back- 
screen  configuration.  The  consistency  of  our  back- screen 
model  with  expected  trends  at  normal  incidence  helps  to  sug¬ 
gest  its  validity.  For  example,  in  accord  with  experimental 
data  in  turbid  phosphors,20  we  demonstrate  that  the  back- 
screen  has  higher  MTF  than  the  front-screen.  Because  opti¬ 


cal  photons  are  predominately  generated  near  the  x-ray  en¬ 
trance  surface  of  the  phosphor,  visible  light  exhibits  less 
lateral  spread  before  reaching  the  photocathode  if  the  back- 
screen  configuration  is  used. 

Consistent  with  prior  authors,  this  paper  demonstrates 
that  the  Swank  factor  of  the  back-screen  is  greater  than  the 
front-screen.  In  experimental  measurements  on  Gd202S:Tb 
phosphors  at  20  keV,  Trauernicht  and  Van  Metter  demon¬ 
strated  that  a  back-screen  has  approximately  5%  higher 
Swank  factor  than  a  comparable  front- screen. 24  Such  a  rela¬ 
tive  change  in  As  with  screen  configuration  matches  the 
results  presented  in  Fig.  5.  It  is  evident  from  Fig.  5  that  the 
benefits  of  the  back- screen  over  the  front- screen  should  hold 
with  higher  incident  angles. 

V.  DISCUSSION 

This  study  develops  an  analytical  model  of  the  transfer 
functions  of  turbid  phosphors  for  oblique  x-ray  incidence. 
The  results  of  the  model  are  consistent  with  prior  observa¬ 
tions  on  oblique  incidence  in  a  range  of  detector  types.  Hav¬ 
ing  an  analytical  model  of  the  transfer  functions  has  in  turn 
led  to  the  development  of  optimization  strategies  for  improv¬ 
ing  detector  design  in  DBT. 

One  area  for  optimizing  detector  design  is  the  choice  of  a 
front-  or  back- screen  configuration.  Although  it  is  well 
known  that  the  back-screen  has  greater  MTF  than  the  front- 
screen  at  normal  incidence,  back- screen  transfer  functions 
have  not  yet  been  analyzed  for  oblique  incidence.  In  this 
work,  it  is  suggested  that  the  back-screen  transfer  functions 
should  have  much  less  angular  dependence  than  the  front- 
screen,  and  consequently,  exhibit  less  variation  with  position 
along  the  detector.  As  a  result,  a  back- screen  configuration 
may  be  chosen  to  optimize  the  design  of  a  phosphor  for 
oblique  incidence.  Initially,  this  result  would  seem  to  have 
no  practical  impact,  since  back- screens  are  not  currently 
used  clinically  due  to  the  glass  substrate  of  the  thin-film  tran¬ 
sistor  (TFT)  array  for  digital  signal  readout.25,27  If  these 
detectors  were  operated  as  back- screens,  the  high  attenuation 
of  glass  (~0.7  mm  thick)  would  prevent  a  large  percentage 
of  x-rays  from  reaching  the  phosphor  itself,19  and  hence  both 
QDE  and  DQE  would  be  compromised.  Based  on  new 
research  on  flexible  organic  light-emitting  diode  (OFED) 
displays,  however,  TFT  arrays  may  soon  be  manufactured 
on  a  plastic  substrate.28,34  Because  plastic  is  much  less 
attenuating  than  glass,  operating  such  a  detector  in  a  back- 
screen  configuration  would  not  be  prohibitive.  Detectors 
manufactured  with  plastic  have  many  benefits  such  as  being 
bendable,  light-weight,  and  easy  to  transport.  These  future 
phosphor-based  detectors  should  preferentially  be  operated 
as  back- screens  in  order  to  optimize  detector  performance 
for  oblique  incidence. 

The  analytical  model  of  the  transfer  functions  was  ulti¬ 
mately  used  as  a  platform  for  optimizing  the  QDE  of  the  de¬ 
tector  for  oblique  incidence.  This  work  provides  a  method 
for  determining  the  QDE  that  maximizes  DQE  at  any  fre¬ 
quency  of  interest,  such  as  the  frequency  of  small  microcal¬ 
cifications  or  fine  cancerous  lesions  within  the  breast.  To  our 
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knowledge,  this  paper  is  the  first  to  show  that  the  optimal 
QDE  is  projection  angle  dependent  in  a  front- screen,  tending 
toward  lower  values  (thinner  phosphors)  with  increasing  pro¬ 
jection  angle.  Because  the  incident  angle  is  greatest  at  the 
periphery  of  the  detector  opposite  the  chest  wall,  a  corollary 
of  this  finding  is  that  one  beneficial  design  feature  would  be 
to  reduce  thickness  at  the  edges  and  corners  of  the  phosphor. 
Although  it  is  beyond  the  scope  of  this  paper  to  determine  a 
single  value  for  the  optimum  thickness  at  each  point  on  a 
DBT  detector,  future  work  should  be  directed  at  modeling 
the  transfer  functions  of  the  reconstruction35  and  optimizing 
thickness  to  maximize  the  3D  DQE  for  a  fixed  frequency  of 
interest.  In  experimental  practice,  the  optimal  thickness 
should  be  calculated  on  a  case-by-case  basis  for  each  detec¬ 
tor  under  consideration,  taking  into  account  the  unique  char¬ 
acteristics  of  the  imaging  system. 

In  this  paper,  it  was  demonstrated  that  the  back- screen 
DQE  is  not  optimized  by  a  single  QDE  but  instead  attains 
reasonably  high  values  over  a  broad  range  of  phosphor  thick¬ 
nesses.  Hence,  the  back-screen  configuration  has  an  impor¬ 
tant  benefit  over  the  front- screen  configuration:  future 
detectors  with  a  plastic  TFT  substrate  can  be  manufactured 
with  very  large  thickness  without  degradation  in  high  fre¬ 
quency  DQE  at  any  projection  angle. 

A  discussion  of  the  limitations  of  this  paper  and  directions 
for  future  modeling  are  now  noted.  One  assumption  made  in 
the  Results  section  is  that  the  incident  x-ray  beam  is  monoe- 
nergetic.  Since  the  phosphor  attenuation  coefficient  fi  is 
energy-dependent,  it  is  important  to  model  polyenergetic 
x-ray  spectra36-38  when  studying  the  phosphor  thickness 
which  maximizes  DQE  in  a  front-screen  (Fig.  11).  Future 
work  should  be  aimed  at  determining  if  the  angular  depend¬ 
ence  of  the  optimal  thickness  becomes  more  or  less  pro¬ 
nounced  upon  varying  the  kVp  and  the  target-filter 
combination.  Since  Fig.  1 1  was  generated  assuming  a  monoe- 
nergetic  x-ray  beam,  the  reader  should  take  caution  against 
directly  applying  the  optimal  thickness  calculations  to  any 
real  imaging  system. 

Although  a  relatively  low  energy  (^20  keV)  is  convention¬ 
ally  chosen  to  maximize  contrast  between  glandular  tissue 
and  cancer,4-6  it  would  be  useful  to  simulate  higher  energies 
found  in  dual  energy  contrast-enhanced  DBT  (DE  CE- 
DBT).39^4  In  DE  CE-DBT,  low  and  high  energy  images  are 
acquired  below  and  above  the  K  edge  of  iodine  (33.2  keV)  af¬ 
ter  iodinated  contrast  is  injected  into  the  blood  stream  45,46 
and  contrast  uptake  is  determined  using  weighted  logarithmic 
subtraction.  Contrast  uptake  can  be  used  to  quantify  blood 
flow  at  the  site  of  a  tumor,  which  exceeds  healthy  tissue.47 

Another  extension  of  this  work  would  be  to  model  detec¬ 
tor  pixilation  due  to  the  TFT  array  in  which  the  phosphor  is 
placed  in  optical  contact.25-27  In  this  setting,  the  PSF  is  the 
convolution  of  the  phosphor  blurring  function  with  the  de¬ 
tector  element  sampling  function.48-50  An  additional  subtlety 
that  may  be  modeled  is  the  blurring  of  the  focal  spot.  The 
MTF  of  the  focal  spot  is  degraded  with  increasing  focal  spot 
size,51  as  well  as  with  increasing  focal  spot  motion  during  a 
continuous  DBT  scan.35,52  Our  model  can  also  be  refined  by 
taking  into  account  detector  lag  and  ghosting.53-55  In  addi¬ 


tion,  one  can  incorporate  the  possibility  for  an  NPS  that  is 
not  quantum-limited  but  possesses  outside  noise  sources.13 

A  final  point  of  investigation  would  be  to  extend  this 
work  to  a  structured  phosphor,  such  as  CsI:Tl.  In  structured 
phosphors,  needlelike  crystals  transmit  the  optical  photons  to 
the  photocathode  by  total  internal  reflection  to  minimize  the 
lateral  spread  of  visible  light.  Structured  phosphors  have 
superior  spatial  resolution  to  nonstructured  phosphors  for 
this  reason.19,23  To  model  a  structured  phosphor,  the  bound¬ 
ary  conditions  for  the  secondary  carrier  currents  would  no 
longer  be  determined  exclusively  by  the  reflectivities  of  the 
backing  and  the  photocathode.  Instead,  boundary  conditions 
would  also  exist  for  the  reflectivities  of  the  walls  of  the  nee¬ 
dlelike  crystals.  Recently,  Freed  et  al.  have  proposed  an  ana¬ 
lytical  model56,57  of  the  PSF  of  CsI:Tl  using  a  competing 
approach;  their  model  was  later  generalized  to  direct- 
converting  detectors.58  A  key  step  in  the  derivation  of  their 
PSF  formula  is  considering  three  different  functional  forms 
(Gaussian,  exponential,  Lorentzian)  to  quantify  the  spread  of 
secondary  carriers  or  electron-hole  pairs  at  each  depth  of  the 
detector  material.  Empirically,  the  authors  find  that  the  Lor¬ 
entzian  provides  the  best  match  to  data  generated  from 
Monte  Carlo  simulations.  The  authors  then  determine  opti¬ 
mal  values  for  the  parameters  in  their  model  by  minimizing 
the  normalized  differences  between  the  analytical  technique 
and  Monte  Carlo  simulations.  Although  useful  in  providing  a 
closed  form  solution  for  the  PSF,  their  model  is  limited  in 
not  being  derived  from  first  principles.  By  contrast,  in  our 
current  work,  all  results  are  derived  from  first  principles, 
and  the  value  of  each  parameter  has  physical  significance.  It 
would  be  intriguing  to  determine  whether  our  current 
paper  could  be  generalized  to  model  the  boundary  condi¬ 
tions  of  a  structured  phosphor  and  have  agreement  with 
Freed  et  al. 

VI.  CONCLUSION 

This  work  develops  analytical  models  of  OTF,  NPS,  and 
DQE  for  a  turbid  phosphor  irradiated  obliquely.  Our  analysis 
differs  from  much  of  the  previous  work  on  oblique  incidence 
in  that  closed  form  solutions  are  obtained  from  first  princi¬ 
ples,  thereby  providing  greater  insight  into  the  underlying 
detector  physics. 

Ultimately,  the  model  provides  a  platform  for  optimiz¬ 
ing  the  design  of  DBT  detectors.  For  example,  in  a  conven¬ 
tional  front-screen  configuration,  the  model  is  a  useful  tool 
for  optimizing  phosphor  thickness  at  various  angles  of  inci¬ 
dence.  Because  the  incident  angle  is  spatially  variant  across 
the  detector  area,  the  potential  merit  of  designing  a  phos¬ 
phor  with  reduced  thickness  near  its  periphery  has  been 
proposed. 

This  work  demonstrates  that  the  transfer  functions  of  the 
back-screen  have  less  angular  dependence  than  the  front- 
screen,  and  that  high  DQE  is  supported  over  a  broader  range 
of  thicknesses  for  all  incident  angles.  As  a  result,  future 
DBT  detectors  manufactured  on  a  plastic  substrate  instead  of 
glass  should  preferentially  be  operated  in  the  back-screen 
configuration. 
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NOMENCLATURE 

d  =  partial  derivative  operator. 

V2  =  Laplacian  operator. 

E  =  set  membership. 

C  =  Laplace  transform  operator. 

JV'  =  set  of  natural  numbers. 

a  =  polar  angle  of  the  2D  spatial  frequency  vector  v. 

ft  =  a  term  defined  by  Eq.  (21)  to  simplify  intermedi¬ 
ate  calculations. 

y±  =  a  term  defined  by  Eq.  (22)  to  simplify  intermedi¬ 
ate  calculations. 

A  =  nominal  projection  angle. 

5  =  delta  function. 

9  =  angle  of  x-ray  incidence  relative  to  the  normal  to 
the  detector. 

/i  =  x-ray  linear  attenuation  coefficient  of  the 
phosphor. 

v  =  spatial  frequency  vector  with  components  vx  and 

Pj  =  a  quantity  defined  by  Eq.  (13)  that  is  related  to 
surface  reflectivity.  The  subscripts  j  =  0,  1  corre¬ 
spond  to  the  phosphor  backing  and  photocathode, 
respectively. 

a  =  reciprocal  of  the  mean  diffusion  length  of  optical 
photons. 

t  =  inverse  relaxation  length  (the  diffusion  limit 
occurs  with  z  — ►  oo). 

</>( r)  =  product  of  photon  density  and  the  diffusion 
constant. 

i/rk(z)  =  Fourier  transform  of  </>(r)  in  a  plane  of  constant  z. 

Aq  =  x-ray  quantum  detection  efficiency. 

Aq0  =  x-ray  quantum  detection  efficiency  at  normal 
incidence. 

As  =  Swank  information  factor. 

B  =  back- screen  configuration  (often  used  as  a 
subscript). 

Ci,  C2  =  constants  of  integration  used  in  intermediate 
calculations. 

CE  =  contrast-enhanced 

COR  =  center-of-rotation  of  x-ray  tube,  corresponding  to 
the  midpoint  of  the  chest  wall  side  of  the  detector. 

D  =  symbolic  abbreviation  for  detective  quantum 
efficiency. 

DBT  =  digital  breast  tomosynthesis. 

DE  =  dual  energy 


DM  =  digital  mammography. 

DQE  =  detective  quantum  efficiency. 

F  =  front- screen  configuration  (often  used  as  a 
subscript). 

FOV  =  field-of-view. 

G(v)  =  optical  transfer  function  found  by  summing  the 
contributions  of  each  depth  z0  of  the  phosphor. 

G(v,  z0)  =  optical  transfer  function  associated  with  the  depth 
0  <  <  T  of  the  phosphor. 

h  =  source-to-COR  distance  for  rotating  x-ray  tube. 
i  =  imaginary  unit  given  as 

j(z)  =  photon  currents  across  any  plane  of  constant  z. 
k  =  a  quantity  equivalent  to  27iv. 

MTF  =  modulation  transfer  function. 

n  =  iteration  number  for  Newton’s  Method. 

NPS  =  noise  power  spectra. 

N(z0)  =  relative  x-ray  signal  at  the  depth  z0  of  the 
phosphor. 

OLED  =  organic  light-emitting  diode. 

OTF  =  optical  transfer  function. 

p  =  independent  variable  of  the  Laplace  transform  of 
a  function. 

PSF  =  point  spread  function. 

q  =  a  quantity  defined  in  Eq.  (6)  to  simplify  interme¬ 
diate  calculations. 

QDE  =  quantum  detection  efficiency. 

Tj  =  reflectivity  of  a  surface,  where  j  =  0,  1  correspond 
to  the  phosphor  backing  and  photocathode, 
respectively. 

*c(v)  =  Lubberts  fraction. 

Rn(v)  =  r^tio  of  the  quantum  noise  power  to  the  total 
noise  power. 

S( r)  =  source  function,  modeled  as  point-like. 

T  =  phosphor  thickness. 

TFT  =  thin-film  transistor. 

W(y)  =  noise  power  spectra. 

x  =  position  along  the  chest  wall  side  of  the  detector. 
y  =  position  perpendicular  to  the  chest  wall  side  of 
the  detector. 

zo  =  depth  within  a  phosphor  (with  or  without  the 
subscript). 
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ABSTRACT 

45  Purpose:  Digital  breast  tomosynthesis  (DBT)  is  a  3D  x-ray  imaging  modality  in  which  tomographic 
sections  of  the  breast  are  generated  from  a  limited  range  of  tube  angles.  Because  oblique  x-ray 
incidence  shifts  the  image  of  an  object  in  sub-pixel  detector  element  increments  with  each 
increasing  projection  angle,  it  is  demonstrated  that  DBT  is  capable  of  super-resolution  ( i.e .  sub¬ 
pixel  resolution). 

50  Methods:  By  convention,  DBT  reconstructions  are  performed  on  planes  parallel  to  the  breast 

support  at  various  depths  of  the  breast  volume.  In  order  for  resolution  in  each  reconstructed  slice  to 
be  comparable  to  the  detector,  the  pixel  size  should  match  that  of  the  detector  elements;  hence,  the 
highest  frequency  that  can  be  resolved  in  the  plane  of  reconstruction  is  the  alias  frequency  of  the 
detector.  This  study  considers  reconstruction  grids  with  much  smaller  pixilation  to  visualize  higher 

55  frequencies.  For  analytical  proof  of  super-resolution,  a  theoretical  framework  is  developed  in  which 
the  reconstruction  of  a  high  frequency  sinusoidal  input  is  calculated  using  both  simple 
backprojection  (SBP)  and  filtered  backprojection  (FBP).  To  study  the  frequency  spectrum  of  the 
reconstruction,  its  Fourier  transform  is  also  calculated.  The  experimental  feasibility  of  super¬ 
resolution  was  investigated  using  a  commercial  DBT  system  (Selenia  Dimensions,  Hologic, 

60  Bedford,  MA)  and  a  commercial  prototype  reconstruction  solution  (Briona™,  Real  Time 

Tomography,  Villanova,  PA).  Images  of  bar  patterns  with  frequencies  higher  than  the  detector  alias 
frequency  were  acquired  and  subsequently  reconstructed.  To  demonstrate  the  clinical  importance  of 
super-resolution,  the  visibility  of  microcalcifications  was  compared  among  reconstruction  grids  of 
varying  sizes. 

65  Results:  Using  analytical  modeling,  it  is  shown  that  the  central  projection  cannot  resolve 

frequencies  exceeding  the  detector  alias  frequency.  The  Fourier  transform  of  the  central  projection 
is  maximized  at  a  lower  frequency  than  the  input  as  evidence  of  aliasing.  By  contrast,  SBP 
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reconstruction  can  resolve  the  input,  and  its  Fourier  transform  is  correctly  maximized  at  the  input 
frequency.  Incorporating  filters  into  the  reconstruction  smoothens  pixilation  artifacts  in  the  spatial 
70  domain  and  reduces  spectral  leakage  in  the  Fourier  domain.  In  addition,  it  is  demonstrated  that  the 
feasibility  of  super-resolution  is  dependent  on  the  directionality  of  the  input  frequency.  Super¬ 
resolution  is  achievable  over  a  broad  range  of  positions  for  frequencies  parallel  to  the  chest  wall 
side  of  the  breast  support,  but  is  feasible  only  at  certain  positions  for  frequencies  perpendicular  to 
the  chest  wall.  At  the  latter  orientation  of  the  input  frequency,  super-resolution  is  achievable  at 
75  positions  sufficiently  displaced  from:  (1)  the  chest  wall,  and  (2)  the  mid  plane  perpendicular  to  the 
chest  wall  and  to  the  breast  support.  Consistent  with  the  analytical  results,  experimental 
reconstructions  of  bar  patterns  showed  visibility  of  frequencies  greater  than  the  detector  alias 
frequency.  For  the  two  orientations  of  the  bar  patterns,  super-resolution  was  present  at  positions 
predicted  from  analytical  modeling.  Ultimately,  super-resolution  was  found  to  improve  the 
80  visibility  of  fine  morphological  details  of  microcalcifications  in  select  examples  of  clinical  DBT 
images. 

Conclusions:  This  work  demonstrates  the  feasibility  of  super-resolution  in  DBT,  and  has 
applications  in  the  visualization  of  microcalcifications  and  other  subtle  signs  of  breast  cancer. 

85  Key  words:  Digital  breast  tomosynthesis  (DBT),  super-resolution,  bar  pattern  phantom, 
microcalcifications,  filtered  backprojection  (FBP). 
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I.  INTRODUCTION 

Digital  breast  tomosynthesis  (DBT)  is  a  3D  imaging  modality  in  which  low  dose  x-ray 
projections  are  acquired  over  a  limited  angular  range  about  the  breast.  Using  digital  image 
95  reconstruction  techniques,  tomographic  sections  at  all  depths  of  the  breast  volume  are  subsequently 
generated.  Unlike  2D  digital  mammography  (DM),  DBT  can  filter  out  overlapping  anatomical 
structures  which  may  hide  a  tumor.  Early  clinical  trials  indicate  that  DBT  has  greater  sensitivity 
and  specificity  for  early  cancer  detection  relative  to  DM.1 

In  conventional  practice,  the  reconstructed  slices  are  generated  on  planes  parallel  to  the 
100  breast  support.  In  order  to  have  the  same  in-plane  resolution  in  the  reconstruction  as  the  detector, 
the  pixel  size  in  each  reconstructed  slice  should  match  that  of  the  detector  elements.  Using  this 
approach,  the  highest  frequency  that  can  be  resolved  in  the  plane  of  reconstruction  is  the  alias 
frequency  of  the  detector.  This  study  considers  the  possibility  for  reconstruction  grids  with  much 
smaller  pixilation  so  that  higher  frequencies  can  be  visualized.  Because  non-normal  x-ray  incidence 
105  causes  the  image  of  an  object  to  be  translated  in  sub-pixel  detector  element  increments  with  each 
increasing  projection  angle,  it  is  demonstrated  in  this  work  that  DBT  is  capable  of  super-resolution 
( i.e .,  sub-pixel  resolution). 

Super-resolution  has  been  well-described  in  a  number  of  applications  involving 
reconstruction  from  projections,2  including  forensics,  satellite  imaging,  computed  tomography  (CT), 
110  and  magnetic  resonance  imaging  (MRI);  however,  to  our  knowledge,  its  potential  in  DBT  has  not 
yet  been  explored.  An  understanding  of  super-resolution  and  an  analysis  of  how  to  optimize  its 
presence  may  prove  to  be  useful  for  designing  the  highest  quality  DBT  systems.  Although  it  is 
possible  to  improve  spatial  resolution  simply  by  reducing  the  pixel  size  of  the  detector,  there  are 
practical  lower  limits  on  the  sizes  that  can  be  manufactured.  In  addition,  one  drawback  of  reducing 
115  the  pixel  size  is  decreasing  the  mean  number  of  photons  incident  on  each  detector  element  and 


4 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


hence  decreasing  the  signal-to-noise  ratio  (SNR)  per  pixel  according  to  Poisson  statistics  for  x-ray 
distributions.  Clinically,  super-resolution  should  be  beneficial  to  diagnostic  radiologists  by 
improving  the  visibility  of  microcalcifications  and  other  subtle  signs  of  breast  cancer  with  no 
increased  radiation  dose  to  the  patient. 

120  In  this  study,  a  theoretical  framework  for  investigating  super-resolution  in  DBT  is  developed 

by  calculating  the  reconstruction  of  a  sine  input  whose  frequency  is  greater  than  the  alias  frequency 
of  the  detector.  For  optimal  visualization  of  high  frequencies  in  the  3D  image,  an  infinitesimally 
fine  (i.e.,  non-pixilated)  reconstruction  grid  is  considered.  The  reconstruction  techniques  include 
both  simple  backprojection  (SBP)  and  filtered  backprojection  (FBP).  In  order  to  investigate  the 
125  experimental  feasibility  of  super-resolution  using  a  commercial  DBT  system,  images  of  a  bar 

pattern  phantom  with  frequencies  higher  than  the  alias  frequency  of  the  detector  were  acquired  and 
subsequently  reconstructed.  In  addition,  clinical  images  of  microcalcifications  were  analyzed,  and 
differences  in  the  visibility  of  microcalcifications  were  examined  on  various  reconstruction  grid 
sizes. 

130 

II.  METHODS 

II.  A.  Input  Object  and  Acquisition  Geometry 

An  analytical  framework  for  investigating  the  potential  for  super-resolution  in  DBT  is  now 
developed  by  calculating  the  reconstruction  of  a  high  frequency  sinusoidal  input.  Accordingly, 

135  suppose  that  a  rectangular  prism  with  infinite  extent  in  the  x  and  y  directions  has  a  linear  attenuation 
coefficient  p(x,  y,  z)  which  varies  sinusoidally  along  the  x  direction  with  frequency  [).  Throughout 
the  remainder  of  this  manuscript,  the  input  object  will  be  termed  a  “sine  plate”.  With  the  xz  plane 
defining  the  chest  wall,  the  frequency  vector  is  therefore  oriented  parallel  to  the  chest  wall  side  of 
the  breast  support.  Figure  1  illustrates  a  cross  section  of  the  sine  plate  in  the  xz  plane.  As  shown, 
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140  the  rectangular  prism  is  positioned  between  z  =  zo  +  ell  and  z  =  zo  -  ell,  where  zo  is  the  central 

height  of  the  prism  and  e  is  the  prism’s  thickness.  Defining  the  origin  O  as  the  midpoint  of  the  chest 
wall  side  of  the  detector,  the  attenuation  coefficient  may  be  written  as 


H(x,y,z) 


C  •  cos  \ln  /0(x  -  x0)]  •  rect 


(1) 


where  C  is  a  constant  denoting  the  amplitude  of  the  waveform,  xo  is  a  translational  shift  in  the 
waveform  relative  to  the  origin,  and  the  rect  function  is  defined  by  the  expression 


rect(w) 


1  ,  |  w  |<  1  /  2 

0  ,\u\>\/l' 


By  setting  the  amplitude  C  to  1/e,  ju(x,y,  z)  may  be  normalized4  so  that  the  total  attenuation  found 


by  integrating  along  the  z  direction  is  simply  cos[27r/o(x  -  xo)]  for  all  e.  Provided  that  |z  -  zo|  <  e,  the 


ID  Fourier  transform  (IF\ )  of  Eq.  (1)  along  the  x  direction  peaks  at  the  frequencies  fx  =  ±/0,  and 


150  vanishes  at  all  other  frequencies.5 

,  y,  z)  =  f  y,  z )  •  e~2*ifxXdx  (3) 

J  —oo 

=  £[e-^.Wx  _/„)  +  S(f,  +  /„)]•  rectf  ^£.1  (4) 

Typically,  only  the  positive  frequency  fx  =  +/o  is  of  interest  in  a  physical  measurement.  Thus, 
although  it  is  non-physical  for  an  attenuation  coefficient  to  vary  between  negative  and  positive 
155  values,  formulating  ju(x,  y,  z)  by  Eq.  (1)  is  helpful  for  a  thought  experiment  in  interrogating  the 

reconstruction  of  a  single  input  frequency.  An  analysis  of  the  case  for  which  the  input  frequency  is 
oriented  along  the  y  direction  (i.e.,  perpendicular  to  the  chest  wall)  is  considered  separately  in 
Appendix  A. 

The  most  general  DBT  acquisition  geometry  with  a  divergent  x-ray  beam  and  a  rotating 
160  detector  is  now  modeled.  In  acquiring  the  nth  projection,  the  focal  spot  emits  x  rays  at  the  nominal 
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projection  angle  y/n  relative  to  the  center-of-rotation  (COR)  of  the  DBT  system.  The  COR  and  the 
focal  spot  lie  in  the  plane  of  the  chest  wall.  In  addition,  the  detector  rotates  about  the  y  axis  at  the 
angle  yn  relative  to  the  x  direction.  The  two  parameters  y/n  and  yn  are  determined  from  the  nominal 
angular  spacing  Ay/  and  the  detector  gear  ratio  g  by  the  relations 


165 

ll 
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(5) 
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1  P 

(6) 
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For  an  odd  number  of  A  total  projections,  the  index  n  varies  between  -( N -  l)/2  and  ( N -  l)/2,  and 
the  special  case  n  =  0  defines  the  central  projection. 

As  a  final  step  in  this  section,  it  is  useful  to  calculate  the  incident  angle  at  each  point  on  the 
170  detector.  Following  Figure  2,  the  vector  from  O  to  an  arbitrary  point  C  on  the  detector  for  the  nth 
projection  is 
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OC  =  uli'n+u2l  (7) 

=  («!  cos  yn  )i  +  u2j  +  ( Wj  sin  yn  )k .  (8) 

The  matrix  transformation  between  the  primed  and  unprimed  coordinate  systems  supports  the 
transition  from  Eq.  (7)  to  Eq.  (8). 


rq 

'  cosr„ 

0 

siny„" 

rn 

t 

= 

0 

1 

0 

j 

(9) 

kJ 

v-sin/„ 

0 

cos 

Additional  vectors  from  O  to  the  COR  at  point  B  and  from  the  COR  to  the  focal  spot  at  point  A  are 


r* 

II 

0Q 

O 

(10) 

BA  =  (-h  sin y/n )i  +  ( h cos y/n) k  , 

(11) 

180  where  /  is  the  distance  between  O  and  the  COR.  In  Eq.  (1 1),  it  is  assumed  that  for  positive  values  of 
y/n ,  the  x  component  of  the  displacement  of  the  focal  spot  at  A  relative  to  the  COR  at  B  is  negative. 
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This  sign  convention  is  chosen  so  that  positive  values  of  y/„  cause  the  x  component  of  the  trajectory 
from  A  to  C  to  be  positive  for  positive  values  of  u\  (Figure  2).  By  the  summation  rules  for  vectors, 
the  net  vector  from  point  C  on  the  detector  to  the  focal  spot  at  A  is 

185  CA  =  -OC  +  OB  +  BA  (12) 

=  — («j  cosy„  +  hsmi//n)i-u2\  +  (l  +  hcosi//n  -  ut  sinyH)k  .  (13) 

Thus  the  angle  of  incidence  is  found  from  the  expression 


COS0n 


cA.k; 
CA  k: 


(14) 


giving 


6„  =  arccos 


_ h  cos(^„  ~y„)  + 1  cos  yn _ 

/(«,  cos/,,  +  hsmiyn)2  +  u;  +(/  +  /?cos^;(  -u,  siny(!f 


The  dot  product  in  Eq.  (14)  has  been  computed  using  Eq.  (9)  to  write  k',  in  terms  of  the  unprimed 


unit  vectors. 


II.  B.  Detector  Signal 

195  To  calculate  the  detector  signal  for  each  projection,  it  is  useful  to  perform  ray  tracing 

through  the  input  object.  We  begin  by  defining  the  line  from  the  focal  spot  at  A  to  the  incident 
point  on  the  detector  at  C  for  the  nth  projection.  This  line  can  be  expressed  as  the  parametric 
equation 

xj  f  Uj  cosy„  +hsini//rt  ^  (  -h sin^„  N 

y  =w  u2  +  0  ,  (16) 

z)  lvM1siny„ -l-hcosy/nJ  [^l  +  h  cosy/nJ 

200  where  (x,  y,  z)  is  a  point  in  M3  and  w  is  a  free  parameter.  The  focal  spot  at  A  has  been  defined  to 
correspond  with  w  =  0,  while  the  incident  point  at  C  has  been  defined  to  correspond  with  w  =  1 . 
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The  x-ray  path  length  C„  through  the  input  object  for  the  nth  projection  image  is  determined  from 


the  intersection  of  Eq.  (16)  with  the  planes  z  =  zo  +  el 2  and  z  =  zo  -  e/2.  The  values  of  w  for  these 
two  points  are 


205 


+  _  z0  +  (£■  /  2)  —  /  —  /z  COS  \j/n 
Mj  sin  yn-l  —  h  cos  y/n 


(17) 


z0  -  (e  /  2)  - 1  -  h  cos  y/n 
ulsmyn  -  l -hcosy/n 


(18) 


where  wf  and  wn  correspond  to  the  entrance  and  exit  points  of  the  x-ray  through  the  input, 
respectively.  For  the  nth  projection  image,  total  x-ray  attenuation  Aju(n)  is  now  found  by  integrating 


ju(x,  y,  z)  along  £„. 

210  Aju(n)  =  J  juds  (19) 

The  differential  arc  length  ds  along  C„  is 


if  dx 
ds  =  .  — 

VI  dw 


fdyf] 

2 

(dz\ 

+ 

) 

v  dw  j 

+ 

KdW  j 

dw 


(20) 


=  (i/j cos yn  +hs>imf/n)2  +u2  +(l  +  hcost//n  -u^myj2  dw  (21) 

=  [h  cos (y/n  -  yn )  +  /  cos  yn  ]  sec(0„  )-dw.  (22) 

215  Eq.  (22)  follows  from  Eq.  (15).  Substituting  Eq.  (22)  into  Eq.  (19)  yields  the  total  x-ray  attenuation 

Aju(n)  =  Kn^  "  cos[2^-/0(m1cos7h  +  hsim//n)w  +  An\dw  (23) 

_  Kn  (Sin  [2^/o  (M1  C0S  Yn  +  h  Sin  Wn  )' Wn  +  K  ]  “  Sil1  [2*/o  (W1  C0S  Yn  +  k  Sin  Yn  )W+n  +  K  ]) 

2^/o(“i  C0SYn  +hsiny/J 


(24) 
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where 


220 


225 


230 


Kn=C[h cos(^„  -  Yn )  + 1  cos 7„ ] sec 9n 
K  =-2xf0(hsim//n+x0). 

Using  a  sum-to-product  trigonometric  identity  for  real  numbers  b\  and  £>2 
sin(hj )  -  sin(62 )  =  2  cos  [ 
one  may  rewrite  Eq.  (24)  as 


(25) 

(26) 


6  h  +  b2 

' bx  -b2^ 

sin 

l  2  J 

{  2  J 

(27) 


A/J(n)  =  Kn  (vU  -  w+n )  cos  [>r/0 (n,  cos  +  h  sin  ^ )«  +  w; )  +  Xn  ]  sine  [/0 (ut  cos  +  A  sin  y/n  ){w~  -  w+n 

(28) 


£K.  COS 


2nf0(l  +  h  cos  y/n  -  z0  ){ux  cos  yn+h  sin  y/n )  , 

emp 

sfM  cos  y„  +  h  sin  y/n ) 

1  An 

l  +  h  cos  y/n  —  Mj  sin  yn 

/  +  hcosYn  -Mj sin yn 

l  +  hcosy/n  -Mj  sinyH 


(29) 


where 


.  .  .  sin (ku) 

smc(u)  - 


7VU 


(30) 


The  transition  from  Eq.  (28)  to  (29)  follows  from  Eqs.  (17)  and  (18).  Eq.  (29)  possesses  a 
singularity  at  ux  =  (/  +  h  cos  y/n )  esc  yn ,  the  point  at  which  the  denominator  vanishes.  For  typical 
acquisition  geometries,  this  singularity  is  not  expected  to  correspond  to  a  position  on  the  detector, 
since  neither  the  attenuation  coefficient  ju(x,  y,  z )  nor  the  path  length  C„  should  have  an  infinity. 


235  Eq.  (29)  provides  an  expression  for  signal  intensity  versus  position  along  the  detector, 

assuming  that  the  detector  is  non-pixilated  and  possesses  an  x-ray  converter  whose  modulation 
transfer  function  (MTF)  is  unity  at  all  frequencies.  An  amorphous  selenium  (a-Se)  photoconductor 
operated  in  drift  mode  is  a  good  approximation  for  an  x-ray  converter  with  these  properties.6  In  a 
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clinical  setting,  a-Se  is  placed  in  contact  with  a  plate  of  amorphous  silicon  (a-Si)  in  which  a  thin- 
film  transistor  (TFT)  array  samples  detector  signal  in  pixels  (i.e.,  detector  elements).7'9  The 
logarithmically-transformed  signal  in  the  mth  detector  element  for  the  nth  projection  is 


Vju(m 


ay(my+ 1)  rax(mx  + 1/2)  du,  dll-, 

A/u{n) — 1 — 2-. 

iymy  J  ax(mx-V2)  a *  Qy 


In  Eq.  (31),  mx  and  my  are  integers  used  for  labeling  detector  elements,  and  ax  and  ay  denote  detector 
element  lengths  in  the  directions  parallel  and  perpendicular  to  the  chest  wall,  respectively.  In  the 
special  case  of  square  detector  elements,  it  is  assumed  that  ax  =  ay  =  a.  Detector  elements  are 
centered  on  u\  =  mxax  and  112  =  (my  +  1/2 )ay,  where  mx  e  Z  and  m  e  Z* . 

It  is  important  to  note  that  the  integrand  in  Eq.  (3 1)  is  dependent  on  both  u\  and  112  due  to  the 
dependency  of  Kn  [Eq.  (25)]  on  the  incident  angle  6n  [Eq.  (15)].  However,  because  6n  should  vary 
minimally  within  the  area  of  a  single  detector  element,  total  attenuation  can  be  well  approximated 
by  the  expression 

Aju(n)  =  Aju(n)\  ,  (32) 

un  umn 

where  dmn  is  the  evaluation  of  6n  at  the  centroid  of  the  mth  detector  element 

a  =  f)  \  1331 

111/7  n \(ui,u1)=(mKax,[my+\l2]ay} ?  ' 

so  that 


^  .  .  rax(mx+ 1/2)  ~  du . 

D//(m,n)=|  ^(«) - 1 

Jax(mx- 1/2)  q 


Because  it  would  be  difficult  to  evaluate  Eq.  (34)  in  closed  form,  it  is  appropriate  to  apply 
approximate  integration  techniques.  One  such  method  is  the  midpoint  formula10 


V/u(m  ,«)=  lim—  Y  Aju(jx,n)., 

J^cc  Jr  77, 


where 
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260 


AjU(jx,n)  =  Aju(n) 


,  1  1  ' 
U]=aJ  — - 1-  mY  — 

JY  2 


(36) 


The  raw  signal  S/i(u\  ,ih)  across  the  detector  can  now  be  determined  for  the  nth  projection  as 


Sju(uvu2 )  =  y^Z>//(m,n)  •  rect 


■rect 


V  “x  J 


(37) 


Using  this  expression  for  raw  signal,  it  is  now  possible  to  calculate  the  x-ray  transform11  X/i{t\  ,C). 


Xju(t j  ,t2)  =  X!  ^//(m,  n)  •  rect 


^  t,  sec6>m„  -mar  ^ 

1  mn  x  x  |  red 


V 


a. 


rt2^c0m m-(mv+\/2)a^  N 


) 


(38) 


265  To  justify  the  transition  from  Eq.  (37)  to  Eq.  (38),  one  must  determine  the  affine  parameters  t\  and  ti 
in  terms  of  u\  and  112  by  considering  a  line  segment  OD  which  is  orthogonal  to  AC  and  which 
connects  the  origin  with  the  x-ray  beam  (Figure  2).  From  trigonometry,  the  length  |t|  =  yjt }  +  ^  of 


OD  is 


|t|  =  |u|cos#„ .  (39) 

270  By  generalizing  Eq.  (39)  to  components,  one  finds  t\  =  u\cos0n  and  h  =  i/2cos61„.  In  Eq.  (38),  the 
incident  angle  across  the  area  of  the  mth  detector  element  for  the  nth  projection  has  been 
approximated  by  its  value  at  the  centroid. 


II.  C.  Filtered  Backprojection  (FBP)  Reconstruction  from  the  Projections 

275  The  reconstructed  attenuation  coefficient  can  now  be  determined  by  filtering  the  x-ray 

transform  with  the  function  <jit  \,  to)  and  backprojecting  the  result  along  the  ray  of  incidence.  It  is 
customary  to  apply  filtering  exclusively  to  frequencies  within  the  plane  of  the  x-ray  tube  motion,  so 
that  the  filter’s  2D  Fourier  transform  Jr1(f)(  f\,  f1)  is  independent  of/2 
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•Wp/2)  =  Wi).  (4°) 

280  and  hence 

Wi  ,t2)  =  <!>(h)8(t2) .  (41) 


The  specific  formula  for  <fi(t\)  will  be  addressed  in  the  next  section.  Assuming  that  the 
reconstruction  grid  is  infinitesimally  fine  ( i.e .,  non-pixilated),  the  filtered  backprojection  (FBP) 
reconstruction  is 


285 


/^FBP 


Vju(m,n) 

N 


)  *  rect 


r  ^  Sec  “  mxa:  ' 


av 


t\  =x'n  COS  0mn  +z’n  cos( rmw )  sin(  0mn ) 


rect 


^  t2sec0mn-(m+\/2)a,  A 


(42) 


t2  =y'n  cos  6mn  +z'n  sin(rmw )  sin(  0mn ) 


where  //fbp  is  the  reconstructed  attenuation  coefficient  and  *  is  the  convolution  operator.  Within 
the  plane  of  the  detector,  backprojection  of  signal  in  the  mth  detector  element  for  the  nth  projection  is 


directed  azimuthally  along  the  angle  rm„  relative  to  the  in'  axis  (Figure  2).  As  shown  in  Eq.  (42), 

backprojection  may  be  performed  for  each  of  the  N  projections  using  the  primed  coordinate  system. 
290  To  evaluate  Eq.  (42)  at  the  point  (x,  y,  z )  in  the  unprimed  coordinate  system,  one  applies  the  matrix 
transformation  given  in  Eq.  (9). 

It  is  now  important  to  illustrate  how  the  azimuthal  backprojection  angle  Tm„  is  calculated. 
Begin  by  considering  an  arbitrary  point  E  along  the  x-ray  beam  at  which  signal  is  backprojected 


from  the  incident  point  C  (Figure  2).  A  line  segment  along  the  k;/  direction  may  then  be  drawn 


295  from  E  to  the  point  F  on  the  detector  for  the  nth  projection.  As  a  result,  within  the  plane  of  the 


detector,  backprojection  is  directed  from  C  to  F  at  the  angle  Tn  relative  to  the  in'  axis.  Point  G  may 
now  be  defined  as  the  position  at  the  chest  wall  side  of  the  detector  which  is  collinear  with  points  C 
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and  F.  A  derivation  of  the  formula  for  Yn  requires  knowledge  of  the  distance  dn  between  G  and  O, 


which  is  now  calculated. 

GO  =  dni'n 

(43) 

=  (dn  cos  y„ )i  +  (dn  sin  )k 

(44) 

Since  ACG  and  ECF  are  similar  triangles,  GA  is  parallel  to  FE  and  is  in  turn  parallel  to  k/  . 
Denoting  x  as  the  cross  product  operator,  it  follows  that 

GAxk/=0,  (45) 

305  where 

GA  =  GO  +  OA  (46) 

=  (dn  cos  yn  -  h  sin  y/n  )i  +  (dn  sin  yn+l  +  h  cos  y/n  )k .  (47) 


To  calculate  OA  in  Eq.  (46),  Eqs.  (10)  and  (1 1)  have  been  summed.  Substituting  Eqs.  (9)  and  (47) 
into  the  cross  product  of  Eq.  (45)  gives 


310 


GAxk/ 


i 

dn  cos  yn~h  sin  y/n 

-sin  Yn 


j  k 

0  dn  smy„  +1 -\-hcosi//n 
0  cos  Yn 


(48) 


=  ~[d«+l sin Yn~h sin (x//n  -  yn )]  j  .  (49) 

By  combining  Eqs.  (45)  and  (49),  one  can  solve  for  dn. 

dn=h  sm(i//n  -  Yn )  -  /  sin  yn  (50) 

Using  this  result,  it  follows  from  trigonometry  that 

315  cosT  =  .  (51) 

^jdjw2 


sinr„=  .  2 

y(wi  +d„)2  +u22 


(52) 
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Substituting  the  coordinates  of  the  detector  element  centroid  into  Eqs.  (51)-(52),  one  finds  that  the 
azimuthal  backprojection  angle  for  the  mth  detector  element  in  the  nh  projection  satisfies  the 
properties 


320 


cos  r . = 


m  ar  +d„ 


\J(mxax  +  dnf  +  (mv  +  1  Ufa2 


(53) 


sinf 


m n 


_ (my  +  1  /  2)Qj, _ 

\j(mxax  +  dnf  +  (mv  + 1  /2)2aJ 


(54) 


These  relations  are  the  expressions  needed  for  FBP  reconstruction  in  Eq.  (42).  One  special  case  of 
Eq.  (42)  is  simple  backprojection  (SBP)  reconstruction 


»)=E— •"  ct 


N 


x'  +  z'n  cos(Tm„ )  tan(6>m„ )  -  mxax 


a „ 


•rect 


ry'n+Zn  sin(rm„  )  ~  (m  +  1  /  2)fl  ^ 


(55) 


325  where  B  denotes  the  backprojection  operator.  With  SBP,  the  filter  <f>(tvt2)  effectively  becomes  the 


product  5{tx)5(t2) .  According  to  Eq.  (55),  backprojection  in  the  primed  coordinate  system  occurs 
by  translating  x'n  and  y'  by  -z'  cos(rmn)tan(6»m„)  and  -z'  sin(  r,n;i )  tan(^m;i ) ,  respectively,  where 
z'  is  the  height  of  the  backprojected  point  (E)  above  the  plane  of  the  detector.  These  translational 
shifts  are  illustrated  in  Figure  2. 

II.  D.  Formulation  of  the  Reconstruction  Filter 

1 3 

Following  Zhao’s  linear  systems  theory  for  DBT,  a  ramp  (RA)  filter  should  be  applied  to 
the  x-ray  transform  of  each  projection  to  reduce  the  low  frequency  detector  response.14  The  filter  is 
truncated  at  the  spatial  frequencies  f\  =  -c  and  f\  =  +c  in  the  Fourier  domain. 
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335 


M 

[0  ,\fy>Z 


(56) 


The  spatial  representation  ^ra(6)  of  the  RA  filter  is  determined  by  its  inverse  Fourier  transform.5 


4A(h)  =  f^4A(/i)-eWl^ 


(57) 


=  4:2[2sinc(24:t1)-  sinc2(#,)] 


(58) 


Using  this  result,  the  convolution  in  Eq.  (42)  can  be  calculated. 

ax  cos(6»mn )  [cos^a^  cos  0mn )  cos  [2^(f,  -  mxax  cos  Omn )]  - 1] 
+2 (h  -  m xax  cos  0m )  si njxaj  cos  6>mn )  sin  [2 gg,  -  cos  6^ )]_ 
2^2  [h  -  K  - 1  /  2)«,  COS  <9mn  ]  [f,  -  (mv  + 1/  2)a,  cos  Omn  ] 

(59) 

Since  noise  tends  to  occur  at  high  frequencies,  a  spectrum  apodization  (SA)  filter  is  often  applied  in 
addition  to  the  RA  filter  in  order  to  reduce  the  high  frequency  detector  response.  Following  Zhao’s 
approach,  a  Hanning  window  function  is  the  SA  filter. 


340  ^M0*rect 


a 


345 


^sa(./D  = 


1  +  COS 


V  4 


,l/i££ 

,|/il>£ 


(60) 


4a  (0 


^sinc(2^() 


l-4^2 


(61) 


According  to  the  convolution  theorem,  the  net  filter  is  thus 


The  convolution  of  the  net  filter  in  Eq.  (62)  with  the  rect  function  in  Eq.  (42)  can  be  performed  in 
350  closed  form  similar  to  Eq.  (59).  This  expression  is  omitted  as  it  is  lengthy. 


16 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


II.  E.  Fourier  Transform  of  the  DBT  Images 

According  to  Eq.  (4),  the  Fourier  transform  of  the  input  along  the  x  direction  peaks  at  the 
frequencies  fx  =  ±/0.  To  determine  whether  the  frequency  spectra  of  the  DBT  images  possess  this 
expected  dependency  on  fo,  their  continuous  Fourier  transforms  may  be  considered.  Within  the 
plane  of  the  detector,  the  2D  Fourier  transform  of  the  nb  projection  is 


T2(Sju)(f,,f2)  =  f”  P  Sp(uvu2)-e^m{f'u'+f^}dii,dii2 

J  -OO  J  -00 


=  «  A 


,  •  sinc(ax/1  )sinc (ayf2 )  •  ^  Vju(m,  i 


-2 7d\  mxaxfl+(my+\/2)ayf2 


360  The  2D  Fourier  transform  of  the  reconstruction  along  the  x  andy  directions  may  now  be  calculated 
by  considering  a  fixed  height  z.  Because  this  study  only  considers  input  frequencies  parallel  to  the 
xy  plane,  it  is  unnecessary  to  transform  along  the  z  direction.  Although  filtered  backprojection 
reconstruction  in  Eq.  (42)  is  performed  in  the  primed  coordinate  system,  it  is  important  to  take  the 
Fourier  transform  in  the  unprimed  coordinated  system.  As  such,  the  reconstructed  attenuation 
365  coefficient  can  be  written  in  the  form 


/w  T,  z) = ■  [  Px  {tx )][ 


■  [/^2  (^2  )]  tr 


f2-cr3m«x+cr4in«T+cr5m«z 


where 


px(tx)  =  ) * rect  -  SCC  — — 

l  a. 


p2{t2)  =  rect 


t2sec0mH-(mv  +  l/2)ay 


370  and 


^ \  m;;  =  cos(y„ )  cos(<9m„ )  -  cos(Tmn )  sin(y„ )  sin(6»mn ) 


^2  mn  =  sin  (jn )  cos((9mn )  +  cos(Tm„ )  cos  (yn )  sin((9m„ ) 
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375 


390 


^3m„  =  -  sin(rm„ )  sm(yn )  sin(£m„ ) 


^4mn  ~  COS  @mn 

=  sin(rm„)cos(y„)sin(^m„) . 

According  to  the  convolution  theorem,  the  Fourier  transform  of  Eq.  (66)  under  the  frequency 
variable  f\  is 


(70) 


(71) 


(72) 


•^A(/i)  =  ^(/i)  •  ax  cos^nJsincK/l  cos6>m„)e 


-2  7rimxaxfl  cos  6m 


(73) 


In  the  special  case  of  SBP  reconstruction,  the  filter  in  Eq.  (73)  is  unity.  In  a  similar  fashion,  the 
380  Fourier  transform  of  Eq.  (67)  may  be  written 


APi(fi)  =  ay  C0S(6L)sinc (a  f2  cos0mn)e 


-2n:i{my+\l  2)ayf2  cos  6^ 


(74) 


The  2D  Fourier  transform  of  Eq.  (65)  at  the  fixed  depth  z  is  now  determined  from  the  expression 


N 

where  Iymn(x)  is  given  by  the  integral 
385  /vm„  W  =  Pi (AmA  +  +  Am,,2) '  e2Kfy*dy . 


(76) 


To  evaluate  Eq.  (76),  one  can  make  the  change  of  variables  rj  =  cr3mnx  +  cr4mny  +  oimnz .  Since 
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where  Jrlp1  has  been  previously  calculated  in  Eq.  (74).  Using  Eq.  (79),  Eq.  (75)  can  now  be 
rewritten  as 


395 
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To  evaluate  Eq.  (81),  it  is  helpful  to  perform  the  substitution  rjxmn  =  <xlmnx  +  crlmnz  . 
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is  given  by  Eq.  (73). 
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400  where  is  given  by  Eq.  (73).  The  final  expression  for  the  2D  Fourier  transform  of  the 
reconstruction  can  now  be  derived  by  combining  Eqs.  (80)  and  (84). 
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405  A  special  case  of  this  result  is  important  to  consider. 
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Eq.  (86)  is  useful  for  analyzing  the  reconstruction  of  an  input  frequency  oriented  along  the  x 
direction;  that  is,fy  =  0. 


410  III.  RESULTS 

III.  A.  Input  Frequency  Directed  Parallel  to  the  Chest  Wall  Side  of  the  Breast  Support 

Image  acquisition  is  now  simulated  for  a  Selenia  Dimensions  integrated  multi-mode 
mammography  and  tomosynthesis  x-ray  system  (Hologic  Inc.,  Bedford,  MA)  having  15  projections, 
an  angular  spacing  (At//)  of  1 .07°  between  projections,  a  source-to-COR  distance  (h)  of  70.0  cm,  a 
415  COR-to-origin  distance  (/)  of  0  cm,  and  square  detector  element  length  (ax  =  ay  =  a)  of  140  pm.  In 
addition,  the  sine  plate  has  a  thickness  (e)  of  0.5  mm,  a  translational  shift  (xo)  of  0  mm  along  the 
direction  of  the  chest  wall  side  of  the  breast  support,  and  a  frequency  (/o)  of  0.7a'1  (5.00  lp/mm) 
parallel  to  the  x  axis.  To  illustrate  DBT’s  potential  for  super-resolution,  the  input  frequency  is 
specified  to  be  higher  than  the  detector  alias  frequency  0.5a'1  (3.57  lp/mm).  The  sine  plate  is  placed 
420  at  a  height  corresponding  to  the  mid-thickness  of  a  typical  breast  size  (50.0  mm  thick)  under 

compression.  With  the  breast  support  positioned  25.0  mm  above  the  origin  of  the  detector,  the  sine 
plate  is  therefore  positioned  at  the  height  zo  =  50.0  mm. 

FBP  reconstructions  are  performed  with  either  the  RA  filter  alone  or  the  RA  and  SA  filters 
together,  assuming  a  truncation  frequency  (<f)  of  2a'1  (14.3  lp/mm).  Although  <f  is  typically  chosen 
425  to  be  the  detector  alias  frequency  0.5a'1,  it  is  necessary  to  choose  a  higher  value  to  achieve  super¬ 
resolution.  The  specified  value  of  c,  corresponds  to  the  second  zero  of  the  MTF  of  the  sampling 
process  for  frequency  measurements  along  the  f\  direction  (fi  =  0). 15-17 

MTF(/p/2)  =  sinc(a;c/1  )sinc(a>,/2 )  (87) 
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Figure  3  shows  a  plot  of  the  reconstruction  filters  versus  frequency.  The  two  filters  almost  perfectly 
430  match  each  other  at  low  frequencies  but  diverge  at  high  frequencies,  since  the  SA  filter  is  intended 
to  suppress  high  frequency  noise. 

III.  A.  1.  Individual  Projections 

At  a  fixed  distance  (112)  of  30.0  mm  from  the  chest  wall,  Figure  4(a)-(b)  shows  a  cross 
435  section  of  signal  versus  detector  position  u\  for  the  central  projection  (n  =  0)  and  an  oblique 

projection  (n  =  7).  The  112  displacement  lies  between  the  chest  wall  and  nipple  of  a  typical  breast. 

In  the  recent  development  of  a  physical  3D  anthropomorphic  phantom  for  image  quality  assessment 
in  DM  and  DBT,18’ 19  Carton  et  al.  modeled  a  distance  of  65.0  mm  between  the  chest  wall  and 
nipple  for  an  average  breast  size  of  450  mL.  The  uo  displacement  considered  in  Figure  4(a)-(b)  thus 
440  corresponds  to  a  position  approximately  halfway  between  the  chest  wall  and  nipple  of  this  phantom. 

In  Figure  4(c)-(d),  the  modulus  of  the  Fourier  transform  of  detector  signal  is  plotted  versus 
frequency  f\,  assuming/?  =  0  [Eq.  (64)].  The  central  and  oblique  projections  are  similar  in  that  they 
both  represent  a  high  frequency  input  as  if  it  were  a  lower  frequency.  The  Fourier  transform  of 
either  projection  does  not  occur  at  the  input  frequency  5.00  lp/mm  but  instead  occurs  at  a  lower 
445  frequency  as  evidence  of  aliasing. 

The  two  projections  and  their  frequency  spectra  are  also  plotted  in  Figure  4  for  an  infinite 
source-to-COR  distance  (/?)  with  no  other  changes  in  the  acquisition  parameters.  This  limiting  case 
transforms  the  divergent  beam  geometry  into  a  parallel  beam  geometry.  Consequently,  the  x-ray 
angle  relative  to  the  normal  to  the  detector  does  not  vary  with  position  (u\,  uj)  by  Eq.  (15)  but 
450  instead  is  always  yjn  -  yn  for  the  nth  projection. 

In  the  parallel  beam  geometry,  the  central  projection  represents  the  input  frequency  as  if  it 
were  a1  -  f  h  or  0.3a"1.  As  a  result,  the  Fourier  transform  has  a  major  peak  at  0.3a"1  (2.14  lp/mm), 
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and  has  minor  peaks  at  0.7a'1  (5.0  lp/mm),  1.3a'1  (9.29  lp/mm),  and  1.7a'1  (12.14  lp/mm).  Unlike 
the  parallel  beam  geometry,  the  divergent  beam  geometry  magnifies  the  input  so  that  it  projects 
455  onto  the  x-ray  converter  with  the  frequency  fo/M 

(88) 

h-z0 

where  M  denotes  the  magnification.  With  a  source-to-COR  distance  (h)  of  70.0  cm  and  an  object- 
to-detector  distance  (zo)  of  50.0  mm,  M  is  1 .077.  As  a  result  of  the  magnification,  the  peaks  in  the 
Fourier  transform  of  detector  signal  occur  at  different  frequencies  than  the  parallel  beam  geometry. 
460  Accordingly,  these  Fourier  peaks  occur  at  a'1  -  fo/M  (2.50  lp/mm),  fo/M  (4.64  lp/mm),  2a'1  -  fo/M 
(9.64  lp/mm),  and  a'1  +  fo/M  (1 1.78  lp/mm).  The  Fourier  transform  of  the  most  oblique  projection 
peaks  at  similar  frequencies  as  the  central  projection. 

III.  A.  2.  SBP  Reconstruction 

465  Figure  5(a)  shows  SBP  reconstruction  versus  position  (x)  measured  parallel  to  the  chest  wall 

side  of  the  breast  support,  performed  at  the  distance  y  =  30.0  mm  from  the  chest  wall  and  at  the 
height  z  =  zo  =  50.0  mm  above  the  breast  support.  Unlike  an  individual  projection,  SBP 
reconstruction  can  resolve  the  input  frequency  5.00  lp/mm.  This  property  arises  because  the 
oblique  projections  give  information  about  the  input  which  is  not  present  in  the  central  projection 
470  alone  [Figure  4(b)],  Although  not  explicitly  plotted  in  Figure  5(a),  it  can  be  shown  that  super¬ 
resolution  is  present  across  a  broad  range  of  x  and  y  positions  in  the  reconstructed  volume. 

The  SBP  Fourier  transform  [Eq.  (86)]  correctly  possesses  its  major  peak  at  5.00  lp/mm.  The 
major  peak  of  an  individual  projection,  occurring  at  2.50  lp/mm,  is  now  highly  suppressed  in 
magnitude  [Figure  5(c)], 

475 
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III.  A.  3.  FBP  Reconstruction 

FBP  reconstructions  are  now  performed  with  either  the  RA  filter  alone  or  the  RA  and  SA 
filters  together.  In  the  spatial  domain,  these  reconstructions  are  plotted  versus  position  (x)  parallel 
to  the  chest  wall  side  of  the  breast  support,  assuming  y  =  30.0  mm  and  z  =  50.0  mm  [Figure  5(b)], 
480  Figure  5(b)  demonstrates  that  reconstruction  filters  smoothen  pixilation  artifacts  found  in  the  SBP 
reconstruction.  In  addition,  Figure  5(b)  shows  that  reconstructing  with  the  RA  filter  alone  yields 
greater  modulation  than  reconstructing  with  the  RA  and  S A  filters  together.  The  modulation  for 
reconstruction  with  the  RA  filter  alone  is  41 .0%,  yet  the  modulation  for  reconstruction  with  the  RA 
and  SA  filters  together  is  29.8%.  It  is  expected  that  reconstruction  with  the  RA  filter  alone  has 
485  greater  modulation,  since  the  amplitude  of  this  filter  exceeds  that  of  the  RA  and  SA  filters  together 
at  the  input  frequency  (Figure  3).  Importantly,  the  modulation  of  either  FBP  reconstruction 
technique  is  well  above  the  limit  of  resolution  for  typical  imaging  systems,  which  is  often  taken  to 
be  5%.  In  addition,  the  modulation  of  either  FBP  reconstruction  technique  is  greater  than  that  of 
SBP  reconstruction  (18.4%). 

490  Although  reconstruction  with  the  RA  filter  alone  has  the  benefit  of  greater  modulation  than 

reconstruction  with  the  RA  and  SA  filters  together,  the  trade-off  is  greater  spectral  leakage  at  very 
high  frequencies.  In  fact,  the  amplitude  of  the  high  frequency  spectral  leakage  is  greater  with  the 
RA  filter  alone  than  with  SBP.  In  experimental  practice,  reconstruction  with  the  RA  filter  alone 
also  increases  the  presence  of  noise,  which  tends  to  occur  at  high  frequencies. 

495 

III.  B.  Input  Frequency  Directed  Perpendicular  to  the  Chest  Wall 

It  is  now  demonstrated  that  the  existence  of  super-resolution  is  dependent  on  the 
directionality  of  the  input  frequency.  Super-resolution  arises  because  of  sub-pixel  detector  element 
shifts  in  the  image  of  an  object  with  each  increasing  projection  angle.  In  order  to  investigate  the 
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Assuming  that  z  =  50.0  mm,  Figure  6  shows  the  magnitude  of  this  translational  shift  versus  position 
y  within  two  planes,  x  =  0  and  x  =  -30.0  mm,  comparing  the  central  projection  and  an  oblique 
510  projection  ( n\  =  0,  m  =  7)  as  well  as  two  oblique  projections  (n\  =  -7,  tn  =  7).  Throughout  the 

remainder  of  this  work,  a  plane  defined  by  a  fixed  value  of  x  will  be  termed  an  AP/SS  plane  since  it 
has  extent  in  both  the  anteroposterior  (AP)  and  source-to-support  (SS)  directions.  Although  the  SS 
direction  technically  varies  with  position  on  the  breast  support  due  to  the  divergence  of  the  x-ray 
beam,  it  is  assumed  to  be  equivalent  to  the  z  direction  for  the  purpose  of  this  work.  As  such,  the  SS 
515  direction  lies  along  the  same  axis  as  the  source-to-image  distance  (SID),  or  the  length  between  the 
focal  spot  and  the  origin  O  for  the  central  projection  (Figure  1).  In  a  cranial-caudal  (CC)  view,  an 
AP/SS  plane  is  thus  a  sagittal  plane  through  the  breast.  By  contrast,  in  a  mediolateral  oblique 
(MLO)  view,  the  same  plane  is  at  an  approximately  45°  angle  relative  to  the  sagittal  and  transverse 
planes  through  the  breast. 

520  In  the  mid  AP/SS  plane  (x  =  0),  translational  shifts  between  projections  are  minimal  in  a 

typical  sized  breast.  For  example,  with  x  =  0,  y  =  30.0  mm,  and  z  =  50.0  mm,  the  translational  shift 
between  the  central  projection  and  an  oblique  projection  is  0.009  mm  (6.52%  of  detector  element 
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length),  and  the  translational  shift  between  the  two  most  oblique  projections  is  zero.  For  this 
reason,  super-resolution  along  the  y  direction  is  simply  not  achievable  within  the  mid  AP/SS  plane. 
525  As  illustrated  in  Figure  7(a),  SBP  reconstruction  at  x  =  0  in  the  region  ye  [29.4  mm,  30.6  mm] 
resembles  a  single  projection. 

In  Figure  6,  it  is  demonstrated  that  the  112  translational  shift  between  projections  increases  as 
the  magnitude  of  the  distance  x  increases.  For  example,  with  x  =  -30.0  mm  and  y  =  30.0  mm,  the 
translational  shift  between  the  central  projection  and  the  oblique  projection  is  0.047  mm  (33.4%  of 
530  detector  element  length),  and  the  translational  shift  between  the  two  most  oblique  projections  is 
0.1 12  mm  (80.0%  of  detector  element  length).  Because  these  translational  shifts  are  sufficiently 
large,  SBP  reconstruction  [Figure  7(c)]  shows  super-resolution  at  x  =  -30.0  mm  over  the  region 
ye  [29.4  mm,  30.6  mm].  Unlike  SBP  reconstruction  for  an  input  frequency  oriented  along  the  x 
direction  [Figure  5(a)],  the  amplitudes  of  the  peaks  in  Figure  7(c)  are  noticeably  different  from  each 
535  another,  indicating  the  presence  of  reconstruction  artifacts.  It  can  be  shown  that  these  artifacts  are 
minimized  by  increasing  the  distance  y  from  the  chest  wall,  since  the  U2  translational  shifts  between 
projections  increase  withy  (Figure  6). 

The  SBP  reconstructions  in  Figure  7  can  be  analyzed  further  by  computing  their  ID  Fourier 
transform  along  the  y  direction.  To  show  differences  in  these  Fourier  transforms  at  fixed  values  of 
540  x,  we  choose  not  to  transform  over  both  x  and  y  as  given  by  Eq.  (85). 
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For  additional  proof  that  super-resolution  is  not  achievable  within  the  mid  AP/SS  plane,  Figure  7(b) 
shows  that  the  major  Fourier  peak  at  x  =  0  occurs  well  below  the  input  frequency.  By  contrast,  the 
545  major  Fourier  peak  at  x  =  -30.0  mm  [Figure  7(d)]  matches  the  input  frequency,  5.00  lp/mm. 
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Although  not  shown  in  the  plot,  it  can  be  demonstrated  that  spectral  leakage  is  reduced  by 
increasing  the  magnitude  of  the  distance  x  relative  to  the  mid  AP/SS  plane. 

As  a  final  point  in  this  section,  it  is  important  to  note  that  by  applying  the  filters  in  Figure  3 
to  the  SBP  reconstructions  of  Figure  7,  the  modulation  effectively  vanishes  (graph  not  shown).  This 
550  finding  arises  because  filtering  is  applied  only  within  the  plane  of  the  chest  wall  [Eq.  (40)].  An 

input  frequency  oriented  along  they  direction  contributes  a  component  of  0  lp/mm  within  the  plane 
of  the  chest  wall;  since  the  reconstruction  filters  vanish  at  0  lp/mm  (Figure  3),  FBP  reconstructions 
are  expected  to  have  no  modulation.  For  this  reason,  future  research  on  filter  optimization  is 
merited  as  described  in  the  Discussion  section. 

555 

IV.  EXPERIMENTAL  AND  CLINICAL  VALIDATION 

Using  a  high  contrast  bar  pattern  phantom,  we  have  experimentally  verified  the  existence  of 
super-resolution  in  DBT.  The  phantom  was  taped  beneath  the  compression  paddle  (24  cm  x  29  cm) 
of  the  Selenia  Dimensions  system,  and  placed  2.5  cm  above  the  breast  support.  With  the  alternating 
560  light  and  dark  bands  of  the  phantom  spanning  a  6.0  mm  length,  the  line  pairs  ranged  in  frequency 
from  1.0  lp/mm  to  10.0  lp/mm.  To  match  the  simulation  of  Section  IIIA,  the  frequency  5.0  lp/mm 
was  oriented  along  the  x  direction  parallel  to  the  chest  wall  side  of  the  breast  support.  Also, 
following  the  simulation,  the  bar  patterns  at  5.0  lp/mm  covered  the  region  xe  [-0.6  mm,  0.6  mm], 
and  the  edge  of  the  bar  patterns  near  the  numeral  “5”  (Figure  8)  was  positioned  slightly  greater  than 
565  y  =  30  mm  from  the  chest  wall.  Using  the  large  (0.3  mm  nominal)  focal  spot  and  a  CC  view,  15 
projections  were  acquired  at  30  kVp  and  14  mAs  with  a  W/Al  target-filter  combination.  The 
technique  for  determining  the  optimal  mAs  with  photo-timing  is  described  in  our  previous  work.4 

Reconstruction  was  subsequently  performed  using  a  backprojection  filtering  (BPF) 
commercial  prototype  reconstruction  solution  (Briona™,  Real  Time  Tomography,  Villanova,  PA). 
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570  Although  it  is  possible  to  reconstruct  on  a  non-pixilated  grid  using  analytical  modeling,  a  pixilated 
grid  was  required  for  the  experimental  data.  In  order  to  ensure  that  high  frequencies  can  be 
resolved  in  the  plane  of  the  reconstruction,  the  pixel  size  of  the  reconstruction  grid  (20.44  pm)  was 
chosen  to  be  significantly  smaller  than  that  of  the  detector  elements  (140  pm).  Consequently,  the 
alias  frequency  of  the  reconstruction  grid  (24.46  lp/mm)  was  substantially  higher  than  the  alias 
575  frequency  of  the  detector  (3.57  lp/ mm) . 

Figure  8  shows  that  the  central  projection  correctly  resolves  frequencies  below  the  detector 
alias  frequency,  3.57  lp/mm.  At  the  next  highest  frequency  (4.0  lp/mm),  one  would  expect  to  see  24 
line  pairs  spanning  a  6.0  mm  length.  Instead,  less  than  24  line  pairs  are  visible,  and  Moire 
patterns  are  present.  Finally,  at  5.0  lp/mm,  only  16  line  pairs  are  evident  within  a  6.0  mm  length, 
580  indicating  that  the  pattern  is  incorrectly  represented  as  a  frequency  between  2.0  and  3.0  lp/mm. 

Unlike  the  central  projection,  BPF  reconstruction  can  resolve  frequencies  higher  than  the 
detector  alias  frequency  (Figure  9).  In  fact,  up  to  6.0  lp/mm  (36  line  pairs  spanning  6.0  mm)  can  be 
observed  at  the  correct  orientation  with  no  Moire  patterns.  At  7.0  lp/mm,  the  signal  becomes  too 
faint  to  distinguish  bar  patterns.  This  finding  arises  because  the  MTF  of  the  reconstruction  is 
585  reduced  with  increasing  frequency;  recall  from  Eq.  (87)  that  the  MTF  of  the  sampling  process 

vanishes  at  the  frequency  f\  =  a 1  (7.14  lp/mm),  assuming =  0.  As  expected  from  the  analytical 
modeling,  it  should  be  noted  that  super-resolution  along  the  x  direction  was  observed  over  many 
different  x-ray  acquisitions  in  which  the  bar  pattern  phantom  was  placed  at  various  positions  in  the 
imaging  volume. 

590  By  rotating  the  bar  pattern  phantom  90°,  the  potential  for  super-resolution  orthogonal  to  the 

chest  wall  was  also  analyzed  (Figure  10).  To  orient  the  reader  with  the  positioning  of  the  phantom, 
it  is  important  to  note  that  the  left  edges  of  the  even  numerals  “4”  and  “6”  were  aligned  on  the  mid 
AP/SS  plane  (x  =  0),  and  that  the  separation  between  4.0  and  5.0  lp/mm  was  positioned  at  a 
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displacement  y  =  30  mm  from  the  chest  wall.  As  expected  from  the  analytical  modeling  (Section 
595  IIIB),  the  extreme  left  regions  of  the  bar  patterns  show  aliasing  of  high  frequencies  due  to  their 

proximity  to  the  mid  AP/SS  plane.  Super-resolution  is  only  present  at  the  extreme  right  of  the  bar 
patterns  (5.0  and  6.0  lp/mm),  where  the  magnitude  of  the  distance  x  relative  to  the  mid  AP/SS  plane 
is  approximately  30  mm  or  greater. 

It  is  important  to  note  that  the  extreme  right  of  the  bar  patterns  at  4.0  lp/mm  does  not  display 
600  super-resolution  as  cleanly  as  the  extreme  right  of  the  bar  patterns  at  5.0  and  6.0  lp/mm  due  to  the 
presence  of  Moire  patterns.  Recall  that  the  phantom  is  positioned  so  that  lower  frequencies  are 
closer  to  the  chest  wall.  Because  the  U2  translational  shift  between  projections  is  minimized  with 
decreasing  distance  from  the  chest  wall  (Figure  6),  it  is  expected  that  super-resolution  along  they 
direction  should  not  be  achievable  at  positions  too  close  to  the  chest  wall. 

605  In  breast  imaging,  super-resolution  has  applications  in  improving  the  visualization  of  small 

microcalcification  clusters,  which  are  an  early  indication  of  cancer.  This  concept  is  illustrated  in 
Figure  1 1 .  In  the  figure,  the  left  image  is  generated  by  bilinearly  interpolating  a  BPF  reconstruction 
performed  using  pixels  of  length  140  pm.  The  bilinear  interpolation  yields  a  final  result  (a)  with 
pixels  of  length  35  pm  (hence,  four- fold  magnification).  By  contrast,  the  right  image  (b)  is 
610  generated  by  performing  a  BPF  reconstruction  using  pixels  of  length  35  pm.  Super-resolution  is 
present  in  (b),  not  (a).  A  few  small  microcalcifications  that  cannot  be  discerned  in  image  (a)  are 
indeed  evident  in  image  (b),  and  finer  detail  in  their  overall  structural  morphology  can  be  observed 
in  image  (b).  Having  a  clearer  image  of  the  structural  features  of  microcalcifications  is  useful  as  it 
yields  greater  diagnostic  information  about  their  pathology,22  and  can  help  to  diagnose  breast  cancer 
615  in  its  early  stages  of  development.  For  this  reason,  super-resolution  in  DBT  has  enormous  potential 
for  future  applications  in  breast  imaging. 
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V.  DISCUSSION 

In  DBT  reconstructions  using  grids  with  the  same  pixel  size  as  the  detector  elements,  the 
620  highest  frequency  that  can  be  resolved  in  each  reconstructed  slice  is  the  detector  alias  frequency. 
This  study  demonstrates  that  reconstruction  grids  with  much  smaller  pixilation  display  super¬ 
resolution,  or  visibility  of  higher  frequencies.  Super-resolution  arises  because  the  image  of  the 
object  is  shifted  in  sub-pixel  detector  element  increments  with  each  increasing  projection  angle. 

Super-resolution  was  first  demonstrated  analytically  by  calculating  the  reconstruction  of  a 
625  sinusoidal  input  whose  frequency  was  oriented  along  the  x  direction  parallel  to  the  chest  wall  side  of 
the  breast  support.  Using  an  infinitesimally  fine  reconstruction  grid,  it  was  shown  that  both  SBP 
and  FBP  can  resolve  higher  frequencies  than  a  single  projection  over  a  broad  range  of  positions  in 
the  imaging  volume.  FBP  reconstructions  were  performed  either  with  the  RA  filter  alone  or  with 
the  RA  and  SA  filters  together.  Although  reconstruction  with  the  RA  filter  alone  has  the  benefit  of 
630  greater  modulation  in  the  spatial  domain,  it  presents  the  trade-off  of  increased  noise  and  spectral 

leakage  at  high  frequencies.  In  rotating  the  sine  plate  by  90°,  the  visibility  of  high  frequencies  was 
found  to  be  dependent  upon  position  within  the  reconstructed  volume.  Positions  with  super¬ 
resolution  along  the  y  direction  must  be  sufficiently  displaced  relative  to  the  chest  wall  (y  =  0)  and 
to  the  mid  AP/SS  plane  (x  =  0). 

635  We  have  experimentally  observed  super-resolution  in  images  of  bar  patterns  using  the 

Selenia  Dimensions  system  and  a  commercial  BPF  reconstruction  algorithm.  A  single  projection 
showed  classical  signs  of  aliasing,  including  Moire  patterns  and  the  visibility  of  fewer  line  pairs 
than  expected.  By  contrast,  reconstructions  using  very  fine  grids  resolved  frequencies  higher  than 
the  alias  frequency  of  the  detector.  For  the  two  orientations  of  the  bar  patterns,  the  presence  of 
640  super-resolution  was  verified  at  positions  predicted  from  analytical  modeling.  In  select  clinical 
examples,  super-resolution  was  found  to  improve  the  visibility  of  microcalcifications,  which  are 
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early  indicators  of  breast  cancer  in  many  women.  As  a  result  of  super-resolution,  the  morphology 
of  microcalcifications  was  seen  more  clearly,  and  small  microcalcifications  that  were  not  otherwise 
visible  became  apparent. 

645  Super-resolution  is  a  particularly  useful  property  for  x-ray  systems  that  employ  binning 

when  switching  from  2D  to  3D  imaging  modes.  For  example,  in  the  Selenia  Dimensions  system, 
the  DM  detector  element  dimensions  are  70  pm  x  70  pm,  whereas  the  DBT  detector  element 
dimensions  are  140  pm  x  140  pm.  Binning  has  the  benefit  of  lowering  the  readout  time,  but 
presents  the  drawback  of  reducing  the  alias  frequency  of  the  detector.  Initially,  it  would  seem  that 
650  binning  should  make  DBT  less  capable  of  resolving  high  frequency  information,  such  as 

microcalcifications.  However,  the  existence  of  super-resolution  in  the  reconstruction  may  counter 
the  trade-offs  of  binning. 

Some  of  the  limitations  of  this  study  and  directions  for  future  investigation  are  now  noted. 

In  calculating  detector  signal,  this  paper  assumes  that  the  MTF  of  a-Se  in  drift  mode  is  unity.  While 
655  this  assumption  is  valid  for  normal  x-ray  incidence,  it  is  less  justifiable  with  oblique  x-ray 
incidence.23"27  Que  and  Rowlands  proposed  the  first  analytical  model  of  the  optical  transfer 
function  (OTF)  of  a-Se  in  drift  mode  for  all  incident  angles.23  Their  work  was  later  validated  by 
Hajdok  and  Cunningham  with  Monte  Carlo  simulations.24  Denoting  //se  as  the  attenuation 
coefficient  of  Se  and  L  as  the  thickness  of  the  photoconductor,  the  OTF  at  each  frequency /is 


The  MTF  is  the  normalized  modulus  of  the  OTF.  For  more  thorough  modeling,  signal  in  the  x-ray 
converter  should  be  convolved  with  the  point  spread  function  (PSF)  of  a-Se  before  detector  element 
sampling  is  performed,  where  the  PSF  is  determined  from  the  OTF  using  Fourier  theory.  It  is 
important  to  model  MTF  degradation  for  measurements  near  the  edge  of  the  detector  opposite  the 
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665  chest  wall,  as  the  incident  angle  deviates  considerably  from  the  normal.  Upon  examining  6n  across 
multiple  projections  in  the  Selenia  Dimensions  detector,  it  can  be  shown  that  the  maximum  incident 
angle  is  approximately  25°.  Assuming  200  pm  thick  a-Se  and  20  keV  x  rays28'30  for  which  pse  is 
20.5  mm'1,31  the  corresponding  MTF  at  5.0  lp/mm  is  85.8%. 

While  it  is  important  to  consider  MTF  degradation  at  positions  distal  to  the  chest  wall,  it  is 
670  less  critical  for  positions  close  to  the  chest  wall.  For  example,  in  the  central  projection  at  the 
position  ui  =  U2  =  30.0  mm,  the  incident  angle  is  3.47°,  and  the  MTF  at  5.0  lp/mm  is  99.7%. 
Consequently,  for  the  purpose  of  this  work,  an  x-ray  converter  with  MTF  of  unity  was  assumed. 

In  addition  to  modeling  the  MTF  of  the  x-ray  converter,  the  analytical  model  of  super¬ 
resolution  can  be  refined  by  modeling  the  MTF  of  the  focal  spot.  Although  this  paper  assumes  a 
675  point- like  focal  spot,  future  work  should  consider  the  MTF  degradation  with  increasing  focal  spot 
size  and  increasing  focal  spot  motion  during  a  continuous  scan  of  the  projections.  ’  Detector 
lag  and  ghosting34'36  should  also  be  simulated,  and  the  presence  of  shot  noise37  should  be  modeled 
at  various  dose  levels.  Because  this  work  considers  a  high  contrast  input  frequency  either  with  the 
analytical  simulation  or  with  the  experimental  bar  patterns,  it  was  not  necessary  to  model  the 
680  presence  of  noise  at  different  dose  levels.  Future  studies  on  super-resolution  with  low  contrast  input 
frequencies  will  require  a  noise  simulation,  as  the  visibility  of  the  patterns  should  be  influenced  by 
dose.  Finally,  because  the  linear  attenuation  coefficient  of  an  input  is  energy  dependent, 
polyenergetic  x-ray  spectra38'40  should  also  be  simulated  in  the  analytical  model.  This  work 
implicitly  assumes  a  monoenergetic  x-ray  beam. 

685  In  CT,  the  conventional  low  frequency  filter  is  the  RA  filter13, 14  which  increases  linearly 

with  frequency  from  zero  (Figure  3).  Since  filtering  was  only  applied  in  the  plane  of  the  x-ray  tube 
motion,  this  work  has  demonstrated  that  the  RA  filter  is  not  suited  for  imaging  frequencies 
perpendicular  to  the  chest  wall,  since  the  modulation  of  the  reconstruction  vanishes  (Section  IIIB). 
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Lauritsch  and  Harer  have  shown  that  replacing  the  conventional  filter  with  a  polynomial  filter 
690  improves  the  quality  of  the  DBT  reconstruction.  Their  filter  rises  sharply  at  low  frequencies, 
plateaus  to  unity  at  intermediate  frequencies,  and  falls  in  value  at  high  frequencies.41,42  Future 
work  should  consider  a  polynomial  filter  for  analytical  modeling  of  super-resolution,  since  the  non¬ 
zero  offset  will  never  yield  zero  modulation  for  any  orientation  of  the  input  frequency.  It  should  be 
noted  that  the  filters  used  in  the  experimental  reconstructions  of  bar  patterns  oriented  perpendicular 
695  to  the  chest  wall  (Figure  10)  have  a  non-zero  offset  at  0  lp/mm,  unlike  the  RA  filter  used  in  the 
analytical  modeling. 

Because  super-resolution  has  important  clinical  applications  in  improving  the  visibility  of 
microcalcifications,  future  work  should  ultimately  transition  from  modeling  a  sinusoidal  input  to 
simulating  microcalcifications  in  a  breast  background.43  Using  model  observers,  improvements  in 
700  the  visibility  of  microcalcifications  should  be  assessed  with  image  reconstructions  at  varying  grid 
sizes.  It  would  be  useful  to  determine  that  coarsest  grid  size  at  which  the  benefits  of  super¬ 
resolution  are  achieved  among  observers,  as  reconstructions  on  coarser  grids  require  less  memory 
for  data  storage. 

705  VI.  CONCLUSION 

To  our  knowledge,  this  work  is  the  first  to  demonstrate  the  existence  of  super-resolution  in 
DBT.  An  analytical  model  of  super-resolution  was  developed  by  calculating  the  reconstruction  of  a 
high  frequency  sinusoidal  input.  While  a  single  projection  cannot  resolve  frequencies  higher  than 
the  alias  frequency  of  the  detector,  a  reconstruction  on  a  very  fine  grid  can  resolve  these 
710  frequencies.  Super-resolution  is  made  possible  by  the  sub-pixel  detector  element  shifts  in  the  image 
of  the  object  between  projections. 
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Using  a  bar  pattern  phantom,  we  have  experimentally  verified  the  existence  of  super¬ 
resolution  in  DBT.  In  considering  an  input  frequency  oriented  both  parallel  and  perpendicular  to 
the  chest  wall  side  of  the  breast  support,  the  experimental  images  confirmed  the  presence  of  super- 
715  resolution  at  positions  predicted  by  analytical  modeling.  Ultimately,  super-resolution  was  shown  to 
improve  the  visibility  of  fine  structural  details  of  microcalcifications  in  clinical  DBT  images. 
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APPENDIX  A:  DETECTOR  SIGNAL  FOR  AN  INPUT  FREQUENCY  DIRECTED 
PERPENDICULAR  TO  THE  CHEST  WALL 


740 


This  appendix  calculates  detector  signal  for  an  input  frequency  perpendicular  to  the  chest 
wall.  Under  this  assumption,  the  input  rectangular  prism  of  thickness  s  has  a  linear  attenuation 
coefficient  ju(x,  y,  z )  which  varies  sinusoidally  along  the  y  direction  with  frequency  fo 


ju(x,  y ,  z)  =  C  •  cos  [2nf0  (y  -  y0 )]  •  rect 


[(Al):  cf.  (1)] 


where  yo  is  a  translational  shift  in  the  waveform  relative  to  the  origin.  The  amplitude  C  of  the 
waveform  is  equivalent  to  1/e  upon  normalizing  total  attenuation  along  the  z  direction.  The  ID 
Fourier  transform  of  Eq.  (Al)  along  the  y  direction  peaks  at  the  frequencies  fy  =  ±/o  and  vanishes  at 
745  all  other  frequencies,  following  a  formula  similar  to  Eq.  (4)  with  the  exchange  of  xo  for  yo  and  fx  for 
fy.  Using  Eqs.  (17)-(19)  and  Eq.  (22),  total  x-ray  attenuation  versus  position  (u\,  U2 )  along  the  plane 
of  the  rotated  detector  is  calculated  for  the  nth  projection  as 


Aju(n )  =  Kn  f  "  cos(2;r  f0u2w  +  A)dw 

J  w* 

Kn  ^sin  [2  +  A J  -  sin  |^2 n f0u2w +  +  Aj  j 

27ifu2 


750  where 


[(A2):  cf.  (23)] 

[(A3):  cf.  (24)] 


A  =  -2 nfy,.  [(A4):  cf.  (26)] 

Following  the  sum-to-product  trigonometric  identity  given  in  Eq.  (27),  one  may  rewrite  Eq.  (A3)  as 
A  An)  =  Kn{wn  -  w„+)cos[^/0w2(wn+  +  w~)  +  A]sinc[/0w2(wn_  -  w„+)]  [(A5):  cf.  (28)] 


£K„  COS 


2^/0m2(/  +  /?cos^-Zo) 
l  +  hcosy/n  -ux  sin yn 


+  A 


sine 


sfu2 


l  +  h  cos  y/ n  -ux  siny„ 


l  +  h  cos  y/n  -  ux  sin  yn 


~  •  [(A6):  cf.  (29)] 
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755 


The  logarithmically-transformed  signal  in  the  mth  detector  element  for  the  nh  projection  is  now 
determined  from  Eq.  (31).  The  midpoint  formula44  for  approximating  this  double  integral  is 


V/u(m,n 


)=  lim  —  V 


y  h=X 


lmi— X^/d/A,./v,«) 


Jx^  J 


jx=  1 


[(A7):  cf  (35)] 


where 


Afu{jxJy,n)  =  Afu{n\uu2)_ 


j  -1/2  1 

• v - — 

Jr  2 


7y_l/2 


[(A8):  c/  (36)] 


760  FBP  reconstruction  now  follows  from  Eq.  (42). 


APPENDIX  B:  NOMENCLATURE 


Symbol 

* 

X 

E 

Aju(n) 

Aju(n) 

8 

P//(m,«) 

T 

4 

R3 

S/u(ux,u2) 


Meaning 

Dot  product  operator. 

Convolution  operator. 

Cross  product  operator. 

Set  membership. 

Total  attenuation  for  the  nth  projection. 

A  useful  approximation  for  total  attenuation  [Eqs.  (32)-(33)]. 

Backprojection  operator. 

Signal  in  the  mth  detector  element  for  the  nh  projection. 
Fourier  transform  operator  (subscript  denotes  dimension). 
Path  length  through  the  input  for  the  »th  projection. 

Euclidean  3 -space. 

Raw  signal  at  coordinate  (u\,  112 )  on  the  rotated  detector. 
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X 

z 

z* 

Yn 

r 

m  n 

S 

A  y/ 

AM2(«j,n2) 


On 


K 

A 

M 


MSe 


X-ray  transform  operator. 

Set  of  integers. 

Set  of  non-negative  integers. 

Angle  of  rotation  of  the  detector  relative  to  the  x  axis  for  the  nth  projection. 

Angle  of  backprojection  within  the  plane  of  the  detector  [Eqs.  (53)-(54)]. 

Delta  function. 

Angular  spacing  between  projections. 

Translational  shift  in  w2  coordinate  of  incident  ray  comparing  projection  numbers  n\ 
and  «2- 

Thickness  of  sine  plate  (Figure  1). 

Angle  of  x-ray  incidence  relative  to  the  normal  to  the  detector  ( 6mn  denotes  the 

special  case  at  the  centroid  of  the  mth  detector  element  for  the  »,h  projection). 

A  quantity  defined  by  Eq.  (25). 

A  quantity  defined  by  Eq.  (26). 

A  quantity  defined  by  Eq.  (A6). 

X-ray  linear  attenuation  coefficient  of  input. 

X-ray  linear  attenuation  coefficient  of  a-Se  photoconductor. 


4 


P\  9  Pi 


(J 

jmn 


<!> 


Truncation  frequency  of  reconstruction  filter. 

Quantities  defined  by  Eqs.  (66)-(67). 

Terms  defined  by  Eqs.  (68)-(72)  used  to  simplify  intermediate  calculations,  where  j 
varies  from  1  to  5. 

Reconstruction  filter. 
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Vn 

AP 


AP/SS 


ax,av 


bx,b2 


C 

cc 

COR 

CT 

dn 

DBT 

DM 

/ 

FBP 

g 

h 

i 

d xmn 

dymni*) 


Nominal  projection  angle. 

Anteroposterior  (in  breast  x-ray  imaging,  the  direction  perpendicular  to  the  chest 
wall). 

Descriptive  acronym  for  a  plane  with  extent  along  the  anteroposterior  (AP)  and 
source-to-support  (SS)  directions. 

Detector  element  dimensions  in  the  x  andy  directions;  if  the  x  and  y  subscripts  are 
removed,  the  detector  element  is  square  ( ax  =  ay  =  a). 

Real  numbers  used  to  illustrate  a  sum-to-product  trigonometric  identity  [Eq.  (27)]. 

Amplitude  of  sine  input  taken  to  be  1/e. 

Cranial-caudal. 

Center-of-rotation  of  x-ray  tube  motion. 

Computed  tomography. 

Distance  between  points  G  and  O  (Figure  2). 

Digital  breast  tomosynthesis. 

Digital  mammography. 

Spatial  frequency  (fo  denotes  the  input  frequency). 

Filtered  backprojection. 

Gear  ratio  of  detector. 

Source-to-COR  distance  for  rotating  x-ray  tube. 

Imaginary  unit  given  as  V-T . 

An  integral  defined  by  Eq.  (81). 

An  integral  defined  by  Eq.  (76). 
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/  Distance  between  the  COR  and  the  midpoint  of  the  chest  wall  side  of  the  detector. 

Ip  Line  pairs. 

m  A  doublet  with  coordinates  ( mx ,  my)  used  for  labeling  detector  elements. 

M  Magnification. 

MLO  Mediolateral  oblique. 

MRI  Magnetic  resonance  imaging. 

MTF  Modulation  transfer  function. 

n  Projection  number. 

N  Total  number  of  projections. 

OTF  Optical  transfer  function. 

RA  Ramp  filter. 

SA  Spectrum  apodization  filter. 

SBP  Simple  backprojection. 

SID  Source-to-image  distance  (commonly  measured  between  the  focal  spot  and  the 

midpoint  of  the  chest  wall  side  of  the  detector  in  the  central  projection). 

SNR  Signal-to-noise  ratio. 

SS  Source-to-support  (defined  to  be  synonymous  with  the  z  direction). 

tx ,  t2  Affine  parameters  of  the  x-ray  transform. 

TFT  Thin-film  transistor. 

Wj ,  u2  Position  in  the  plane  of  the  rotated  detector  (parallel  and  perpendicular  to  the  chest 

wall,  respectively). 

w  Parameter  ranging  between  0  and  1  in  the  equation  of  the  x-ray  beam  between  the 

focal  spot  and  the  incident  point  on  the  detector  [Eq.  (16)]. 
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Value  of  w  at  the  entrance  ( w+n  )  and  exit  (w~)  points  of  the  x-ray  beam  through  the 
sine  plate  (Figure  1)  for  the  nth  projection. 

Position  parallel  to  the  chest  wall  side  of  the  breast  support;  rotation  by  the  angle 
yn  about  the  y  axis  yields  x'n . 

Translational  shift  in  the  input  waveform  along  the  x  direction. 

Position  perpendicular  to  the  chest  wall;  it  is  equivalent  to  y'n . 

Translational  shift  in  the  input  waveform  along  the  y  direction. 

Position  perpendicular  to  the  plane  of  the  breast  support;  rotation  by  the  angle  yn 
about  the  y  axis  yields  z' . 

Central  height  of  the  input  relative  to  the  midpoint  of  the  chest  wall  side  of  the 


detector. 
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X-Ray  Beam 


Cross  Section 
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l l ■■■ ■■ ll 

\  \  k: 


Figure  1.  The  3D  input  object  is  a  rectangular  prism  whose  linear  attenuation  coefficient  varies 
sinusoidally  with  position  x  parallel  to  the  chest  wall  side  of  the  breast  support.  A  2D  cross  section 
of  the  input  object  through  the  plane  of  the  chest  wall  is  shown  (figure  not  to  scale).  In  acquiring 
projection  images,  the  x-ray  tube  rotates  within  the  xz  plane  about  point  B,  and  the  detector 


780  simultaneously  rotates  about  they  axis.  The  primed  unit  vectors  ij  and  jn'  define  the  coordinate 


axes  of  the  plane  of  the  detector  for  the  nlh  projection. 
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Figure  2.  A  schematic  diagram  of  the  DBT  acquisition  geometry  is  shown  (figure  not  to  scale). 
The  x-ray  beam  strikes  point  C  at  the  angle  6n  relative  to  the  normal  to  the  detector.  In  FBP 
reconstruction,  signal  at  C  is  backprojected  to  an  arbitrary  point  E  along  the  incident  ray.  Within 
the  plane  of  the  detector,  backprojection  is  directed  toward  point  F  along  the  angle  Yn  relative  to  the 


795  in'  axis. 
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Figure  3.  Reconstruction  is  performed  with  either  the  ramp  (RA)  filter  alone  or  the  RA  and 
805  spectrum  apodization  (SA)  filters  together.  The  SA  filter  is  a  Hanning  window  function. 
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x  (mm) 

- .  Sinusoidal  Input  (Normalized) 

- Central  Projection,  h  =  70.0  cm 

- Central  Projection,  h  —  oo 


- Central  Projection,  h  =  70.0  cm 

- Central  Projection, 


Sinusoidal  Input  (Normalized) 


- Oblique  Projection,  h  =  70.0  cm 

- Oblique  Projection,  h  =  oo 


- Oblique  Projection,  h  =  70.0  cm 

— -  -  Oblique  Projection,  h  =  qo 


Figure  4.  At  a  distance  y  of  30.0  mm  from  the  chest  wall,  cross  sections  of  signal  in  the  central 
projection  ( n  =  0)  and  the  most  oblique  projection  (n  =  7)  are  plotted  versus  position  x.  In  addition, 
Fourier  transforms  are  shown  versus  frequency.  The  major  Fourier  peaks  do  not  occur  at  the  input 
825  frequency  5.00  lp/mm,  illustrating  the  presence  of  aliasing.  Reducing  the  source-to-COR  distance 
( h )  magnifies  the  input  frequency  projected  onto  the  detector. 
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Figure  5.  Unlike  a  single  projection  (Figure  4),  simple  backprojection  (SBP)  reconstruction  can 
resolve  a  high  frequency  input  oriented  along  the  x  direction.  Applying  filters  to  the  reconstruction 
smoothens  pixilation  artifacts  in  the  spatial  domain  and  reduces  low  frequency  spectral  leakage  in 
the  Fourier  domain.  Reconstructing  with  the  ramp  (RA)  filter  alone  has  the  benefit  of  greater 
835  modulation  than  reconstructing  with  the  RA  and  spectrum  apodization  (SA)  filters  together.  The 
drawback  of  reconstructing  with  the  RA  filter  alone  is  increasing  the  amplitude  of  high  frequency 
spectral  leakage. 
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840  Figure  6.  At  a  reconstruction  depth  (z)  of  50.0  mm,  the  magnitude  of  the  translational  shift  in  the 
U.2  coordinate  of  the  image  [Eq.  (90)]  is  plotted  versus  position  y  measured  perpendicular  to  the 
chest  wall.  In  the  mid  AP/SS  plane  (x  =  0),  translational  shifts  are  minimal  comparing  the  central 
projection  and  an  oblique  projection  (n\  =  0,  m  =  7),  and  are  zero  comparing  the  two  most  oblique 
projections  ( n\  =  -7,  m  =  7).  Increasing  the  magnitude  of  the  distance  x  relative  to  the  mid  AP/SS 
845  plane  yields  a  noticeable  change  in  the  translational  shift. 
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SBP  (x  —  0  mm) 
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Figure  7.  (a)  Within  the  mid  AP/SS  plane  (x  =  0),  SBP  reconstruction  resembles  a  single  projection 
over  the  region  ye  [29.4  mm,  30.6  mm]  for  an  input  frequency  oriented  along  the  y  direction 
860  perpendicular  to  the  chest  wall,  (b)  The  ID  Fourier  transform  of  the  SBP  reconstruction  is  plotted 
versus  frequency  measured  along  the  y  direction.  Within  the  mid  AP/SS  plane  of  a  typical  sized 
breast,  the  major  Fourier  peak  occurs  at  a  frequency  lower  than  the  input  frequency,  5.00  lp/mm. 

(c)  With  x  =  -30.0  mm,  super-resolution  in  an  SBP  reconstruction  is  indeed  achievable  over  the 
region  y  e  [29.4  mm,  30.6  mm],  (d)  For  additional  proof  of  super-resolution  at  x  =  -30.0  mm,  the 
865  major  peak  of  the  corresponding  Fourier  transform  occurs  at  the  input  frequency,  5.00  lp/mm. 
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Central  Projection 


18  15  12  9  6  3  0  -.3 

Distance  .v  (in  mm)  Parallel  to 
Chest  Wall  Side  of  Breast  Support 

Figure  8.  The  central  projection  of  a  bar  pattern  phantom  misrepresents  frequencies  higher  than  the 
detector  alias  frequency,  3.57  lp/mm  for  140  pm  detector  elements.  For  example,  at  4.0  lp/mm, 

870  Moire  patterns  are  present.  At  5.0  lp/mm,  fewer  than  30  line  pairs  are  observed  over  a  6.0  mm 
length. 
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Figure  9.  Unlike  the  central  projection  (Figure  8),  BPF  reconstruction  can  clearly  resolve  high 


frequencies  along  the  x  direction  parallel  to  the  chest  wall  side  of  the  breast  support.  Frequencies 


up  to  6.0  lp/mm  are  resolved  with  no  Moire  patterns  or  other  evidence  of  aliasing. 


890 


895 


900 


48 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


BPF  Reconstruction 
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Figure  10.  Super-resolution  along  the  y  direction  is  analyzed  with  bar  patterns  using  a  BPF 
reconstruction.  The  left  edges  of  the  even  numerals  (“4”  and  “6”)  were  aligned  on  the  mid  AP/SS 
plane  (x  =  0),  and  the  separation  between  4.0  and  5.0  lp/mm  was  positioned  y  =  30  mm  from  the 
905  chest  wall.  At  the  extreme  left  of  the  bar  patterns,  less  line  pairs  are  visible  than  expected, 

illustrating  that  super-resolution  is  not  achievable  near  the  plane  x  =  0.  In  addition,  Moire  patterns 
at  4.0  lp/mm  indicate  that  super-resolution  is  not  possible  too  close  to  the  chest  wall  (y  =  0).  Super¬ 
resolution  is  evident  only  at  positions  sufficiently  displaced  from  the  planes  x  =  0  and  y  =  0;  see  the 
extreme  right  of  the  bar  patterns  at  5.0  and  6.0  lp/mm. 
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(a)  (b) 


Figure  11.  Clinical  images  of  microcalcifications  are  shown.  In  (a),  BPF  reconstruction  is 
915  performed  with  pixels  of  length  140  pm,  and  the  result  is  bilinearly  interpolated  to  produce  a  final 
image  with  pixels  of  length  35  pm  (hence,  four- fold  magnification).  In  (b),  BPF  reconstruction  is 
performed  with  pixels  of  length  35  pm.  As  a  result  of  super-resolution,  image  (b)  displays  the 
structural  morphology  of  the  microcalcifications  more  clearly,  and  even  shows  a  few  small 
microcalcifications  that  cannot  be  discerned  in  (a). 

920 


925 


930 


50 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


REFERENCES 

'Rafferty  EA.  Tomosynthesis:  New  weapon  in  breast  cancer  fight.  Decisions  in  Imaging 
Economics.  2004;  17(4). 

Park  SC,  Park  MK,  Kang  MG.  Super-resolution  Image  Reconstruction:  A  Technical  Overview. 

935  IEEE  Signal  Processing  Magazine.  2003;20(3):21-36. 

'J 

Bushberg  JT,  Seibert  JA,  Edwin  M.  Leidholdt  J,  Boone  JM.  Chapter  10:  Image  Quality.  In:  John  J- 
R,  Snyder  A,  DeGeorge  T,  editors.  The  Essential  Physics  of  Medical  Imaging.  2  ed.  Philadelphia, 
PA:  Lippincott  Williams  &  Wilkins;  2002.  p.  255-91. 

4Acciavatti  RJ,  Maidment  ADA.  Investigating  the  Potential  for  Super-Resolution  in  Digital  Breast 
940  Tomosynthesis.  In:  Pelc  NJ,  Samei  E,  Nishikawa  RM,  editors.  SPIE;  2011;  Lake  Buena  Vista,  FL: 
SPIE;  2011.  p.  79615K-1  -  K-12. 

5Barrett  HH,  Myers  KJ.  Chapter  3:  Fourier  Analysis.  In:  Saleh  BEA,  editor.  Foundations  of  Image 
Science.  New  York,  NY:  John  Wiley  &  Sons;  2004.  p.  95-174. 

6Lee  DL,  Cheung  LK,  Rodricks  B,  Powell  GF.  Improved  imaging  performance  of  a  14  x  17-inch 
945  Direct  Radiography  (TM)  System  using  Se/TFT  detector.  SPIE  Conference  on  Physics  of  Medical 
Imaging;  1998:  SPIE;  1998.  p.  14-23. 

7Jing  T,  Goodman  CA,  Drewery  J,  Cho  G,  Hong  WS,  Lee  H,  et  al.  Amorphous  silicon  pixel  layers 
with  cesium  iodide  converters  for  medical  radiography.  IEEE  Trans  Nucl  Sci.  1994;41(4):903-9. 
8Cowen  AR,  Kengyelics  SM,  Davies  AG.  Solid-state,  flat-panel,  digital  radiography  detectors  and 
950  their  physical  imaging  characteristics.  Clin  Radiol.  2008;63:487-98. 

9Nagarkar  VV,  Gupta  TK,  Miller  SR,  Klugerman  Y,  Squillante  MR,  Entine  G.  Structured  CsI(Tl) 
scintillators  for  x-ray  imaging  applications.  IEEE  Trans  Nucl  Sci.  1998;45(3):492-6. 

10Stewart  J.  Chapter  7:  Techniques  of  Integration.  In:  Pirtle  B,  Green  S,  editors.  Calculus:  Early 
Transcendentals.  5e  ed.  Belmont,  CA:  Brooks/Cole  -  Thomson  Learning;  2003.  p.  474-545. 

955  "Barrett  HH,  Myers  KJ.  Chapter  10:  Energy  Transport  and  Photons.  In:  Saleh  BEA,  editor. 
Foundations  of  Image  Science.  New  York,  NY:  John  Wiley  &  Sons;  2004.  p.  551-630. 

"Barrett  HH,  Myers  KJ.  Chapter  4:  Series  Expansions  and  Integral  Transforms.  In:  Saleh  BEA, 
editor.  Foundations  of  Image  Science.  New  York,  NY:  John  Wiley  &  Sons;  2004.  p.  175-214. 
13Zhao  B,  Zhao  W.  Three-dimensional  linear  system  analysis  for  breast  tomosynthesis.  Medical 
960  Physics.  2008;35(12):5219-32. 

14Mertelmeier  T,  Orman  J,  Haerer  W,  Dudam  MK.  Optimizing  filtered  backprojection 
reconstruction  for  a  breast  tomosynthesis  prototype  device.  In:  Flynn  MJ,  Hsieh  J,  editors.  Physics 
of  Medical  Imaging;  2006;  San  Diego:  SPIE;  2006. 

"Albert  M,  Maidment  ADA.  Linear  response  theory  for  detectors  consisting  of  discrete  arrays. 

965  Medical  Physics.  2000;27(10):2417-34. 

16Acciavatti  RJ,  Maidment  ADA.  An  Analytical  Model  of  NPS  and  DQE  Comparing  Photon 
Counting  and  Energy  Integrating  Detectors.  In:  Samei  E,  Pelc  NJ,  editors.  Physics  of  Medical 
Imaging;  2010;  San  Diego,  CA:  SPIE;  2010.  p.  762201-1  - 1-12. 

17Acciavatti  RJ,  Maidment  ADA.  A  Comparative  Analysis  of  OTF,  NPS,  and  DQE  in  Energy 
970  Integrating  and  Photon  Counting  Digital  X-ray  Detectors.  Medical  Physics.  2010;37(12):6480-95. 
18Carton  A-K,  Bakic  P,  Ullberg  C,  Maidment  ADA.  Development  of  a  3D  high-resolution  physical 
anthropomorphic  breast  phantom.  In:  Samei  E,  Pelc  NJ,  editors.  Medical  Imaging  2010:  Physics  of 
Medical  Imaging;  2010;  San  Diego,  CA:  SPIE;  2010.  p.  762206-1  - -8. 

19Carton  A-K,  Bakic  P,  Ullberg  C,  Derand  H,  Maidment  ADA.  Development  of  a  physical  3D 
975  anthropomorphic  breast  phantom.  Medical  Physics.  201 1;38(2):891-6. 

20Wolbarst  AB.  Chapter  30:  Resolution  and  Magnification.  Physics  of  Radiology.  2nd  ed.  Madison, 
WI:  Medical  Physics  Publishing;  2005.  p.  333-40. 


51 


980 

985 

990 

995 

1000 

1005 

1010 

1015 

1020 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


2 1 

Albert  M,  Beideck  DJ,  Bakic  PR,  Maidment  ADA.  Aliasing  effects  in  digital  images  of  line-pair 

phantoms.  Medical  Physics.  2002;29(8):  1716-8. 

22 

Lanyi  M.  Chapter  7:  Differential  Diagnosis  of  Microcalcifications.  Diagnosis  and  Differential 
Diagnosis  of  Breast  Calcifications.  Berlin:  Springer- Verlag;  1988.  p.  193-231. 

rs-> 

Que  W,  Rowlands  JA.  X-ray  imaging  using  amorphous  selenium:  Inherent  spatial  resolution. 
Medical  Physics.  1995;22(4):365-74. 

24Hajdok  G,  Cunningham  IA.  Penalty  on  the  detective  quantum  efficiency  from  off-axis  incident  x 
rays.  In:  Yaffe  MJ,  Flynn  MJ,  editors.  Medical  Imaging  2004:  Physics  of  Medical  Imaging;  2004; 
San  Diego:  SPIE;  2004.  p.  109-18. 

25Mainprize  JG,  Bloomquist  AK,  Kempston  MP,  Yaffe  MJ.  Resolution  at  oblique  incidence  angles 
of  a  flat  panel  imager  for  breast  tomosynthesis.  Medical  Physics.  2006;33(9):3 159-64. 

26Acciavatti  RJ,  Maidment  ADA.  Calculation  of  OTF,  NPS,  and  DQE  for  Oblique  X-Ray  Incidence 
on  Turbid  Granular  Phosphors.  Lecture  Notes  in  Computer  Science.  2010;6136:436-43. 

27Acciavatti  RJ,  Maidment  ADA.  Optimization  of  phosphor-based  detector  design  for  oblique  x-ray 
incidence  in  digital  breast  tomosynthesis.  Medical  Physics.  201 1;38(1 1):6188-202. 

28  Johns  PC,  Yaffe  MJ.  X-ray  Characterization  of  normal  and  neoplastic  breast  tissues.  Physics  in 
Medicine  and  Biology.  1987;32(6):675-95. 

29Yaffe  MJ.  Chapter  5:  Digital  Mammography.  In:  Beutel  J,  Kundel  HL,  Metter  RLV,  editors. 
Handbook  of  Medical  Imaging  Volume  1  Physics  and  Psychophysics.  Bellingham,  WA:  SPIE  -  The 
International  Society  for  Optical  Engineering;  2000.  p.  329-72. 

on 

Bushberg  JT,  Seibert  JA,  Edwin  M.  Leidholdt  J,  Boone  JM.  Chapter  8:  Mammography.  In:  John  J- 
R,  Snyder  A,  DeGeorge  T,  editors.  The  Essential  Physics  of  Medical  Imaging.  2  ed.  Philadelphia, 
PA:  Lippincott  Williams  &  Wilkins;  2002.  p.  191-229. 

31Berger  MJ,  Hubbell  JH,  Seltzer  SM,  Chang  J,  Coursey  JS,  Sukumar  R,  et  al.  XCOM:  Photon 
Cross  Sections  Database,  [cited  2011  November  10];  Available  from:  http  ://physics .nist.  gov/xcom 
32Johns  HE,  Cunningham  JR.  Chapter  16:  Diagnostic  Radiology.  The  Physics  of  Radiology.  4th  ed. 
Springfield,  IL:  Charles  C  Thomas;  1983.  p.  557-669. 

33Ren  B,  Ruth  C,  Stein  J,  Smith  A,  Shaw  I,  Jing  Z.  Design  and  performance  of  the  prototype  full 
field  breast  tomosynthesis  system  with  selenium  based  flat  panel  detector.  In:  Flynn  MJ,  editor. 
SPIE;  2005;  San  Diego,  CA:  SPIE;  2005.  p.  550-61. 

34Bloomquist  AK,  Yaffe  MJ,  Mawdsley  GE,  Hunter  DM.  Lag  and  ghosting  in  a  clinical  flat-panel 
selenium  digital  mammography  system.  Medical  Physics.  2006;33(8):2998-3005. 

35Siewerdsen  JH,  Jaffray  DA.  A  ghost  story:  Spatio-temporal  response  characteristics  of  an  indirect- 
detection  flat-panel  imager.  Medical  Physics.  1999;26(8):  1624-41. 

36Zhao  W,  DeCrescenzo  G,  Kasap  SO,  Rowlands  JA.  Ghosting  caused  by  bulk  charge  trapping  in 
direct  conversion  flat-panel  detectors  using  amorphous  selenium.  Medical  Physics.  2005;32(2):488- 
500. 

Barrett  HH,  Myers  KJ.  Chapter  12:  Noise  in  Detectors.  In:  Saleh  BEA,  editor.  Foundations  of 
Image  Science.  New  York,  NY:  John  Wiley  &  Sons;  2004.  p.  701-800. 

38Tucker  DM,  Barnes  GT,  Chakraborty  DP.  Semiempirical  model  for  generating  tungsten  target  x- 
ray  spectra.  Medical  Physics.  1991;18(2):21 1-8. 

39Boone  JM,  Fewell  TR,  Jennings  RJ.  Molybdenum,  rhodium,  and  tungsten  anode  spectral  models 
using  interpolating  polynomials  with  application  to  mammography.  Medical  Physics. 

1997;24(12):  1863-74. 

40Blough  MM,  Waggener  RG,  Payne  WH,  Terry  JA.  Calculated  mammographic  spectra  confirmed 
with  attenuation  curves  for  molybdenum,  rhodium,  and  tungsten  targets.  Medical  Physics. 
1998;25(9):  1605-12. 


52 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


1025  41Lauritsch  G,  Harer  WH.  A  theoretical  framework  for  filtered  backprojection  in  tomosynthesis. 
Imaging  Processing;  1998;  San  Diego:  SPIE;  1998.  p.  1127-37. 

42 

Kunze  H,  Haerer  W,  Orman  J,  Mertelmeier  T,  Stierstorfer  K.  Filter  determination  for 
tomosynthesis  aided  by  iterative  reconstruction  techniques.  9th  International  Meeting  on  Fully 
Three-Dimensional  Image  Reconstruction  in  Radiology  and  Nuclear  Medicine;  2007  July  9-13; 
1030  Lindau,  Germany;  2007.  p.  309-12. 

43Bakic  PR,  Zhang  C,  Maidment  ADA.  Development  and  characterization  of  an  anthropomorphic 
breast  software  phantom  based  upon  region-growing  algorithm.  Medical  Physics.  201 1;38(6):3165- 
76. 

44Stewart  J.  Chapter  15:  Multiple  Integrals.  In:  Pirtle  B,  Green  S,  editors.  Calculus:  Early 
1035  Transcendentals.  5e  ed.  Belmont,  CA:  Brooks/Cole  -  Thomson  Learning;  2003.  p.  980-1053. 


53 


TUDAD301D07:  Experimental  and  Theoretical  Validation  of  Breast  Tomosynthesis  Reconstructions  along  Oblique  Planes  |  Browse  -  Medical  Physics 


I  Jump  to  Content 


I  Increase  text  size  I 


flSign  In  View  Cart  Help 


J  Volume/Page  |  Keyword  |  DOI  |  Citation  |  Advanced 

Your  Google  Scholar  search  was  "raymond 
acciavatti  -  experimental  and 
theoretical" 


The  Internationa]  Journal  of  Medical  Physics  Research 


Home 

Browse 

About  T 

Authors  T 

Librarians  T 

Track  ▼ 

Advertisers  T 

Scitation  T 

Special  Publications  T 

MedPhys  Home 

AAPM 

Your  access  to  this  publication  is  provided  through  the  subscription  of  Univ  of  Pennsylvania  Library. 

Medical  Physics  /  Volume  38  /  Issue  6  /  2011  JOINT  AAPM/COMP  MEETING  PROGRAM  /  IMAGING  Previous  Article  |  Next  Article 

SCIENTIFIC  SESSION:  ROOM  301 

Med.  Phys.  38,  3746  (2011);  doi:10. 1118/1. 3613097  ( 1  page ) 

TU  A  301  07:  Experimental  and  Theoretical  Validation  of  Breast 
Tomosynthesis  Reconstructions  along  Oblique  Planes 

R  Acciavatti  and  A  Maidment 

University  of  Pennsylvania ,  Philadelphia ,  PA 

View  Map1?1 


O  ScandiDos 

Delta4®1 "■ 
Tine  Volumetric 
Pre  -Treatment 
^Verification 


The  onlv, 
complete 
QA  solution 
for: 


Abstract 

Full  Text:  Download  PDF  |  View  Cart  ©Alerts  &  Tools  £  Share 

Purpose:  Digital  breast  tomosynthesis  (DBT)  is  a  3D  imaging  modality  in  which 
tomographic  sections  of  the  breast  are  generated  from  a  limited  range  of  x  ray 
projections.  Conventional  practice  is  to  perform  reconstructions  possessing 
pitches  within  the  angular  range  of  the  DBT  scan,  since  the  Central  Slice 
theorem  states  that  Fourier  space  is  sampled  within  double  napped  cones 
(DNCs)  whose  opening  angle  matches  that  angular  range.  This  work 
investigates  the  possibility  for  both  resolution  and  super  resolution  (i.e., 
sub  pixel  resolution)  outside  the  angular  range  of  the  DNCs. Methods:  Because 
the  image  of  an  object  is  translated  in  sub  pixel  detector  element  increments 
with  each  projection,  our  prior  work  has  demonstrated  that  DBT  is  capable  of 
super  resolution.  The  previous  study  assumed  a  reconstruction  plane  parallel 
to  the  breast  support;  our  current  work  analyzes  super  resolution  in  oblique 
reconstruction  planes.  Experimentally,  a  bar  pattern  phantom  was  imaged  with 
a  commercial  DBT  system  using  a  goniometry  stand,  and  reconstruction  was 
performed  in  the  oblique  plane  of  the  bar  patterns.  Clinical  images  of 
microcalcifications  were  similarly  reconstructed  in  various  oblique  planes. 

Subsequently,  an  analytical  framework  for  investigating  super  resolution  in 
oblique  reconstructions  was  developed  by  calculating  the  filtered  backprojection 
(FBP)  reconstruction  of  a  high  frequency  sine  input. Results:  Bar  pattern 
reconstructions  showed  visibility  of  frequencies  both  less  than  and  greater  than 
the  alias  frequency  of  the  detector  at  pitches  well  outside  the  angular  range  of 
the  DBT  scan.  Visibility  of  microcalcifications  did  not  differ  considerably  using 
similar  oblique  reconstruction  planes.  For  analytical  proof  of  super  resolution  in 
oblique  reconstruction  planes,  we  demonstrated  that  FBP  could  properly  resolve 
a  high  frequency  sine  input  whose  Fourier  transform  is  non  vanishing  outside 
the  DNCs  of  frequency  space.  Conclusions:  This  work  provides  a  platform  for 
investigating  super  resolution  in  oblique  reconstruction  planes  whose  pitches 
are  outside  the  angular  range  of  the  DBT  scan. 
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Super-resolution  reconstruction  of  microcalcifications  along  a  plane  parallel  to  the  breast  support 

30°  Pitch 


Super-resolution  reconstruction  of  microcalcifications  along  an  oblique  plane 

In  digital  breast  tomosynthesis  (DBT),  tomographic  sections  of  the  breast  are  generated  from  x-ray 
projections  acquired  over  a  limited  range  of  tube  angles.  Because  the  image  of  an  object  is  shifted  in 
subpixel  increments  with  increasing  projection  angle,  we  have  previously  shown  that  DBT  is  capable  of 
super-resolution  (i.e.,  subpixel  resolution).  This  previous  work  assumed  a  reconstruction  plane  parallel  to 
the  breast  support;  it  is  now  clear  that  super-resolution  is  also  achievable  along  an  oblique  reconstruction 
plane.  A  theoretical  framework  for  investigating  super-resolution  has  been  developed  in  which  the 
reconstruction  of  a  high-frequency  sine  input  is  calculated.  With  this  model,  we  show  that  filtered 
backprojection  can  resolve  the  input  over  a  range  of  angular  pitches.  The  feasibility  of  this  concept  was 
verified  with  a  commercial  DBT  system  (Selenia  Dimensions,  Hologic,  Bedford,  Massachusetts)  and  a 
commercial  prototype  reconstruction  solution  (Briona,  Real  Time  Tomography,  Villanova,  Pennsylvania). 
Using  a  goniometry  stand,  images  of  a  bar  pattern  at  oblique  angles  relative  to  the  breast  support  were 
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acquired  and  successfully  reconstructed.  To  demonstrate  the  clinical  importance  of  these  findings, 
reconstructions  of  microcalcifications  were  performed  along  oblique  planes  whose  pixel  sizes  (28  pm)  are 
much  smaller  than  the  detector  elements  (140  pm).  As  shown,  the  visibility  of  microcalcifications  in  the 
lower  half  of  the  figures  is  not  compromised  by  using  a  30°  pitch.  The  upper  right  calcifications  are  not 
visible  with  the  30°  pitch  as  they  are  out  of  the  reconstruction  plane;  visualization  is  achieved  by  simply 
translating  the  plane  of  reconstruction.  In  conclusion,  this  work  demonstrates  the  feasibility  of  super¬ 
resolution  in  oblique  DBT  reconstructions  and  has  applications  in  the  visualization  of  microcalcifications 
and  other  subtle  signs  of  breast  cancer. 
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Optimization  of  Continuous  Tube  Motion  and  Step-and-Shoot 
Motion  in  Digital  Breast  Tomosynthesis  Systems  with  Patient  Motion 

Raymond  J.  Acciavatti  and  Andrew  D.  A.  Maidment 
University  of  Pennsylvania,  Department  of  Radiology,  3400  Spruce  St.,  Philadelphia  PA  19104 
E-mail:  racci@seas.upenn.edu  and  Andrew.Maidment@uphs.upenn.edu 


ABSTRACT 

In  digital  breast  tomosynthesis  (DBT),  a  reconstruction  of  the  breast  is  generated  from  projections  acquired  over  a 
limited  range  of  x-ray  tube  angles.  There  are  two  principal  schemes  for  acquiring  projections,  continuous  tube  motion 
and  step-and-shoot  motion.  Although  continuous  tube  motion  has  the  benefit  of  reducing  patient  motion  by  lowering 
scan  time,  it  has  the  drawback  of  introducing  blurring  artifacts  due  to  focal  spot  motion.  The  purpose  of  this  work  is  to 
determine  the  optimal  scan  time  which  minimizes  this  trade-off.  To  this  end,  the  filtered  backprojection  reconstruction 
of  a  sinusoidal  input  is  calculated.  At  various  frequencies,  the  optimal  scan  time  is  determined  by  the  value  which 
maximizes  the  modulation  of  the  reconstruction.  Although  prior  authors  have  studied  the  dependency  of  the  modulation 
on  focal  spot  motion,  this  work  is  unique  in  also  modeling  patient  motion.  It  is  shown  that  because  continuous  tube 
motion  and  patient  motion  have  competing  influences  on  whether  scan  time  should  be  long  or  short,  the  modulation  is 
maximized  by  an  intermediate  scan  time.  This  optimal  scan  time  decreases  with  object  velocity  and  increases  with 
exposure  time.  To  optimize  step-and-shoot  motion,  we  calculate  the  scan  time  for  which  the  modulation  attains  the 
maximum  value  achievable  in  a  comparable  system  with  continuous  tube  motion.  This  scan  time  provides  a  threshold 
below  which  the  benefits  of  step-and-shoot  motion  are  justified.  In  conclusion,  this  work  optimizes  scan  time  in  DBT 
systems  with  patient  motion  and  either  continuous  tube  motion  or  step-and-shoot  motion  by  maximizing  the  modulation 
of  the  reconstruction. 

Keywords:  Digital  breast  tomosynthesis  (DBT),  continuous  tube  motion,  step-and-shoot  motion,  patient  motion,  image 
reconstruction,  filtered  backprojection,  modulation,  optimization. 


1.  INTRODUCTION 

Digital  breast  tomosynthesis  (DBT)  is  a  3D  imaging  modality  in  which  tomographic  sections  of  the  breast  are  generated 
from  a  limited  range  of  x-ray  projections.  Preliminary  studies  indicate  that  DBT  has  increased  sensitivity  and  specificity 
for  early  cancer  detection  relative  to  conventional  2D  digital  mammography.1  There  are  two  main  schemes  for  acquiring 
projection  images  in  DBT,  step-and-shoot  motion  and  continuous  tube  motion.  Systems  with  continuous  tube  motion 
have  the  benefit  of  shorter  scan  time  and  thus  less  patient  motion;  the  trade-off  is  increased  blurring  due  to  focal  spot 
motion.  Using  a  prototype  DBT  system  (Hologic  Inc.,  Bedford,  MA),  Ren  et  al.  showed  that  blurring  due  to  focal  spot 
motion  increases  with  height  above  the  breast  support.  At  a  4.0  cm  height,  the  projected  distance  traveled  by  the  focal 
spot  during  a  single  exposure  is  approximately  half  the  detector  element  length.2 

According  to  Zhao,  focal  spot  motion  degrades  the  modulation  transfer  function  (MTF)  of  each  projection  by  sinc^^), 
where  ax  is  the  projected  distance  traveled  by  the  focal  spot  and  fr  is  radial  frequency  perpendicular  to  the  ray  of 
incidence.  Because  focal  spot  motion  has  no  effect  on  noise  power  spectra  (NPS),  the  degradation  in  detective  quantum 
efficiency  (DQE)  is  more  pronounced  than  the  degradation  in  MTF  due  to  the  dependency  of  DQE  on  the  square  of 
MTF.  At  the  alias  frequency  of  5.9  line  pairs  per  millimeter  (lp/mm)  in  a  prototype  system,  Zhao  found  that  focal  spot 
motion  degrades  MTF  and  DQE  by  30%  and  50%,  respectively.3 

In  order  to  minimize  the  blurring  due  to  focal  spot  motion  in  a  system  with  continuous  tube  motion,  Bissonnette  et  al. 
proposed  lengthening  the  scan  time.  They  demonstrated  that  a  39  s  scan  time  effectively  eliminated  image  quality 
degradation  due  to  focal  spot  motion  in  a  prototype  Siemens  Novation™  system.4  Unfortunately,  it  is  not  practical  to 
employ  a  long  scan  time  as  it  permits  greater  patient  motion. 
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Figure  1  illustrates  the  blurring  due  to  patient  motion  in  clinical  images  acquired  with  the  Selenia  Dimensions  system 
(Hologic  Inc.,  Bedford,  MA)  and  reconstructed  with  a  commercial  prototype  backprojection  filtering  algorithm  (Briona, 
Real  Time  Tomography,  Villanova,  PA).  A  small  region  of  interest  (ROI)  with  two  microcalcifications  at  the  height 
22.0  mm  above  the  breast  support  is  shown.  This  depth  was  found  to  minimize  the  motion  of  the  microcalcifications  in 
the  15  individual  backprojections,  and  thus  is  the  height  above  the  breast  support  at  which  the  microcalcifications  are  in 
focus.  By  comparing  the  positions  of  the  microcalcifications  relative  to  a  fixed  marker  (x)  among  all  backprojections  at 
the  22.0  mm  depth,  the  net  displacements  of  the  microcalcifications  in  the  mediolateral  direction  (top-to-bottom  in  the 
figure)  and  the  chest  wall-to-nipple  direction  (left-to-right  in  the  figure)  are  approximately  140  pm  and  280  jam, 
respectively,  corresponding  to  one-  and  two-times  the  length  of  a  detector  element.  With  a  3.7  s  scan  time,  the 
microcalcification  velocities  are  thus  38  jam/s  and  76  jam/s  in  these  two  respective  directions.  The  microcalcifications 
appear  blurry  and  artificially  enlarged  in  the  reconstruction  as  a  result  of  patient  motion. 

Although  image  quality  degradation  due  to  continuous  tube  motion  has  been  modeled  by  many  authors,  no  one  has 
incorporated  patient  motion  into  the  analysis.  Because  these  two  types  of  motion  have  competing  influences  on  whether 
scan  time  should  be  very  long  or  very  short,  one  would  expect  image  quality  to  be  optimized  by  an  intermediate  scan 
time.  For  this  reason,  the  purpose  of  this  work  is  to  determine  the  optimal  scan  time  by  maximizing  the  modulation  of 
the  reconstruction  at  various  frequencies.  To  optimize  step-and- shoot  motion  in  a  similar  fashion,  we  calculate  the  scan 
time  for  which  the  modulation  matches  the  maximum  value  achievable  in  a  comparable  system  with  continuous  tube 
motion.  This  scan  time  provides  a  threshold  below  which  the  benefits  of  step-and-shoot  motion  are  justified. 
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Figure  1.  Backprojections  of  identical  ROIs  of  clinical  tomosynthesis  images  at  a  22.0  mm  height  above  the  breast  support  are 
shown.  This  height  was  found  to  minimize  the  motion  of  the  microcalcifications,  ensuring  that  the  microcalcifications  are  in  focus  in 
the  corresponding  reconstructed  slice.  The  microcalcifications  shift  position  relative  to  a  fixed  marker  (x)  in  the  15  individual 
backprojections.  Such  patient  motion  causes  blurring  and  artificial  enlargement  of  the  microcalcifications  in  the  reconstruction. 


2.  METHODS 

2.1  Acquisition  Geometry  for  Continuous  Tube  Motion  (CTM) 

In  order  to  calculate  the  scan  time  which  optimizes  the  modulation  of  a  DBT  reconstruction  at  various  frequencies,  it  is 
first  necessary  to  model  the  acquisition  geometry.  We  simulate  a  DBT  system  in  which  the  detector  rotates  in  synchrony 
with  the  x-ray  tube  during  the  acquisition  of  the  projections.  As  diagrammed  schematically  in  Figure  2,  the  x-ray  tube 
rotates  within  the  plane  of  the  chest  wall  ( i.e .,  the  xz  plane)  about  the  origin  O,  corresponding  to  midpoint  of  the  chest 
wall  side  of  the  detector.  In  addition,  the  detector  rotates  about  the  y  axis,  with  O  acting  as  the  pivot  point.  At  the  x-ray 
tube  angle  y/  relative  to  the  z  axis,  the  detector  rotation  angle  (7)  is  found  from  the  ratio  y//g,  where  g  is  the  gear  ratio  of 
the  detector.  Positive  directionalities  of  y/  and  y  are  defined  as  those  presented  in  Figure  2,  and  in  the  limit  g  — >  00,  a 
stationary  detector  can  be  recovered  (y  —>  0). 

In  a  system  with  continuous  tube  motion  at  a  constant  angular  velocity  co ,  each  projection  is  acquired  over  the  exposure 
time  t  as  the  tube  is  swept  over  the  angular  extent  T*  =  cot.  For  the  nth  projection,  the  x-ray  tube  arc  is  centered  about  the 
angle  \f/n  =  nAy/,  so  that  the  x-ray  tube  angle  varies  between  1//  =  y/n  +  TV2  and  y/  =  y/n-  TV2  during  the  exposure  time  r. 
In  the  literature,  y/n  is  often  termed  the  nominal  projection  angle  and  Ay/  the  angular  spacing  between  projections.5  With 
an  odd  number  of  N  total  projections,  the  index  n  varies  between  -(N  -  l)/2  and  ( N  -  l)/2,  and  the  special  case  n  =  0 
defines  the  central  projection.  Denoting  the  total  scan  time  as  Tt ,  the  total  angular  range  of  the  x-ray  tube  motion  can  be 
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Figure  2.  A  diagram  of  the  acquisition  geometry  is  shown  (not  to  scale).  The  attenuation  coefficient  of  the  input  object  varies 
sinusoidally  along  the  x  direction.  To  model  patient  motion,  the  input  object  has  velocity  v  at  the  angle  f  relative  to  the  x  direction. 

written  as  coTt ,  or  equivalently,  as  the  difference  between  the  initial  x-ray  tube  angle  ((///)  and  the  final  x-ray  tube  angle 

<>/)• 


coT,  =w,~vf  = 


N- 1 


A  ^ 


N- 1 


A  y/- 
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=  (N-l)Ay/  +  x¥ 


(i) 


Substituting  co  =  Wr  in  the  left-hand  side  of  Eq.  (1),  the  angular  sweep  of  the  x-ray  tube  over  the  exposure  time  r  can  be 
expressed  in  terms  of  the  total  scan  time  Tt  instead  of  the  tube’s  angular  velocity  co. 


x¥- 


r  N- 1^ 


tAxj/ 


(2) 


This  formula  for  ¥  is  useful  as  the  total  scan  time  is  more  directly  measurable  than  the  tube’s  angular  velocity. 


2.2  Detector  Signal  for  Sinusoidal  Input  to  CTM  System 

A  framework  for  investigating  tube  motion  and  patient  motion  in  DBT  is  now  developed  by  calculating  the  modulation 
of  the  reconstruction  of  a  sinusoidal  input.  Accordingly,  suppose  that  a  thin  rectangular  plate  with  its  long  axis  parallel 
to  the  breast  support  possesses  a  linear  attenuation  coefficient  ju(x,  z)  which  varies  sinusoidally  with  position  x.  Although 
an  actual  input  to  a  clinical  breast  imaging  system  would  be  3D,  a  2D  construct  is  a  useful  tool  for  simulating 
measurements  in  the  plane  of  the  chest  wall.  The  extension  of  this  framework  to  measurements  made  perpendicular  to 
the  chest  wall  is  reserved  for  future  work. 


As  shown  in  Figure  2,  the  rectangular  plate  is  positioned  between  z  =  z(y/) -s/2  and  z  =  z(y/)  +  s  12  ,  where  z{y/)  is 
the  central  height  of  the  plate  above  the  detector  at  the  x-ray  tube  angle  yj  and  s  is  the  plate’s  thickness.  The  height  z(y/) 

is  taken  to  be  dependent  upon  the  x-ray  tube  angle  yj  in  order  to  model  the  presence  of  patient  motion.  For  an  input 
frequency  f0,  the  attenuation  coefficient  may  be  written 


/j(x , z)-C  •  cos [2n f0  [x  -  x(^)])  • 


rect 


z-z{\f/) 


rect(M) 


1  ,|  u  |<  1/2 

[0  ,  |  w  |>  1  /  2  ’ 


(3) 
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where  C  is  the  amplitude  of  the  sinusoidal  waveform  and  x(y/)  is  its  translational  shift  along  the  x  direction.  Provided 

that  \z  -  z0 1  <  s,  the  Fourier  transform  of  Eq.  (3)  along  the  x  direction  is  a  linear  sum  of  delta  functions5  which  peak  at  the 
frequencies /=  ±f0.  Typically,  only  the  positive  frequency /=  +fo  is  of  interest  in  an  experimental  measurement.  Hence, 
although  it  is  non-physical  for  a  linear  attenuation  coefficient  to  vary  between  negative  and  positive  values,  formulating 
p(x ,  z)  by  Eq.  (3)  is  helpful  for  a  thought  experiment  in  the  reconstruction  of  a  single  input  frequency. 

In  Figure  2,  the  displacements  x\(y/)  and  x2(y/)  determine  the  entrance  and  exit  points  of  the  x-ray  beam  through  the  sine 
plate  for  an  arbitrary  incident  point  on  the  detector  at  a  distance  r  from  O.  Following  our  previous  work,5  x\(y/)  and  x2{y/) 
can  be  written  as  x\(y/)  =  p(y/)-r  -  k+(y/)  and  x2(y/)  =  p(y/)-r  -  k~{y/),  where  p{y/)  =  cos[y(^)]  +  sin[y(^/)]tan[0(^)  +  y(y/)]  and 
A±  (y/)  =  [z<y )  ±  s  /  2]  •  tan \0(y/)  +  y(y/)] .  The  expression  for  the  incident  angle  relative  to  the  normal  to  the  detector  also 
follows  from  our  previous  work 


0(y/)  =  -y(y/)  +  arctan 


r  h  sin  y/  +  r  cos  [y  {y/)\  ^ 
h  cos  y/ -r  sin  [y(^)]  J  ’ 


(4) 


where  h  is  the  source-to-origin  distance  (Figure  2).  Total  attenuation  Ap(y/)  recorded  by  the  x-ray  converter  at  the  tube 
angle  y/  may  now  be  calculated  by  integrating  the  attenuation  coefficient  of  the  sine  plate  over  the  path  length  C(y/). 

=  f  ?/)  fids  =  £2]  C  • cos  {l7rfo  ix  -  ^(^)])  • csc  +  r(v)]dx  (5) 

In  order  to  calculate  the  total  attenuation  Ap{n)  recorded  by  the  x-ray  converter  for  the  nth  projection,  one  must  integrate 
Ap(y/)  over  the  angular  arc  swept  by  the  x-ray  tube  during  the  exposure  time  r. 


Aju(n)  = 


C-Jr  c ' cos  (2' 71  ^  _  *(^)])  • csc  [oty) + dxd[y 

sin  (2^/0  [_p(w)  -r-y x(^)]) 
-  sin  ■  r  -  A+  (ip)  -  x (^)]) 


C  rv'n+'V! 2 


2^/o  ^n-V!2 


CSC 


[0(vs)+r(v)]- 


dy/ 


(6) 

(7) 


Eq.  (7)  provides  an  expression  for  signal  intensity  versus  position  r  along  the  detector,  assuming  that  the  detector  is  non- 
pixilated  and  possesses  a  modulation  transfer  function  (MTF)  of  unity  at  all  frequencies.  An  amorphous  selenium  (a-Se) 
photoconductor  operated  in  drift  mode  is  a  good  approximation  for  a  detector  with  these  properties.6 

Total  attenuation  for  the  nth  projection  can  now  be  simplified  using  a  sum-to-product  trigonometric  identity  for  real 
numbers  a  and  /?;  namely,  sina  -  sin/?  =  2cos[(a  +  /?)/2]sin[(a  -  /?)/2] 

A/u(n)  =  sec  [6(y/)  +  y(i//)] sine  [ef0  tan [<9(^)  +  /(^)]] 

5  (8) 

•  cos  [2 7rf0  [ p(y y)  •  r  -  z{y/)  •  tan  [6(y/)  +  y(y/)]  -  x(y/)~]  J  dy/ 

where  sinc(i/)  =  sin  (mi)/ (mi).  Because  it  is  difficult  to  perform  the  integration  in  Eq.  (8)  in  closed  form,  it  is  necessary  to 
use  approximation  techniques.  One  such  method  is  the  midpoint  formula.  For  the  nth  projection,  the  angular  sweep  of 
the  x-ray  tube  can  be  divided  into  K  intervals  between  y/=y/n  +  TV2  and  y/=  y/n-  TV2.  The  tube  angle  at  the  midpoint  of 
the  kth  interval  is 


Vkn=Vn~  —  \  !■ 


2k~\ 

K 


T 

n - 

2 


N- 1 

\Tt~r  j 


2k -l 
K 


Ay/  ,  k  e  N  . 


Eq.  (8)  can  now  be  evaluated  by  averaging  the  integrand  over  each  of  the  K  intervals  in  the  limit  of  infinite  K 


(9) 
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Cs  K 

A/.i(  n )  =  lim  —  ^  scc(  Qkn  +^„)sinc[f/0  tan(6>,„  +  ykn)]cos\_2/rfa[Pknr-~zkn  tan  (9kn  +ykn)~xkn]\, 

^  k= 1 


(10) 


where  6km  ykn,  and  pkn  are  calculated  by  evaluating  0(y/),  y(y/),  and  p(y/)  aty/  =  y/kn.  In  Eq.  (10),  the  displacements  xkn  and 
z bn  determine  the  position  of  the  sine  plate  at  the  time  point  Tkn.  The  special  case  Tkn  =  0  is  defined  to  occur  at  the  x-ray 
tube  angle  yj  =  0,  so  that  the  scan  time  occurs  between  the  time  points  -TJ2  and  TJ 2. 


T  Wh, 

*"  co  'P 


T 

■<T,  <-L 

kn  2 


(11) 


To  model  the  presence  of  patient  motion,  the  sine  plate  is  taken  to  have  constant  velocity  v  at  the  angle  C  relative  to  the  x 
direction  (Figure  2).  The  displacements  x^  and  zkn  can  be  written  in  terms  of  the  velocity  components  and  v2  as 


*kn  -  X0  +  TknVx  -  X0  + 


T.-Awh.w&G 

N-\)  Ay/ 


=  Z0+7i„Vz  =Z0  + 


T.-AWh,^ 

N-\)  Ay/ 


where  x0  and  z0  are  positions  which  determine  the  location  of  the  sine  plate  at  the  x-ray  tube  angle  yj  =  0. 


(12) 


In  a  digital  detector,  the  a- Se  x-ray  converter  is  placed  in  electrical  contact  with  a  large  area  plate  of  amorphous  silicon 
(<2-Si)  in  which  a  thin-film  transistor  (TFT)  array  samples  the  total  attenuation  in  pixels  (i.e.,  detector  elements).  Using 
Eq.  (10),  the  logarithmically- transformed  signal  in  the  mth  detector  element  for  the  nth  projection  is 


(13) 


Detector  elements  are  taken  to  be  centered  on  r  =  ma,  and  the  detector  element  containing  O  is  the  one  corresponding  to 
m  =  0.  Because  the  incident  angle  varies  minimally  within  each  detector  element,  the  integration  in  Eq.  (13)  can  be 
evaluated  by  approximating  the  incident  angle  6kn  with  its  value  at  the  centroid  of  the  detector  element.  Thus 


sec(^,„  +  )  • sinc  [g/0  tan(^„,„ + r* )] 

2  nf,KPhm 


sin  [2 (aPkmn  (m  + 1  /  2)  -  ~zkn  tan {0kmn  +  yhn )  -  Sck„ )] 
-  sin  [2  71  f0  (aPkmn  (m  - 1  /  2)  -  zfa  tan  {0kmn  +yj-xhl )] 


(14) 


where  Qkmn  and  pkmn  are  calculated  by  evaluating  0kn  and  pkn  at  r  =  ma.  Eq.  (14)  can  be  simplified  further  by  using  the 
sum-to-product  trigonometric  identity  described  previously. 


Cs  K 

VP(m,  n)  =  lim  ■ —  Z  sec(6'*»»  +  ft, )  • sinc  [£fo  tan(^„  +  )] '  sine  [aPkmnf0  ] 


•cos 


[2^/0  [maPkmn  - zkn  tan(6»ta„  +ykn)~ xkn ]] 


(15) 


2.3  Detector  Signal  for  Step-and-Shoot  Motion  (SSM) 

In  a  similar  fashion,  detector  signal  for  a  system  with  step-and-shoot  motion  can  be  calculated.  All  expressions  between 
Eqs.  (4)  and  (8)  continue  to  hold,  so  that  the  total  attenuation  recorded  by  the  x-ray  converter  for  the  nth  projection  is 

Aju(n)  =  Cs  sec(<9„  +  yn  )sinc  [sft  tan (6>„  +  yn )]  cos  [2 nf0  [Pr  -  :(T)  •  tan(6>„  +yn)-  x(J)]]  .  (16) 

Because  the  x-ray  tube  angle  yj  remains  constant  during  a  single  projection  in  a  system  with  step-and-shoot  motion,  it  is 
acceptable  to  simplify  Eq.  (8)  by  evaluating  0(y/),  y(y/),  and  p(y/)  at  the  nominal  projection  angle  y/  =  y/n  as  denoted  by  the 
parameters  0„,  ym  and  pm  respectively.  Although  the  object  coordinates  x(y /)  and  z{y/)  are  dependent  upon  the  x-ray 
tube  angle  yj  in  a  system  with  continuous  tube  motion,  the  same  coordinates  are  now  dependent  upon  time  T. 
Consequently,  the  integral  over  yj  in  Eq.  (8)  can  be  replaced  by  an  integral  over  T  to  take  into  account  the  presence  of 


5 


DRAFT:  The  final  publication  may  differ  slightly  from  this  version. 


patient  motion.  The  integration  limits  in  Eq.  (16)  are  thus  the  initial  and  final  time  points  at  which  the  tube  emits  x  rays 
during  a  single  projection,  where  Tn  is  the  central  time  point  of  the  projection  of  duration  r  and  AT  is  the  time  difference 
between  consecutive  projections.  It  can  be  shown  that  Tn  =  n(r  +  AT)  and  A T=  ( Tt  -  Nz)/(N -  1).  Substituting  the  object 
coordinates  x(T)  =  x0+vT  cos£  and  z(T )  =  z0  +  vT  sin£  into  Eq.  (16)  yields 

Aju(n)  =  Cs  sec(  0  +  yn  )sinc  [sf0  tan(6>„  +  yn )]  sine  \  f0zv  [cos  ^  +  sin(0  tan(<9„  +  yn )]] 

r  r  m  (17) 

•  cos  |_2 7rf0  \pnr-  (x0  +  vTn  cos  Q  -  (z0  +  vTn  sin  Q  tan(6>„  +  yn )]  J 

By  integrating  A/u(n )  over  the  detector  element  length  a  [Eq.  (13)],  the  logarithmically-transformed  signal  in  the  mth 
detector  element  for  the  nth  projection  can  be  determined 

Vu(m,n)  =  Cs  sec(dmn  +  yn  )sinc  [sf0  tan(  6nm  +  yn )]  sine  [/0rv  [cos  6,"  +  sin(  4  )  tan(  dmn  +^)]]sinc[/0apm„] 

•  cos [2 nf0  [mapmn  -  (x0  +  vTn  cos  Q  -  (z0  +  vTn  sin  Q  tan(6>mn  +  yn )]] 

where  6mn  and  pmn  are  calculated  by  evaluating  0n  and  pn  at  r  =  ma. 


2.4  Filtered  Backp rejection  (FBP)  Reconstruction 

The  attenuation  coefficient  can  now  be  reconstructed  using  filtered  backprojection  (FBP).  From  our  previous  work,5  the 
FBP  reconstruction  for  an  infinitesimally  fine  ( i.e .,  non-pixilated)  reconstruction  grid  is  determined  from  the  expression 


//FBP  (X,  z)  =  T  Yj  Z  n)  ■ 


(/)(t)*  rect 

v 


( t  sec  0mn  -  ma 


\t=xcos{yn  +0mn  )+zsin(yn  +6mrl ) 


(19) 


where  //FBp  is  the  reconstructed  attenuation  coefficient  and  *  is  the  convolution  operator.  The  reconstruction  filter  $7) 
follows  from  linear  systems  theory  for  DBT.  A  ramp  (RA)  filter,  given  by  J/f  in  the  Fourier  domain,  is  first  applied  to 
reduce  the  low  frequency  detector  response.  Since  noise  tends  to  occur  at  high  frequencies,  a  spectrum  apodization  (SA) 
filter  is  also  used;  following  Zhao’s  approach,3  we  apply  a  Hanning  window  function  as  the  SA  filter.  In  the  Fourier 
domain,  the  filters  are  truncated  at  the  frequencies  /=  ±£  and  the  net  filter  is  the  product  of  the  RA  and  SA  filters.  As 
shown  in  our  previous  work,5  the  net  filter  can  be  calculated  in  closed  form  using  the  inverse  Fourier  transform. 


3.  RESULTS 

Reconstructions  are  now  simulated  for  a  Selenia  Dimensions  system  with  15  projections  acquired  at  an  angular  spacing 
(Ay/)  of  1.07°,  assuming  C  =  1.0  mm"1,  h  =  70.0  cm,  s  =  0.50  mm,  and  a  =  140  pm.  At  the  x-ray  tube  angle  y/  =  0,  the 
centroid  of  the  sine  plate  is  taken  to  coincide  with  the  midpoint  of  the  chest  wall  side  of  a  50.0  mm  thick  breast.  With 
the  breast  support  positioned  25.0  mm  above  the  origin  of  the  detector,  the  x0  and  z0  coordinates  of  the  input  object  are 
therefore  0  and  50.0  mm,  respectively. 

3.1  Effect  of  Continuous  Tube  Motion  on  Modulation 

The  effect  of  continuous  tube  motion  on  modulation  is  analyzed  in  Figure  3(a)  by  first  simulating  a  system  with  no 
patient  motion.  For  a  30.0  ms  exposure  time,  corresponding  to  the  mean  value  of  r  for  the  Selenia  Dimensions  system, 
reconstructions  of  the  frequency  2.0  lp/mm  are  calculated.  As  expected,  the  modulation  increases  with  scan  time.  For 
example,  with  scan  times  of  2.0,  3.0,  and  4.0  s,  the  modulation  attains  the  values  0.59,  0.70,  and  0.74,  respectively.  This 
trend  arises  because  the  tube’s  angular  sweep  ¥  during  a  single  projection  decreases  with  scan  time  [Eq.  (2)]. 

For  the  same  DBT  system,  the  modulation  in  the  reconstruction  is  also  studied  as  a  function  of  the  input  frequency  f0 
[Figure  3(c)].  At  low  frequencies,  the  modulation  increases  linearly  from  zero,  following  the  ramp  filter.  At  higher 
frequencies,  the  spectrum  apodization  filter  and  the  MTF  of  the  detector  sampling  process  reduce  the  modulation, 
countering  the  ramp  filter;  hence,  there  is  an  intermediate  frequency  at  which  the  modulation  is  maximized.  This 
frequency  dependence  of  the  modulation  matches  Zhao’s  formulation  of  in-plane  MTF  in  DBT  reconstructions,3 
providing  a  built-in  check  on  the  validity  of  Figure  3(c).  Like  Zhao,  we  plot  the  modulation  over  a  frequency  range 
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Continuous  Tube  Motion,  No  Patient  Motion  Step-and-Shoot  Motion,  Patient  Motion 
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Figure  3.  Although  modulation  increases  with  scan  time  in  a  system  with  continuous  tube  motion  and  no  patient  motion,  the  opposite 
trend  holds  in  a  system  with  step-and-shoot  motion  and  patient  motion. 

spanning  at  least  one  zero  of  the  MTF  of  the  sampling  process  in  the  detector.7  The  filter  truncation  frequency  (£)  of  2a  1 
(14.3  lp/mm)  is  simulated,  corresponding  to  the  second  zero  of  the  detector  sampling  MTF  [sine (a/)].  Figure  3(c)  shows 
that  increasing  the  scan  time  increases  the  modulation  of  the  reconstruction,  and  thus  generalizes  the  trend  presented  in 
Figure  3(a)  to  all  frequencies.  By  contrast,  increasing  the  exposure  time  decreases  the  modulation  of  the  reconstruction. 

In  addition,  Figure  3(c)  demonstrates  that  the  modulation  of  the  reconstruction  may  possess  zeros  at  frequencies  different 
from  those  of  the  detector  sampling  MTF,  whose  first  and  second  zero  are  a 4  (7.1  lp/mm)  and  2a'1  (14.3  lp/mm). 
Increasing  the  exposure  time  decreases  these  additional  zeros,  while  increasing  the  scan  time  increases  the  zeros.  For  a 
50.0  ms  exposure  time,  the  first  zero  not  equivalent  to  a'1  or  2a'1  is  3.0,  4.5,  and  6.0  lp/mm  for  2.0,  3.0,  and  4.0  s  scan 
times  (5.0,  7.6,  and  10.1  lp/mm  for  a  30.0  ms  exposure  time).  The  appendix  shows  that  the  formula  for  these  zeros 
follows  from  the  calculation  of  the  MTF  of  focal  spot  motion.  These  zeros  place  an  important  limit  on  the  resolution  of 
the  system. 

3.2  Effect  of  Patient  Motion  on  Modulation 

In  Figure  3(b),  the  effect  of  patient  motion  on  modulation  is  investigated  by  considering  a  system  with  step-and-shoot 
motion.  With  an  exposure  time  of  30.0  ms  and  an  object  velocity  of  60.0  pm/s  oriented  along  the  x  direction,  Figure  3(b) 
demonstrates  that  modulation  decreases  with  scan  time.  The  object  velocity  considered  in  Figure  3(b)  is  comparable  to 
the  value  observed  in  Figure  1  showing  clinical  images  of  microcalcifications.  The  modulation  attains  the  values  0.69, 
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Effect  of  Angle  on  Modulation 


Angle  C Relative  to  x  Direction 
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Effect  of  Scan  Time  on  Modulation  (C=  45°) 
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Figure  4.  (a)  By  increasing  the  angle  f  of  patient  motion  relative  to  the  x  direction,  the  modulation  of  the  reconstruction  increases. 

(b)  With  continuous  tube  motion  (CTM)  and  patient  motion  occurring  simultaneously,  modulation  is  optimized  by  an  intermediate 
scan  time.  By  contrast,  modulation  is  maximized  by  a  short  scan  time  with  step-and-shoot  motion  (SSM).  (c)  The  dependency  of  the 
optimal  CTM  scan  time  on  object  velocity  (v),  exposure  time  (r),  and  frequency  (f0)  is  investigated,  (d)  The  optimal  CTM  scan  time  is 
larger  with  patient  motion  oriented  at  a  45°  angle  relative  to  the  x  direction  than  a  0°  angle  [Figure  4(b)]. 

0.61,  and  0.52  for  2.0,  3.0,  and  4.0  s  scan  times,  respectively.  This  dependency  of  modulation  on  scan  time  is  expected 
since  the  net  displacement  of  the  object  in  the  x  direction,  given  by  vTt cos£  increases  with  scan  time. 

In  Figure  3(d),  modulation  is  plotted  versus  f0  for  a  30.0  ms  exposure  time.  Figure  3(d)  shows  that  the  modulation 
decreases  with  scan  time  over  all  frequencies,  thus  generalizing  the  trend  shown  in  Figure  3(c).  Unlike  a  system  with 
continuous  tube  motion  and  no  patient  motion,  modulation  varies  minimally  with  exposure  time  in  a  system  with  step- 
and-shoot  motion  and  patient  motion.  This  finding  arises  because  the  object  displacement  between  projections  is 
significantly  greater  than  the  corresponding  motion  during  the  exposure  time  of  an  individual  projection. 

It  is  also  demonstrated  in  Figure  3(d)  that  the  frequency  corresponding  to  the  first  zero  of  the  modulation  may  be  less 
than  that  of  the  detector  sampling  MTF.  This  zero  follows  from  the  point  spread  function  (PSF)  of  patient  motion,  which 
is  a  rectangle  function  whose  width  is  the  object  displacement  during  the  scan  time.  Accordingly,  the  MTF  of  patient 
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motion  is  sinc(vr/cosf),  and  the  first  zero  is  (vTt)AsecC  With  scan  times  of  2.0,  3.0,  and  4.0  s,  the  zeros  of  the  MTF  of 
patient  motion  are  8.3,  5.6,  and  4.2  lp/mm,  respectively. 

In  the  same  system,  Figure  4(a)  studies  the  dependence  of  modulation  on  the  directionality  of  patient  motion.  At  any 
fixed  frequency,  the  modulation  increases  with  the  angle  C  relative  to  the  x  direction.  In  addition,  the  difference  in 
modulation  comparing  any  two  scan  times  is  minimized  as  C  increases.  Modulation  is  virtually  independent  of  scan  time 
if  the  patient  motion  is  oriented  along  the  z  direction  (C=  90°). 

3.3  Optimization  of  Scan  Time 

With  both  continuous  tube  motion  and  patient  motion  occurring  simultaneously,  there  is  a  trade-off  in  the  benefits  of 
long  and  short  scan  time,  and  hence  modulation  is  maximized  by  an  intermediate  scan  time  [Figure  4(b)].  For  example, 
with  object  velocities  (v)  of  30.0  and  60.0  pm/s  oriented  along  the  x  direction,  the  optimal  scan  times  for  continuous  tube 
motion  are  3.3  and  2.4  s,  respectively,  assuming  an  input  frequency  of  2.0  lp/mm  and  an  exposure  time  of  30.0  ms. 
Figure  4(c)  demonstrates  that  this  optimal  scan  time  decreases  with  object  velocity  (v)  and  increases  with  exposure  time 
(r).  Exposure  times  between  30.0  ms  (the  mean  value  of  the  Selenia  Dimensions  system)  and  50.0  ms  (the  maximum 
value  of  the  system)  are  considered.  Figure  4(c)  also  shows  that  the  optimal  scan  time  is  frequency  dependent. 

For  any  fixed  scan  time,  step-and-shoot  motion  (SSM)  yields  greater  modulation  than  continuous  tube  motion  (CTM). 
To  optimize  a  step-and-shoot  system,  one  may  calculate  the  scan  time  giving  the  same  modulation  as  the  highest 
achievable  with  continuous  tube  motion  [Figure  4(b)].  This  scan  time  provides  a  threshold  below  which  the  use  of  step- 
and-shoot  motion  is  justified.  For  example,  with  30.0  pm/s  patient  motion,  an  SSM  scan  time  of  4.7  s  yields  the  same 
modulation  as  the  optimal  CTM  scan  time  of  3.3  s.  With  60.0  pm/s  patient  motion,  the  analogous  SSM  and  CTM  scan 
times  are  3.3  and  2.4  s. 

The  dependence  of  the  optimal  CTM  scan  time  on  the  directionality  of  patient  motion  is  investigated  in  Figure  4(d)  by 
considering  patient  motion  along  a  45°  angle  relative  to  the  x  direction.  The  optimal  scan  time  for  continuous  tube 
motion  is  larger  in  Figure  4(d)  than  in  Figure  4(b)  with  £  =  0°.  For  example,  with  an  object  velocity  (v)  of  30.0  pm/s,  the 
optimal  CTM  scan  times  for  f  =  0°  and  45°  are  3.3  and  3.8  s,  respectively  (2.4  and  2.8  s  for  v  =  60.0  pm/s).  In  addition, 
with  C=  45°,  there  is  a  broader  range  of  scan  times  for  which  the  modulation  is  within  the  limit  of  resolution  of  the 
system,  which  is  often  taken  to  be  0.05.  For  example,  with  an  object  velocity  of  60.0  pm/s  and  a  scan  time  of  7.5  s,  the 
modulation  with  £=  45°  is  0.26  (resolvable),  yet  the  modulation  with  C=  0°  is  0.03  (not  resolvable). 


4.  DISCUSSION  AND  CONCLUSION 

To  our  knowledge,  this  work  is  the  first  to  model  both  continuous  tube  motion  and  patient  motion  in  DBT,  and  as  a 
result,  to  develop  a  technique  which  optimizes  scan  time.  Continuous  tube  motion  and  patient  motion  have  competing 
influences  on  scan  time;  we  show  that  the  modulation  of  the  reconstruction  is  optimized  by  an  intermediate  scan  time. 

In  Figure  4,  it  is  demonstrated  that  continuous  tube  motion  and  step-and-shoot  motion  have  nearly  identical  modulation 
in  systems  with  very  long  scan  time.  For  example,  with  a  2.0  lp/mm  input  and  a  30.0  ms  exposure  time,  one  can  show 
that  the  relative  difference  in  modulation  between  the  two  systems  does  not  exceed  1.0%  for  scan  times  of  10  s  or  more. 
This  result  holds  regardless  of  the  object  velocity  studied  in  Figure  4.  For  this  reason,  there  is  effectively  no  difference 
in  image  quality  between  the  two  systems  if  both  operate  at  the  same,  very  long  scan  time;  patient  motion  is  a  much 
more  significant  cause  of  image  quality  degradation  than  continuous  tube  motion  at  these  scan  times. 

To  minimize  patient  motion  in  DBT,  the  system  should  have  a  short  scan  time  comparable  to  2D  digital  mammography. 
In  systems  with  continuous  tube  motion,  the  drawback  of  lowering  the  scan  time  is  substantially  reducing  the  modulation 
relative  to  an  analogous  step-and-shoot  system  (Figure  4).  Since  systems  with  step-and-shoot  motion  tend  to  have  longer 
scan  times  than  those  with  continuous  tube  motion  due  to  mechanical  considerations,  Figure  4  demonstrates  that  it  is  still 
possible  to  operate  a  system  with  continuous  tube  motion  at  a  scan  time  yielding  superior  image  quality  relative  to  a  step- 
and-shoot  system.  An  additional  benefit  of  continuous  tube  motion  might  include  eliminating  microphonic  vibrations 
during  the  exposure  time  of  each  projection. 
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In  conventional  2D  digital  mammography,  it  has  been  demonstrated  that  the  scan  time  should  be  less  than  2.0  s  to 
minimize  patient  motion.  Currently,  no  such  guideline  for  DBT  has  been  developed.  In  order  to  optimize  scan  time  in 
systems  with  continuous  tube  motion,  our  simulations  considered  object  velocities  between  30.0  and  60.0  jam/s.  These 
velocities  were  chosen  to  be  comparable  to  values  seen  in  clinical  images  of  microcalcifications  presented  in  Figure  1. 
Although  Figure  1  illustrates  motion  of  microcalcifications  in  one  clinical  data  set,  the  case  is  not  necessarily 
representative  of  the  most  significant  extent  of  patient  motion,  and  thus  additional  cases  should  be  considered  in 
developing  guidelines  for  scan  time  in  DBT. 

In  systems  with  continuous  tube  motion,  this  paper  demonstrates  that  the  modulation  of  the  reconstruction  may  be  zero  at 
frequencies  smaller  than  the  zeros  of  the  MTF  of  detector  sampling.  This  finding  has  important  implications  on  the 
visibility  of  high  frequencies  in  DBT,  which  was  the  subject  of  our  prior  work.5  In  performing  reconstructions  on  a  grid 
whose  pixel  size  is  much  smaller  than  the  detector  elements,  we  have  previously  demonstrated  that  DBT  is  capable  of 
super-resolution.  Although  the  alias  frequency  of  the  Selenia  Dimensions  detector  is  3.6  lp/mm,  reconstructions  of  bar 
patterns  have  clearly  shown  higher  frequencies  due  to  super-resolution.  In  this  paper,  Figure  3(c)  demonstrates  that  the 
ability  to  achieve  super-resolution  is  influenced  by  focal  spot  motion,  as  the  zeros  of  the  modulation  in  a  system  with 
continuous  tube  motion  vary  with  exposure  time.  A  formula  for  these  zeros  is  derived  in  the  Appendix.  Using  high 
contrast  bar  patterns,5  Figure  1  in  our  previous  work  demonstrated  visibility  of  6.0  lp/mm  in  the  Selenia  Dimensions 
system  with  a  3.7  s  scan  time  and  a  30.2  ms  exposure  time.  For  these  settings,  it  follows  from  the  Appendix  that  the 
modulation  possesses  a  zero  at  9.4  lp/mm.  It  is  worth  noting  that  if  the  same  measurements  were  taken  with  the  system’s 
maximum  exposure  time  of  50.0  ms,  the  limiting  resolution  of  focal  spot  motion  would  be  5.6  lp/mm,  and  one  would  not 
expect  the  same  6.0  lp/mm  bar  patterns  to  be  successfully  resolved. 

There  has  been  recent  interest  in  acquiring  DBT  images  with  less  compression  than  conventional  2D  mammography.8 
The  purpose  of  the  reduced  compression  is  to  spread  out  tissues  and  thus  improve  resolution  in  the  z  direction 
perpendicular  to  the  breast  support.  This  work  argues  against  reduced  compression,  since  it  inherently  leads  to  greater 
patient  motion  and  would  be  expected  to  degrade  the  modulation  of  the  reconstruction.  In  systems  with  long  scan  times, 
the  need  for  full  compression  is  particularly  evident  because  the  net  object  displacement  should  increase  with  scan  time. 

Some  of  the  limitations  of  this  study  and  directions  for  future  modeling  are  now  noted.  Future  work  should  more 
carefully  model  the  MTF  degradation  due  to  non-normal  x-ray  incidence9, 10  as  well  as  the  finite  size  of  the  focal  spot.11 
In  addition,  the  presence  of  noise  at  various  radiation  dose  levels  could  be  simulated.  Although  this  work  implicitly 
assumes  a  high  contrast  input  frequency  whose  visibility  is  independent  of  dose,  future  studies  should  demonstrate  how 
the  optimization  of  scan  time  is  influenced  by  noise12  at  various  dose  levels  for  low  contrast  input  frequencies.  Because 
the  attenuation  coefficient  of  the  input  object  is  energy  dependent,  poly  energetic  x-ray  spectra  should  also  be  simulated. 
Furthermore,  motion  in  the  chest  wall-to-nipple  direction  (y)  should  be  simulated  in  addition  to  the  x  and  z  directions 
(Figure  2).  Finally,  while  this  work  considers  a  constant  object  velocity,  there  are  instances  in  which  the  velocity  is 
expected  to  be  time-dependent.  For  example,  the  velocity  may  be  sinusoidal  with  time  in  order  to  simulate  the  pulsatile 
motion  of  structures  lying  along  blood  vessels. 
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6.  APPENDIX 

This  appendix  calculates  the  MTF  of  the  focal  spot  in  a  system  with  continuous  x-ray  tube  motion  by  analyzing  rays 
emanating  from  different  focal  spot  positions  during  a  single  projection.  For  the  nth  projection,  rays  are  first  drawn 
between  the  point  (0,  z0)  and  the  two  endpoints  of  the  tube  arc.  Subsequently,  following  Eq.  (19),  rays  are  backprojected 
toward  the  nominal  projection  angle  y/=yjn-  The  length  bn  between  backprojected  rays  at  the  height  z  =  z0  in  the 
reconstruction  yields  a  rectangle  function  corresponding  to  the  PSF  of  tube  motion.  The  MTF  of  focal  spot  motion  is 
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Figure  5.  The  point  spread  function  (PSF)  of  focal  spot  motion  is  found  by  ray  tracing. 

then  calculated  from  the  Fourier  transform  of  the  PSF.  To  implement  this  approach,  it  is  first  necessary  to  determine  the 
lines  £\(y/n  +  TV2)  and  £\(i//n  -  TV2)  between  the  point  (0,  z0)  and  the  endpoints  of  the  x-ray  tube  arc  at  ys  =  y/n  +  TV2  and 

yj=\IJn-  TV2.  Since  the  focal  spot  coordinates  are  given  by  x  =  -hsm(yjn  ±  lF/2)  and  z  =  /zcos(^„  ±  ^72),  the  rays  through 
the  point  (0,  z0)  lie  along  the  lines 

4 (yn  +  'P  /  2)  =  {(x, z) :  z  =  - [cot(^„  +  ¥  / 2)  -  (z0  /  h) csc(K  ±'F/2)]x  +  z0} .  (Al) 

The  lines  £\(y/n  +  vP/2)  and  £\(y/n  -  'i'/l)  strike  the  detector  at  the  points  V\(y/n  +  vF/2)  and  V\(y/n  -  TV2),  respectively 

7>(^±'F/2)=£D[(^±'F/2)/g]n/;(^+'F/2),  4[(^„  ±'P/2)/g]  =  {(x,z):z  =  xtan[(^„  ±Y/2)/g]}  (A2) 

where  £v(y/n  +  TV2)  and  £v(y/n  -  *F/2)  are  lines  along  the  length  of  the  detector  at  the  x-ray  tube  angles  y/=  y/n  +  TV2  and 
yj=yjn-  W2.  Combining  Eq.  (Al)  and  (A2)  gives 

Vl(ysn±'i'/2)  =  (x,z):x  =  zoZ*  tan(^  ±  V / 2),  z  =  zoXt  tan [(</„  ±  'P / 2) / g] tan(<//„  ±vF/2) ,  (A3) 

where 


^  =  [l  +  tan  \{if/n  ±  /  2)  /  g]  tan(^/ 7  ±  T*  /  2)  -  (z0  /  /z)  sec ±  'F  /  2) J 


(A4) 


Signal  at  the  detector  positions  V\(i//n  +  'F/2)  and  -  *F/2)  is  in  turn  backprojected  to  the  focal  spot  at  the  x-ray  tube 
angle  yj  =  y/m  forming  the  lines  £2{y/n  +  vF/2)  and  £2(^«  -  vF/2). 


4(^±Y/2)  =  Ux,z):i 


-  \ji  cos  -  zaxt tan  [(¥„  ±  ¥  /  2)  /  g]  tan(^„  +  ¥  /  2)]  [x  +  h  sin  y/n  ] 


h  sin  y/n  +  z0^  tan(^;7  ±  T1  /  2) 


+  /z  cos  y/n 


(A5) 


At  the  reconstruction  depth  z  =  z0,  the  lines  £2(y/n  +  ¥72)  and  £2(y/n  ~  lF/2)  intercept  the  two  points  V2(y/n  +  *F/2)  and 


7M^-TV2). 
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V2(ysn±'i'/2)  =  £1(V,n±'¥/2)n{(x,z):z  =  z0} 

[A sin y/n  +  z0/jA  tanQ//„  ±  T/2)][/jCosi//„ 

*  cos  Yn  -  ZoZ*  tan  [('//„  + 'tv  2)/ g]  tan(i//„  +  'P/2) 


(A6) 

(A7) 


A  rectangle  function  of  length  bn  is  created  between  +  lP/2)  and  ~  vF/2),  thereby  forming  the  effective  PSF  of 

the  focal  spot  at  the  reconstruction  depth  z  =  z0  (Figure  5).  Using  MATLAB,  one  can  show  that  bn  does  not  vary 
significantly  with  projection  number  n  for  the  Selenia  Dimensions  system.  Consequently,  the  special  case  n  =  0  is  a 
useful  approximation  for  the  effective  width  of  the  PSF  of  focal  spot  motion  in  all  projections. 


2  zn 


f  1  -  — 1  tan 

|  rAy/ 

l  h)  [ 

U-rJ 

2 

rAy/ 

U-rJ 

2 

(A8) 


Since  TV: 2  is  well  under  1°  for  typical  acquisition  geometries,  one  can  use  the  approximation  tanfT'  /(2g)\  tan(*F  /  2)  1 

to  derive  this  result.  With  the  PSF  of  focal  spot  motion  given  by  rect(x/Z?0),  the  MTF  of  focal  spot  motion  is  thus 
sinc(Z?o/).  The  zeros  of  this  MTF  are  integer  multiples  of  1  /b0.  This  formula  perfectly  calculates  the  zeros  of  the 
modulation  of  the  reconstruction  in  Figure  3(c)  for  a  system  with  continuous  tube  motion  and  no  patient  motion. 
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